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I.  TECHNICAL  DISCUSSION 


1.  Introduction 

During  the  performance  period  of  the  contract 
reflected  by  this  final  report  (26  June  1979  to  26  August 
1980)  SAI  personnel  addressed  two  broad  areas  of  importance 
in  understanding  radiation  from  high  temperature  plasmas. 

The  first  area  of  investigation  is  the  detailed 
atomic  physics  and  radiation  transport  which  determines  the 
radiating  efficiency  of  a  plasma.  To  accurately  compute  or 
predict  such  radiation,  it  is  clear  that  the  fundamental 
rates  of  emission  from  all  the  significant  energy  levels  as 
well  as  the  photons'  probability  of  escape  from  the  plasma 
must  be  known.  During  the  performance  period,  SAI ' s  efforts 
were  concentrated  on:  1)  adding  additional  and  necessary  (as 
shown  by  experimental  spectra)  detail  to  the  energy  level 
models  of  various  radiating  elements  and  mixtures  of  elements 
2)  Obtaining  more  computationally  efficient  techniques  to  cal¬ 
culate  the  probability  that  a  photon  escapes  the  plasma  without 
reabsorption  once  emitted.  The  results  of  these  efforts  are 
detailed  in  the  next  section  and  in  the  appendices. 

The  second,  equally  important  area  of  technical 
activity  which  SAI  was  engaged  in  is  the  development  of  hydro¬ 
dynamics  models  which  realistically  take  into  account  the  feed¬ 
back  between  plasma  motion  and  plasma  radiation.  When  a  plasma 
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is  compressed,  the  resulting  heating  obviously  increases  the 
plasma's  resistance  to  further  compression.  However,  radiative 
cooling  processes  -  via  x-rays  for  a  high  temperature  plasma  - 
will  mitigate  this  effect.  Just  how  strong  the  mitigating 
effect  is  will  determine  in  large  measure  the  evolution  of  the 
plasma  and  thus  its  radiative  response  to  higher  power  generation. 
Progress  made  by  SAI  in  understanding  this  crucial  problem  is 
detailed  in  Section  3. 


2.  Ionization,  Excitation,  and  Radiation  Transport 
Modeling 

a)  Spectroscopy  and  Atomic  Physics 

The  excitation,  ionization,  and  radiation  physics  of 

a  dense,  high  temperature  plasma  is  of  importance  for  two 

fundamental  reasons.  First,  the  emitted  radiation  alters  the 

energetics  and  thus  influences  the  evolution  of  the  plasma. 

Second,  ratios  of  spectroscopically  detected  lines  are  a 

valuable  tool  for  diagnosing  plasma  temperature  and  density. 

Appendix  A,  "Self-absorption  of  heliumlike  satellite  lines  in 

high  density  fusion  plasmas" ^  details  studies  of  the  density 

dependence  of  argon,  neon  and  aluminum  satellite  lines  in  two- 

electron  (heliumlike)  ions  taking  into  account  the  temperature 

( 2 ) 

and  opacity  variation.  Appendix  B  ,  "Line  emission  from  hot, 
dense  aluminum  plasmas"  focuses  on  the  behavior  of  K-shell  x-ray 
lines  in  aluminum  as  diagnostic  indicators  of  the  plasma 
temperature  and  density. 

In  many  cases,  most  of  the  radiative  output  of  high 
temperature  aluminum  plasmas  appears  as  softer  x-rays  from  the 
L-shell.  Appendix  C,  "Soft  x-ray  and  XUV  radiation  from  high 
density  aluminum  plasmas"  describes  and  quantifies  the  radiative 
output  of  aluminum  cylinders  (imploded  wires)  as  a  function  of 
electron  temperature,  ion  density,  and  opacity  (related  to 
radius)  over  a  wide  energy  range,  for  the  free-bound  continuum 
as  well  as  lines.  A  milticell  escape  probability  model  was 


used  to  transport  both  lines  and  the  continuum.  This  appendix 
contains  numerous  predictions  of  experimental  spectra  of 
cylindrically  symmetric  imploded  wire  arrays. 

One  result  of  SAI's  effort  during  the  performance 
period  was  a  prediction  that  boron,  radiating  optimally,  is 
not  as  efficient  a  radiator  as  aluminum  radiating  optimally. 
Appendix  D  contains  the  text  of  a  letter  to  the  spectroscopy 
group  at  Maxwell  Laboratories,  Inc.  detailing  this  calculation 
and  its  implications  for  experimenters. 

b)  Radiation  Transport  Modeling 

The  work  of  the  performance  period  was  aimed  at  the 
achievement  of  primary  objectives  of  direct  relevance  to  wire 
and  Baccarat  goals.  First,  to  develop  a  quantitative  analysis 
of  the  effects  of  plasma  motion  on  escape  of  line  radiation 
from  exploded-wire  plasmas.  Second,  to  develop  a  capability  to 
treat  plasma  mixtures,  and  thus  provide  analyses  of  the  effect 
of  one  plasma's  radiation  on  another  plasma.  The  final  objec¬ 
tive  was  to  continue  the  synthesis  of  previous  work  on  the 
exact  transfer  equation  on  the  one  hand  and  on  probabilistic 
photon  escape  on  the  other.  The  purpose  of  this  third  effort 
was  to  enable  exact  calculations  to  be  performed  economically 
on  the  myraids  of  spectral  lines  which  are  observed  to  emanate 
from  exploded-wire  plasmas.  This  would  enable  more  accurate 
and  physically  realistic  simulations  of  how  closely  line 
radiation  may  be  expected  to  approach  black-body  values. 
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These  objectives  were  accomplished,  and  the  separate 
sections  below  set  out  in  detail  the  methods  used  and  results 
obtained . 

i.  Effects  of  Plasma  Motion  on  Yield  and  Diagnostics 

Since  plasmas  in  general  can  be  characterized  by  both 
expanding  and  compressing  regions,  it  was  essential  that 
algorithms  developed  to  treat  the  effects  of  motion  on  radiative 
output  be  capable  of  solving  the  transfer  equation  for  such  non¬ 
monotonic  velocity  fields.  The  result  of  the  effort  was  that 
a  computer  code  containing  such  an  algorithm  was  written  and 

successfully  benchmarked  against  previously  published  line 

(3) 

profile  results  of  Kunasz  and  Hummer  .  The  set  of  subroutines 
which  compute  the  motional  broadening  effects  have  been  inter¬ 
faced  with  our  existing  large  multi-ionization  radiative 
transfer  model. 

The  fundamental  physical  effect  of  line  broadening  by 
plasma  mass  motions  is  that  the  line  photons  are  more  easily 
able  to  escape  the  plasma  because  they  are  shifted  away  from 
local  line  center,  where  the  absorption  is  highest.  Thus,  an 
optically  thick  plasma  is  turned  in  the  direction  of  being  an 
optically  thin  one.  This  affects  the  emitted  power  of  an 
exploded-wire  plasma  in  a  positive  direction  -  increasing  the 
power.  Some  photons  which  would  be  bottled  up  near  line  center 
and  quenched  in  a  static  plasma  will  instead  escape  and  thus 
increase  the  yield. 
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To  gauge  the  effect  of  plasma  implosion  on  line 
radiation  we  select  one  of  the  strongest  single  lines  emitted 
xn  typical  exploded-wire  plasmas,  the  helium-like  Is  -  ls2p  P 
line.  An  inward-directed  velocity  profile  V(Z)  =  VQZ/Zo  has 
been  chosen,  where  VQ  =  1.3  x  107  cm/sec,  ZQ  =  plasma  half¬ 
width  (750  microns),  and  the  temperature  and  density  are  typical 
of  a  wire  plasma  (600  eV,  1  x  lO1^  cm-3)  at  600  eV,  1.3  x  107 
is  exactly  twice  the  thermal  velocity.  Of  course,  as  VQ  *  ®, 
the  plasma  would  become  thin.  Much  smaller  velocities  would 
produce  essentially  the  results  for  a  static  plasma. 

In  Fig.  1  are  plotted  the  profiles  for  the  static 
case  and  the  imploding  case  described  above.  Line  power  out¬ 
put  (i.e.,  the  integral  under  the  profile)  does  not  change 
because  the  plasma  in  this  case  is  effectively  thin,  that  is, 
the  photons  generally  escape  after  many  scatterings,  at  this 
density.  The  profile,  however,  charges  quite  dramatically. 

The  static  plasma  produces  a  self-absorbed  line  with  symmetrical 
blue  and  red  peaks,  as  expected.  The  imploding  plasma  line 
profile's  blue  peak  far  exceeds  its  red  in  intensity.  This  is 
due  to  the  fact  that,  in  a  contracting  plasma,  every  emitted 
photon  is  blueshifted  relative  to  the  rest  of  the  plasma.  Thus 
photons  emitted  on  the  blue  side  of  the  line  center  are  'seen' 
as  even  bluer  by  the  rest  of  the  plasma  and  thus  escape  easily. 
Photons  emitted  on  the  red  side,  however,  are  shifted  blueward 


to  line  center  where  they  are  largely  absorbed  due  to  the 
high  opacity  at  line  center.  Hence  the  domination  of  the 
blue  peak  for  the  imploding  case. 

20 

Results  at  a  total  ion  density  of  10  have  been 

obtained  using  the  same  total  mass  (75u  plasma)  and  temperature 

(.600  eV)  .  The  line  profile  phenomenology  is  the  same  as  that 
19  20 

described  at  10  .  However,  at  10  the  line  energetics  are 

,  2  1 
significantly  altered.  The  A1  XII  Is  -  ls2p  P  line  remains 

unchanged  in  output  by  motional  effects.  Since  the  decay 

probability  is  so  high  (2.8  x  103-3  sec-'1')  ,  the  probability 

that  excitation  of  the  2p  level  is  not  followed  by  spontaneous 

decay  to  the  ground  state  is  extremely  low.  Hence,  most  photons 

in  this  line  continue  to  excape  even  in  the  static  case  and 

motions  have  no  effect.  However,  there  are  effects  on  other 

lines.  The  following  table  gives  the  ratio  of  various  lines' 

2 

output  to  that  of  the  Is  -  ls2p,  for  the  static  and  imploding 
2  0 

cases  at  =®  10  ,  Te  =  600  eV,  and  plasma  half-width  of  75u- 


Line 

Static 

Imploding 

Is2  -  ls2p1P(Al  XII) 

1.000 

1.000 

Is2  -  ls3p1P(Al  XII) 

0.093 

0.141 

Is2  -  ls4p1P(Al  XII) 

0.037 

0.053 

ls22p  -  ls2p2  2D(A1  XI) 

0.046 

0.057 

Is2  -  ls2p3P (Al  XII) 

0.071 

0.099 

Is  -  2p  (Al  XIII) 

1.226 

2.102 

Is  -  3p  (Al  XIII) 

0.156 

0.237 

Is  -  4p  (Al  XIII) 

0.175 

0.210 
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Clearly,  every  line  but  the  Is2  -  IsZp1?  of  Al  XII 
has  increased  in  output,  with  the  very  strong  ls-2p  hydrogen¬ 
like  resonance  line  nearly  doubling  because  motional  broadening 
renders  the  line  thinner  than  it  would  be  in  the  static  case. 

ii .  Analysis  of  Plasma  Mixtures 

Developing  an  understanding  of  how  radiation  from  one 
plasma  affects  another  is  important  in  analyzing  the  effects 
of  dissimilar  plasmas  upon  each  other  and  in  developing  a 
diagnostic  capability  for  wire  or  puff  shots  when  two  or  more 
different  plasma  elements  are  actually  present.  Figure  2, 
taken  from  NRL  Memorandum  Report  3934,  by  P.  Burkhalter  and 
J.  Davis,  shows  a  mixed  silicon-aluminum  spectrum  from  PITHON 
shot  PI  6.  Line  ratios  from  both  aluminum  and  silicon  can  be 
used  to  firm  up  the  diagnostic  information  obtained  from  such 
a  spectrum. 

A  part  of  the  contract's  work,  a  capability  to  analyze 
silicon-aluminum  mixtures  using  our  detailed  radiation  transport- 
multiionization  dynamics  code  was  established.  The  model  was 
built  in  the  following  manner.  First,  a  rate  table  for  silicon 
transitions  precisely  analogous  to  the  existing  aluminum  table 
was  set  up.  Eight  prominent  spectral  lines  (the  ones  listed  in 
the  table  contained  in  the  preceding  section)  are  fully  self- 
consistently  transported.  The  model  includes,  of  course,  the 
effects  of  plasma  motion  as  described  in  Section  II.  Second, 
the  emitted  spectrum,  angle-by-angle,  can  be  stored  on  the 
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TI/ASC's  disk  files  for  either  silicon  or  aluminum.  Third, 
when  the  silicon  or  aluminum  model  is  run,  any  stored  spectrum 
on  the. disk  from  either  element  can  be  read  and  self-consistently 
transported  through  the  element  being  run  using  the  disk  spectrum 
as  a  boundary  condition.  The  total  spectrum  from  both  elements 
is  outputted,  and  thus  the  effect  of  mixture  can  be  gauged. 

A  dramatic  example  of  how  radiation  from  a  silicon  plasma 
affects  the  state  and  emitted  spectrum  of  an  aluminum  plasma  can 
be  seen  by  comparing  the  next  two  figures.  Figure  3  is  the  spec¬ 
trum  emitted  by  an  aluminum  plasma  at  a  temperature  of  100  eV, 

19 

total  ion  density  of  5  x  10  ,  and  half-width  of  5.56y.  The  lines 

are  labelled,  and  the  slanted  line  above  the  spectrum  is  the 

emission  from  a  100  eV  blackbody.  The  peaks  of  all  the  lines 

are  comfortably  under  the  Planck  curve.  By  contrast,  Fig.  4 

illustrates  the  spectrum  emitted  by  the  same  aluminum  plasma 

when  it  is  irradiated  by  a  silicon  plasma  surrounding  it  in  a 

'sandwich'  configuration.  The  silicon  plasma  is  at  a  tempera- 

19 

ture  of  400  eV  and  total  ion  density  of  8  x  10  .  There  is 

assumed  a  1500y  layer  of  silicon  on  each  side  of  the  aluminum. 

The  spectrum  emitted  by  the  embedded  aluminum  (with  the  silicon 
lines  removed  for  clarity)  is  radically  different  from  that  of 
the  bare  aluminum.  In  this  case  we  have  assumed  that  the 
silicon  exterior  of  the  plasma  sandwich  is  imploding  toward 
the  aluminum  at  a  velocity  of  3  x  10^  cm/sec.  The  aluminum 
lines  now  are  far  above  the  100  eV  blackbody  curve  in  intensity, 
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even  though  the  kinetic  temperature  of  the  aluminum  plasma 

is  100  eV.  This  effect  is  of  course  due  to  the  fact  that  the 

silicon  radiation,  both  lines  and  continuum,  photoionizes  and 

photoexcites  bound  levels  of  the  aluminum  plasma.  This  raises 

the  state  of  excitation  of  the  aluminum  plasma  far  above  what 

is  obtained  from  a  "bare"  100  eV  aluminum  plasma  -  in  fact  even 

far  above  that  of  an  LTE  100  eV  aluminum  plasma.  The  ratios  of 

lines  are  quite  significantly  changed  also,  with  the  satellite 

line  actually  seen  in  absorption  in  the  pumped  plasma.  The 

2 

pumped  plasma  also  exhibits  an  ultra-strong  A1  XII  Is  -  ls3p 
o 

line  (at  6.635A).  This  effect  is  due  to  the  fact  that  the 
2  o 

Si  XIII  Is  -  ls2p  line  lies  at  6.650  A,  and  therefore  photons 

emitted  in  this  line  by  the  outer  silicon  region  are  inter- 

2 

cepted  in  the  Al  XII  Is  -  ls3p  transition,  greatly  augmenting 
its  output.  The  lesson  is  clear  -  when  two  plasmas  are  inter¬ 
acting,  the  photons  from  one  can  drastically  affect  the 
excitation  state  and  thus  the  diagnostic  interpretation  of 

radiation  emitted  by  the  other  plasma. 

2 

The  pumping  of  the  Al  XII  Is  -  ls3p  line  by  the 
2 

Si  XIII  Is  -  ls2p  line  referred  to  above  is  also  useful  in 
principle  as  a  diagnostic  of  exploding  -  wire  plasma  mixtures. 
Under  certain  conditions,  the  Al  XII  3-2  states  may  become 

o 

inverted  in  population,  resulting  in  the  lasing  of  the  44  A 
lines  coupling  these  states.  The  establishment  of  a  population 
inversion,  and  hence  the  strength  of  this  line,  is  sensitively 


dependent  on  properties  of  the  plasma  mixture.  As  an  example 
of  this.  Figs.  5  and  6  display  the  gain  coefficient  of  the 
A1  XII  ls2p^P  -  ls3d^D  line  {.the  strongest  lasing  line)  as  a 
function  of  aluminum  plasma  density  and  size.  Such  variations 
with  size  and  density  reflect  the  extent  to  which  Si  XIII 
pumping  photons  are  able  to  penetrate  the  aluminum  region,  as 
well  as  the  collisional  quenching  of  the  gain  by  electron 
impact  processes  occurring  in  the  aluminum  region.  Our  studies 
of  these  pumping  effects  also  indicate  considerable  sensitivity 
of  the  gain  to  aluminum  plasma  temperature.  The  silicon  plasma's 
radiation  tends  to  over-ionize  the  aluminum  plasma  and  thus  little 
or  no  gain  is  observed  at  temperatures  in  the  1  Kev  range.  It  is 
to  be  emphasized  that  the  degree  of  gain  present  in  the  3-2  tran¬ 
sitions  will  be  sensitively  reflected  in  the  experimentally 

o 

observed  strength  of  the  44  A  lines.  In  turn,  this  experimental 
l  quantity  reflects  the  state  of  the  component  silicon  and  aluminum 

l  plasmas.  Much  remain?*,  to  be  understood  on  the  extent  to  which 

such  pumping  effects  are  usable  as  plasma  mixture  diagnostics, 
iii.  Probabilistic  Radiation  Transport 

As  reported  previously,  SAI  has  developed  an  efficient 
escape  probability  technique  for  transporting  line  radiation.  In 
the  performance  period  of  the  present  contract,  the  technique  was 
used  to  analyze  effects  of  optically  thick  lines  arising  from 


excited  states  in  an  aluminum  plasma.  A  copy  of  the  completely 
updated  report  on  this  method, ^  "Direct  Solution  of  the 
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Equation  of  Transfer  Using  Frequency  -  and  Angle  -  Averaged 
Photon  Escape  Probabilities,  with.  Application  to  a  Multistage, 
Multilevel  Aluminum  Plasma",  is  included  as  Appendix  E. 

iv.  Effect  of  Experimental  Resolution  on  Spectra 
One  of  the  best  ways  to  diagnose  the  properties  of 
wire-like  and  imploding  pellet  plasmas  is  the  analysis  of 
various  spectral  features.  Every  experimentally  measured 
spectrum  is  affected,  to  a  greater  or  lesser  degree,  by  instru¬ 
mental  and  source  broadening  which  is  generally  assumed  to  be 
Gaussian  in  character.  To  provide  a  quantitative  and  detailed 
capability  to  gauge  the  effects  of  experimental  broadening,  a 
program  was  established  to  catalog  all  the  spectra  calculated 
theoretically  and  convolute  them  with  Gaussian  broadening 
functions.  The  convoluted  spectrum  is  then  suitable  for 
direct  comparison  with  experiment.  Our  results  graphically 
illustrate  the  evolution  of  a  "pure  theoretical"  spectrum 
into  an  observed  spectrum  as  the  experimental  broadening 
increases.  This  evolution  is  seen  in  the  next  section  for 
the  case  of  an  opacity  and  velocity-broadened  aluminum  spectrum. 
This  spectral  convolution  capability  will  provide  a  much  more 
direct  connection  of  our  theoretical  models  with  experiment 
and  thus  aid  in  the  determination  of  the  basic  properties  of 
both  wire  and  pellet  plasmas. 
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V. 


Analysis  of  Motional  Broadening  in  an  Aluminum 
Wire  Plasma  Spectrum 

The  principal  resonance  line  of  helium-like  Al  XII, 

2  1 

Is  -  ls2p  P,  as  observed  during  Maxwell  Laboratories'  shot 

159,  exhibited  an  interesting  feature  on  the  red  side  of  the 

line  profile,  as  seen  in  Fig.  7.  As  discussed  above,  such  a 

red  wing  can  be  produced  by  plasma  implosion  when  the  line 

is  optically  thick  enough  to  be  self  reversed.  In  the  case 

of  implosion,  the  blue  wing  exceeds  the  red  in  intensity. 

Using  both  a  simplified  probability-of-escape  model  and  a 

detailed  ray-trace  aluminum  plasma  model,  a  radius  of  250y, 

18  -3 

temperature  of  500  ev,  and  total  ion  density  of  3  x  10  cm 
were  obtained  as  characteristic  plasma  parameters  from  the 
various  spectral  line  ratios.  The  detailed  model  agreed 
very  well  with  the  simpler  code  in  this  case,  since  the 
plasma  is  "effectively  thin".  The  table  gives  the  measured 
and  modeled  line  ratios  for  the  above  parameters. 


r  =  250y, 

T  =  500  ev,  N.  = 
e  l 

18  - 

8  x  10  cm  - 

Ratio 

(Spherical) 

Measured 

Modeled 

Ljls2  -  ls2ph 

0.14 

0.12 

Lg|ls2  -  ls2p1P 

0.04 

0.031 

2  1,2 
is  -  lsap^Pl is"1  - 

ls2p1P 

0.015 

0.011 

Is2  -  ls4p1p|ls2  - 

1S2P1? 

0.007 

0.007 
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Given  these  plasma  properties,  a  series  of  computations 
with  different  implosion  velocities  were  performed  to  attempt 
to  diagnose  the  characteristic  velocity  in  shot  159  from  the 
resonance  line  profile.  An  implosion  velocity  of  4  x  lO^cmsec- 
at  the  outer  radius  (linearly  ramping  from  r  =  0)  was  found 
to  reproduce  quite  well  the  observed  profile.  This  velocity, 
which  is  obviously  quite  important  to  know  in  order  to  under¬ 
stand  the  conversion  of  plasma  kinetic  energy  to  x-rays,  agrees 
quite  well  with  theoretical  hydrodynamic  calculations.  The 
major  uncertainty  in  the  analysis  is  the  amount  of  experi¬ 
mental  broadening  in  the  spectrum  in  Fig.  7.  The  amount  of 

o 

experimental  broadening,  AA  ,  in  A  (full-width-half-maximum) 
is  given  by 


AA 


e 


/TAA)^  -  ( A  A  >  ( 


where  ( A  A )  and  ( A  A )  .  are  the  measured  and  intrinsic  line 
m  i 

FWHM's.  The  measured  FWHM  of  the  A1  XII  intercombination 

o 

line  is  6.2  +  1.5  m  A.  The  intrinsic  line  width,  obtained 

o 

from  a  radiative  transfer  calculation  is  approximately  4mA. 

o 

Clearly,  the  widths  are  consistent  with  very  little  (^1  m  A) 

o 

or  quite  a  bit  (v7  m  A)  of  broadening.  The  effect  of  broaden¬ 
ing  on  the  observed  spectrum  is  dramatically  illustrated  in 

Figs.  8-11  where  the  spectrum  with  broadening  of  0.,  1., 

o 

1.59,  and  2.73  m  A,  respectively,  is  shown.  The  "washing- 
out"  of  the  motional-broadening  feature  is  clearly  seen. 
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J 


Figure  11 


Interestingly,  the  most  dramatic  agreement  of  the  line  profile 
is  seen  for  no  assumed  experimental  broadening.  The  fact  that 
this  feature  is  observed  would  appear  to  indicate  that  the 
source  and  instrumental  broadening  were  indeed  quite  low  in 
this  case.  It  should  also  be  pointed  out  that  the  Ly  a 
modeled  width  is  smaller  than  the  experimental  width  because 
the  fine  structure  is  not  modeled.  However,  the  integrated 
intensity  of  this  line  agrees  quite  well  with  the  experiment, 
vi.  Properties  of  an  Imploded  Neon  Plasma 
One  technique  which  has  been  successfully  used  to 
obtain  information  characterizing  imploded  pellets  is  to 
utilize  an  inert  filler  gas  such  as  argon  or  neon,  and  analyze 
the  filler  spectrum.  To  aid  in  such  analysis,  a  multi-stage, 
multi-level  neon  model  was  constructed  for  use  in  our  large 
radiative  transfer  code.  Self-consistent  Stark  profiles 
supplied  by  P.C.  Kepple  were  used  as  the  line  opacity  contour. 

To  benchmark  this  study,  we  have  reproduced  remarkably  well 
with  the  model  an  experimental  spectrum  previously  obtained 
at  Rochester.  We  obtain  as  a  fit  for  this  spectrum  a  density 
of  4.5  x  1021  ions  cm-3  and  temperature  of  385  ev.  By  contrast, 
the  original  Rochester  analysis  yielded  Tg  =  300  ev  and 
Ni  =  7. 8  x  1021cm-3.  Figure  12  displays  the  experimental 
spectrum  at  the  top,  followed  by  the  "best  fit"  theoretical 
spectrum.  The  lowest  spectrum  on  the  graph  is  that  which 
would  be  emitted  by  a  plasma  characterized  by  the  conditions 
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8.6  ATM. 


ENERGY  (keV) 


of  the  Rochester  analysis.  The  increase  in  temperature  to 

385  ev  is  necessary  to  remove  the  prominent  helium  continuum 

feature  which  would  occur  at  300  ev  and  is  not  observed. 

21  -3 

The  reduction  in  density  (from  7.8  to  4.5  x  10  cm  )  is 

21  -3 

necessitated  by  the  fact  that  at  7.8  x  10  cm  ,  opacity 
broadening  renders  the  predicted  line  widths  greater  than 
those  experimentally  observed.  This  benchmark  analysis 
strikingly  illustrates  the  great  value  of  analyzing  all  the 
features  in  a  given  spectrum  to  obtain  precise  plasma  parameters 


3.  Radiation  Hydrodynamics  Modeling 

a)  Cylindrical  Geometry 

It  is  generally  accepted  that  cylindrical  geometry 
is  the  most  realistic  simple  geometry  for  modeling  collapsing 
z-pinches.  Significant  effort  has  been  performed  on  the 
establishment  of  a  collapsing  z-pinch  hydrodynamics  code 
WHYRAD,  initialized  by  wires,  gas  puffs,  or  annular  foils, 
employing  a  variety  of  ionization  physics-radiation  input 
subroutines.  A  wide  variety  of  initial  radial  distributions 
of  particle  number  density  and  temperature,  timewise  variation 
of  machine  discharge  voltage,  and  anomalous  resistivity 
functions  may  be  employed;  in  addition,  operational  flexibility 
due  to  general  re-start  and  subroutine  vectorization  as  well 
as  spectral  and  radial  distribution  plotting  routines  is  also 
available. 

i.  One-Dimensional  z-Pinch  Model 
Analysis  of  the  dynamics  of  an  imploding  z-pinch 
plasma  depends  on  magnetohydrodynamic  coupling  of  the  time- 
dependent  equations  of  mass,  momentum,  electron  and  total 
specie  energy,  magnetic  flux  and  the  thermal  and  caloric 
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equations  of  state,  as  well  as  the  specification  of  appropriate 
transport  processes.  Solution  of  these  equations  is  also  de¬ 
pendent  on  the  coupling  into  the  driving  discharge  circuit, 
which  in  turn  is  altered  by  the  changing  plasma  conditions. 

The  appropriate  equations  are  assumed  to  be  one  dimensional, 
and  cy lindncally  symmetric.  The  energy  density  of  the  plasma 
involves  an  adiabatic  exponent  which  contains  the  potential 
ionization  energy  that  is  calculated  from  a  selected  ionization- 
radiation  package.  Thermal  flux  coefficients  are  determined 
from  the  expressions  listed  in  Braginski ;  ' ® ^  thermally  per¬ 
fect  plasma  species  and  charge  neutrality  are  also  assumed. 

Magnetic  field  diffusion  within  the  z-pinch  involves 
n ,  a  cross-field  resistivity  having  a  limiting  classical 
vaiue  of 

c 

n,  =  1.15  x  10'14  zT  inA  (T  )"3'2  i. sec )  :1) 

1  I  e 


where  Zj  is  the  ionic  charge,  in  A  is  the  Coulomb  logarithm 
and  T  ,  the  electron  temperature,  is  in  eV.  Experimental 
results  indicate  that  non-classical  anomalous  turbulent 
values  of  some  of  the  transport  processes,  in  particular 
the  resistivity,  are  required  for  suitable  correlation. 

Plasma  turbulence  would  be  strongly  influential  in  the  outer 
lower  density  coronal  region  of  the  discharge  and  for  tens 
of  nanoseconds  and  above  time  scales.  As  a  first  model  it 
is  assumed  that  only  the  resistivity  is  affected  by  anomalous 
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turbulence  levels,  the  other  transport  pr6perti.es  retaining 
their  classical  magnitudes.  The  empirical  quadratic  relation 
of  the  form 


an 


=  n 


ci 


Jl. 


0  +  (g ( t)  -  1.0)  (r/R) 


(2) 


is  chosen,  where  R  is  the  outer  discharge  radius  and  the 
growth  factor  g  is 


g  (t) 


1.0  +  A  (t/t  .  ) 

pulse 


(3) 


where  A  and  B  are  empirical  coefficients  and  t  .  is  the 

pulse 

total  pulse  time.  The  core  limit  of  n^n  -*■  as  r  -  o 

is  maintained  as  well  as  the  initial  condition  ^an  - 


1  1 


as  t  -  o  for  any  r. 

Solution  of  the  equation  system,  except  for  the  mag¬ 
netic  flux  relation  and  the  thermal  diffusion  is  by  appli¬ 
cation  of  the  explicit  flux  corrected  transport  (FCT) 
i  8  )  f  9  ) 

method^  This  method  effectively  transports  steep  fluid 

gradients  and  does  not  require  the  incorporation  of  artifi¬ 
cial  viscosity.  Initial  condition  distributions  of  the  de¬ 
pendent  variables  are  specified  at  ail  finite  difference 
radial  locations  either  at  t  =  o  or  at  an  appropriately 
chosen  restart  time.  The  magnetic  field  and  thermal  diffusion 
equations  are  solved  by  an  implicit  tndiagonal  .Gaussian 
elimination)  algorithm.  The  radial  finite  difference  repre¬ 
sentation  is  incrementally  uniform  with  zero  gradients  of 
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dependent  variables  assumed  at  the  discharge  centerline  and 
outer  plasma  boundary. 

The  discharge  plasma  within  its  return  current  cage 

is  modeled  as  a  plasma  cylinder  in  contact  at  one  end  with 

a  grounded  perfectly  conducting  metal  return  current  path 

and  at  the  other  end  with  an  electrode  at  voltage  V  ,  ft). 

pi 

Specification  of  Ohm's  law  and  Faraday's  equation  within  the 
discharge  yields  the  radial  electric  field  distribution  in 
terms  of  V  ,  and  the  3Bq/3t  flux  of  the  azimuthal  magnetic 
field . 

The  ionization-radiation  packages  employed  are  ATPHYS, 
a  relatively  simple  ground  state  only  ion  balance  with 
optically  thick  radiation,  and  MCELL,  a  package  which  solves 
for  the  ground  and  excited  state  population  densities.  Line, 
radiation  recombination  and  Bremsstrahlung  radiation  trans¬ 
port  are  considered,  with  the  ionization  balance  determined 

by  a  number  of  distinct  processes  within  the  coilisional- 

.  1  1 ' 

radiation  equilibrium  method.  Radiation  transport  is 

modeled  with  emission  and  absorption  balanced  out  ever  the 
multi-cell  radial  increment  zones. 

To  check  the  feasibility  of  the  WH V RAD -MCELL  linkage, 
a  series  of  optically  thin  aluminum  k-shell  computations 
were  performed.  While  physical  validity  may  be  questioned, 
expecialiy  for  10  eV  initial,  wire  plasma  temperatures, 
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useful  information  regarding  parameter  dependencies,  such 
as  due  to  anomalous  resistivity,  may  be  obtained.  The  first 
series  of  runs  assume  a  representative  model  voltage-time 
distribution  illustrated  in  Figure  13,  having  a  maximum 
value  of  1.5  x  10^V  and  symmetric  within  a  pulse  time  of 
2  x  10  7sec.  The  radial  extent  is  broken  up  into  33  incre¬ 
ments,  and  300  time  steps  are  employed  with  an  initial 
Gaussian  distribution  of  ion  number  density 


Nr  =  Nn 


max 


exp 


{-lr-ro)2/R2} 


(4) 
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where  NT  =  5.0  x  10  cm  ,  r  ,  the  initial  wire  radius 

I ,  max  o 

being  0.75  cm  and  R  =  0.25  cm.  This  combination  of  para¬ 
meters  together  with  a  discharge  length  of  3.0  cm  and  an 

outer  plasma  edge  of  1.32  cm  yields  a  total  mass  of 
-  4 

4.72  x  10  gm.  This  is  of  the  same  order  of  magnitude  as 
experimental  multi-wire  configurations.  Initial  ion  and 
electron  temperatures  are  assumed  to  be  10  eV;  as  the  opti¬ 
cally  thin  radiation  quickly  cools  the  plasma  to  values 
below  i  eV,  which  is  in  a  range  with  poor  representation 
of  continuum  radiation,  the  radiation  power  feedback  from 
MCELL  to  WHYRAD  is  cutoff  for  an  initial  limited  number  of 
time  steps  (30  or  10%  of  the  total).  Such  a  procedure 
seems  to  be  adequate  for  the  various  input  data. 


Figure  14  illustrates  T  at  four  different  times 

0 
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FIGURE  14:  Optically  Thin,  K-Shell  , 

A£  T  (r,t)  Distribution.  B=10  , 

Radiilly  Uniform  n 


with  A  =  99.0  and  B,  the  timewise  resistivity  power  factor 

-4 

in  Eq.  (3),  equalling  10  ,  essentially  the  always  completely 

turbulent  extreme.  The  radial  dependency  of  n  ,  given  by 
Eq.  (2)  ,  is  suppressed  by  letting  r  =  R  so  as  to  emphasize 
the  timewise  influence.  It  is  seen  that  for  later  times 
and  larger  radii  the  electron  temperature  drops  below  the 
1.0  eV  cutoff.  It  is  to  be  noted  that  there  is  a  minimum 
amount  of  radial  compression  through  the  run  period  due  to 
the  lack  of  radiant  absorption.  Figure  15  shows  Tg  at 
essentially  the  same  four  times  with  B  =  1.0,  modelling 
a  linear  increase  between  the  classical  and  100  times 
classical  limits.  The  Tg  distributions  fall  off,  radially, 
somewhat  more  steeply  than  the  previous  case  beyond  the 
maximum  voltage  point  (at  10  7  sec.);  up  to  there  the 
distributions  are  virtually  identical.  This  illustrates 
an  interesting  transient  effect  due  to  the  resistivity 
magnitude;  when  there  is  a  factor  of  4  and  2  in  resistivity 
difference,  respectively,  at  50  and  100  nsec  the  temperature 
distributions  are  very  similar.  However,  at  later  times, 
when  the  resistivities  are  closer  in  magnitude,  the  temper¬ 
atures  are  much  farther  apart.  Accurate  representation 
of  the  resistivity  growth  rate  is  therefore  important. 
Analagous  effects  on  the  radial  distribution  of  axial 
current  and  azimuthal  magnetic  field  should  also  be  noted. 
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The  more  complete  representation,  involving  a 
radial  dependency  (Eq.  (2))  with  a  laminar  core-turbulent 
corona  is  tested  by  correlation  with  a  set  of  Ar  gas  puff 
data,  shot  #1440  at  Physics  International.  The  ionization- 
radiation  option  is  the  ATPHYS  ground  state  model.  Even 
though  good  correlation  of  the  magnitude  of  the  radiant 
emission  as  well  as  any  information  regarding  the  spectral 
line  distribution  cannot  be  made  with  the  simple  radiation 
option,  it  is  informative  to  compute  the  timewise  variation 
of  the  line  and  continuum  radiant  emission  versus  the  ex¬ 
ternal  circuit  tube  voltage  distribution  and  compare  with 
the  data.  Furthermore,  the  initial  ion  radial  number 
density  distribution  furnishes  a  superior  estimate  than 
the  usual  Gaussian. 

Figure  16  shows  the  timewise  voltage  distribution 
employed;  linear  fits  over  26  time  segments  closely  model 
the  experimental  input  out  to  the  400  nsec  cutoff.  Figure 
17  illustrates  the  N^(r)  @  t  =  0  distribution;  values  at 
radial  distances  beyond  1.7  cm  is  a  linear  estimate.  Cathode 
to  anode  discharge  length  is  5.0  cm;  it  is  also  assumed 
that  initial  radially  uniform  Tg  and  T^  magnitudes  of  2.0 
eV  are  suitable.  Figures  18  and  19  illustrate  the  electron 
temperature  and  average  Ar  ion  number  density  at  4  separate 
time  steps  over  a  390  time  step  calculation,  with  the  flow 
field  divided  into  39  radial  increments.  It  is  assumed 
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FIGURE  19:  Radial  Distribution  of  Ar  ION  Number  Density 


that  the  maximum  time  step  over  the  first  170  nsec  is  2.0 
nsec  and  after  this  the  timestep  re-adjusts  to  accomodate 
local  changes  proportional  to  the  net  radial  plasma  velocity, 
plasma  sound  speed,  etc.  Indeed,  time  steps  near  the  region 
of  maximum  beam  collapse  can  get  to  as  low  as  10  ^  sec. 

As  is  seen  from  Figure  18  the  total  zonal  width  de¬ 
creases  somewhat  with  time;  a  maximum  electron  temperature 
of  over  10^  eV  is  achieved  at  ~230  nsec  which  is  close  to 
the  initial  collapse  time.  The  illustrated  temperatures 
are  high  due  to  the  neglect  of  higher  ionization  states 
in  the  radiation-chemistry  subroutine;  nevertheless,  refer¬ 
ence  to  the  Nj  distributions  illustrated  in  Figure  19  show 
that  the  temperature  peaks  are  close  to  the  steep  density 
dropoff  regions.  This  is  representative  of  a  core-corona 
model  with  the  high  density,  somewhat  cooler  region  sur¬ 
rounded  by  a  higher  temperature,  strongly  radiating,  lower 

density  annular  zone.  The  radially  quadratic  distri- 

-4 

bution  (here  with  B,  the  time  factor,  equalling  10  )  is 

seen  to  have  some  physical  validity.  Maximum  N ^  density 

18  —  3 

levels  exceeding  10  cm  correspond  to  values  between 
19  20  -3 

10  and  10  cm  ,  which  seem  to  be  present  m  the  P.I. 
experiment  (private  communication  from  P.  Burkhalter) . 

Figure  20  illustrates  the  simple  Ar  model  radiative 
power  (in  arbitrary  units) ,  versus  time  results  with  separ¬ 
ate  distributions  for  line,  continuum,  and  hv  >1  keV  line 
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radiation.  Computation  these  /a.  .es  initiate  at  tr.e 
one  hundredth  time  step  .  The  timewise  iistr  muti  .r.  .  :  tr.e 
radiative  power  can  be  compared  with  experiments .  r'-rauiat. ..  r. 
emission,  as  shown  in  Fiuure  21.  It  is  seen  that  tr.e 
numerical  results  occur  somewhat  earlier  than  the  experi¬ 
mental  output;  furthermore,  the  .numerical  ustrii  -t.  .ns 
flatten  out  for  over  50  nsec,  which  is  not  fm_.ru  expe:  .  r  to¬ 
tally.  In  reality,  plasma  dynamic  instabilities  may  serve 
to  break  up  and  therefore  extinguish  the  radiatmc  pmer. 
at  times  much  beyond  the  maximum  collapse  point, •  why  PAL  *s 
hydrodynamically  insufficiently  sophisticated  to  acco-nt 
for  this. 


The  latest  employment  of  WHY  RAD  is  the  ccuplmu  to 
a  full  L-shell  representation  of  A<  within  the  MCELL  routines 
140  ground  and  excited  electronic  states  with  12  possible 
ions  are  considered.  Line  radiation  transport  between  all 
the  ionic  states,  balanced  across  the  radial  cells  as  is 
seen  from  the  above  solutions  using  the  mere  approximate 
ATPHYS  package  employed  for  Ar,  about  40  spatial  increments 
are  necessary  to  adequately  represent  the  collapsed  dis¬ 
charge  steep  radial  gradients:  requires  a  formidable  amount 
of  computation.  Asanexample,  inclusion  of  full  line  and 
continuum  transport  costs  about  310  per  time  step  or.  tr.e 
Texas  Instruments  A SC  vectorized  computer,  while  an  entire 
500  time  step  computation  using  the  ATPHYS  -round  state 
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radiation-ionization  package  costs  around  340.  We  feel 
that,  at  present,  it  is  mere  cost  effective  to  consider 
the  appreciacle  simplifications  of  local  radiation  transport, 
i.e.,  line  radiation  is  either  locally  absorbed  within  the 
same  emission  cell  or  completely  escapes  the  plasma,  or 
employ  a  simpler  restricted  excited  state  version  of  the 
MCELL  pacxaue.  Effort  has  therefore  been  expended  on  a 
somewhat  limited  package  that  omits  a  number  of  the  excited 
1-shell  states  but  adequately  represents  the  bound-bcund 
transport  at  the  c-pinch  plasma  densities  and  at  temperatures 
above  3  00  e'.'.  It  is  assumed  that  continuum  radiation 
dominates  at  lower  temperatures.  Similar  considerations 
have  beer,  initiate.!  with  anon  sc  that  :as  puff  solutions 
can  be  made  with  a  form  of  the  MCELL  package . 

ii.  Jcnc »  us i on 

The  core i nation  of  numerically  coupled  z-p inch  and 
mu.ti-ce*.  radiation  codes  together  with  close  correlation 
witr.  availaole  experimental  data  is  we  1  *  underway .  The 
wc  r.<  i  s  a.sc  c.:selv  integrated  with  other  hvdrodvnamic- 


crr.bir.ee  codes  so  that  si  mi:  leant  ccrrcr.alitv 
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a  self-consistent  spherical  numerical  radiation-hydrodynamic 
model  for  use  in  predicting  emission  from  a  moving,  high 
density  plasma.  The  advances  in  the  radiation-hydrodynamic 
model  include  the  capability  to  treat  opacity  effects,  the 
creation  of  an  argon  plus  DT  model,  the  extension  of  the 
aluminum  model  to  lower  ionization  stages,  and  a  restart 
and  graphics  ability. 

The  radiation  hydrodynamics  model  solves  the  time 
dependent  hydrodynamic  equations  for  the  conservation  of 
matter,  momentum,  and  energy  for  a  single  temperature, 
single  fluid  plasma  using  a  flux-corrected  transport  (FCT ) 
algorithm,  with  electron  thermal  conduction  effects  being 
taken  into  account  through  a  time-splitting  of  the  energy 
equation.  The  energy  equation  also  includes  the  effects 
of  heating-cooling  due  to  radiation  transfer.  The  plasma 
is  assumed  to  be  in  ionization  equilibrium,  with  the  ground 
and  excited  levels  being  calculated  with  a  collisional- 
radiative  equilibrium  (CRE)  model.  During  the  performance 
period  the  ionization  model  for  aluminum  was  extended  from 
its  original  level  structure  containing  excited  states  only 
for  hydrogen-like,  helium-like  and  lithium-like  ions  to 
include  exicted  states,  and  hence  line  emission,  from  all 
K-shell  and  L-shell  ions.  In  addition,  an  ionization  model 
for  argon,  including  excited  states  and  line  emission  for 
hydrogen-like,  helium-like,  and  lithium-like  ions  was  developed. 
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The  ability  to  take  opacity  effects  into  account  has 
been  incorporated  into  the  model.  Radiation  transfer 
of  both  line  and  continuum  photons  are  treated  with  the 
escape  probability  model  of  reference  5.  This  model  allows 
for  radiation  interactions  between  all  grid  points  used  in 
the  code.  The  escape  probability  model  gives  an  accurate 
description  of  radiation  transfer  and  is  much  faster  than 
an  angle  and  frequency  dependent  radiation  transfer  scheme. 

The  original  numerical  escape  probability  code  was 
vectorized  and  optimized  with  a  resulting  increase  in  speed 
on  the  ASC  computer  of  roughly  a  factor  of  10.  The  run 
time  for  the  code  has  been  further  reduced  by  treating  only- 
selected  lines  with  the  escape  probability  model,  leaving 
the  other  lines  to  be  treated  in  the  optically  thin  limit. 
For  lines  which  are  always  optically  thin  or  for  lines 
which  are  not  important  energetically  or  for  diagnostic 
purposes,  this  is  a  reasonable  procedure.  For  comparison 
purposes,  a  model  for  the  continuum  transport  using  a 
Rosseland  mean  opacity  was  also  developed. 

Numerous  model  calculations  have  been  made  during 
the  course  of  the  year  for  both  aluminum  and  argon  plus  DT. 
One  example  of  this  is  the  study  of  the  implosion  and  blow- 
off  of  a  spherical  solid  aluminum  shell  Figure  22  shows 
the  aluminum  K  and  L  shell  lines  used.  For  this  problem, 
line  cooling  was  insignificant  due  to  the  high  densities, 


28 


Figure 

22 

ATOM  in  LEVEL 

STRUCTURE 

* 

A1  XIII 

AT  XII 

'  A1  XI 

Al  X 

A1  IX 

n  *  2 

2s  3S 

2p 

2s2p  3P 

2s2p2  4P 

n  »  3  1 

2p  3P 

3s 

2s2p  h 

2s2 p2  2S,2P,2D 

n  *  4 

2s  h 

3p 

2  p2  ^  P 

2s23s  2S 

n  *  5 

2p  4P 

3d 

2p2  4D 

2s23p  2P 

3  (sing) 

4  (doub) 

2p2  h 

2s23 d  2D 

4  (sing) 

5  (doub) 

2 s 3 s  3S 

2s2p3d  4P,4D 

5  (sing) 

6  (doub) 

2s3p  3P 

2s2p(3P)3d  2P,2D,2F 

2s2p  3P 

1s2p2  20 

2s3d  30 

2s2p(1P)3d  2P,20,2F 

2p2  2P 

2s3s  h 

'  2s2p4d  4P,4D 

.  ^ . 

2 s3 p  4P 

2 s 3 d  ]Q 

A1  VIII 

A I  VII 

AT  VI 

2p2  h,h 

2 s22 p3  2P,20 

• 

2sZ2p4  4S,4D 

2s2p3  5S 

2 s2 p4  4P 

2s2p5  3P 

2s2p3  3S,3P,3L) 

2s2p4  2S,2P,2D 

Zs2p5  h 

2s2p3  ]?,]Q 

2s22p23s  4P. 

2s22p33s  3S,3P,3D 

1. 

2s2Zp3d  3Pf3D, 

3f 

2s22p23s  2D 

2s22p33s  VjD 

t 

2s22p3d  ]P,h, 

’f 

2s22p23d  4P,4D 

2s22p3(4S)3d  3D 

*  VI 

2s2p2(20)3s  30 

2s22p(3P)3d  2P 

,2o,2f 

2sz2p3(z0)3d  3S,3P,3D,3F 

'  e 

2s2p23d  5P 

2s22p(1D)3d  2S,2P,2D,2F 

2 s22 p3 (2D)3d  ^Jd^F 

2s2p2(4P)3d  3P 

,3d,3f 

2s22p3(2P)3d  3P,3Q,3F 

2s2p2 (2D)3d  3S,3Pf30,3F 

2s22p3(2P)3d  WV 

Ail 

A 1  IV 

AT  III 

2s2p6  2S 

2s22p53s  3P 

2s22p63s  2S 

2 s2 2 p4 3 s  2P 

Zs22p53s  V 

2s2Zp63p  2P 

2s2Zp43s  20 

2s22p53p  ]D 

2s22p64s  2S 

2s22p43s  2S 

2s22p53p  3P 

Zs22p64p  2P 

2s22p4(3P)3d 

2p,2o,2f 

2s22p53d  3P,1P 

Zs22p64d  20 

2s22p4 ( 1 D)3d 

2s,2p, 2d,2f 

2s22p53d  3D 

Zs22p64f  2F 

2s22p54s  3P 

L_ 

2s22p 

54s  4P 

and  line  emission  was  post-processed.  Figures  23  and  24 

show  temperature  and  density  profiles  respectively  at 

various  times.  Figure  25  shows  net  cooling  profiles. 

The  solid  shell  compressed  the  low  density  aluminum  filler 

in  this  calculation  to  a  peak  central  density  of  roughly 

five  times  solid  aluminum.  Figures  26  and  27  show  spectra 

at  T  =  4.4  x  10  s  (early  in  the  evolution  of  the  system 

before  the  shell  has  significantly  been  affected  by  the 

radiation  absorbed)  and  Figures  28  and  29  show  spectra  at 
-9 

T  =  2.34  x  10  s  (peak  compression).  It  is  clear  from  the 
comparison  of  Figure  26  and  28,  and  of  27  and  29  that  the 
peak  compression  temperature  is  considerably  higher  than 
that  for  the  earlier  time.  The  figures  highlight  the  grea 
amount  of  detail  carefully  included  in  the  model  of  the 
radiation  field. 
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Self-absorption  of  heliumlike  satellite  lines  in  high-density  fusion  plasmas 

D  Duston*  und  J  Davis 

P-u\ma  PjJulutn  i iroup.  P'.j^ma  Phv no  Dins. on  kjiut  Roi'U'rli  Larti>rjitirv  14 -list'll;:.  r  [)  C  2'/*"' 

iRc-eneii  14  September  I^'Si 

The  high  densities  dbiained  in  revent  User  fusion  experiment:.  nave  erejted  a  need  tor  additional  plasma 
density  diagnostics  The  ratio  ot  '.no  heliumlike  satellite  lines,  the  l>2p  P  -ip  '  P  and  ’.-.Hi  S  -li  Ip  P 
transitions,  has  shown  promise  as  a  spectral  Jiagn.istis  when  iaser-imptoded  microbailuonv  are  seeded  with 
..  edium- atomic- weigh [  gases  In  this  study  results  are  presented  trom  a  colhsional-radiatise  .duration 
dynamics  model  with  photoeacitatton  processes  included  that  indicate  that  the  emission  from  these  satellite 
lines  is  strongly  affected  by  opacity  in  density  and  temperature  regimes  common  to  plasmas  to  which  this 
diagnostic  could  be  applied  effectively  The  radiation  emission  attenuated  by  photon  reabsorption  is 
presented  for  neon,  aluminum,  and  argon  plasmas  and  compared  with  results  predicted  wnen  the  calculation 
ts  undertaken  in  an  optically-thm-plasma  approximation  The  jpacitv  effects  are  seen  to  .ause 
multivaluedness  in  the  line  ratio  at  veveral  temperatures,  and  an  overall  ;oss  .'f  -envitintv  f  'he  ratio  with 
density  is  predicted 


In  recent  years  much  work  has  been  devoted  to 
understanding  radiation  spectra  emitted  from 
dense  high-temperature  plasmas.  Although  most 
of  the  available  experimental  spectra  have  been 
obtained  m  the  past  from  laser-produced  plas¬ 
mas.1  exploding  wire  plasmas,1  or  gas-puff  2 
pinches.1  data  have  recently  been  obtained  by- 
seeding  glass  microballoons  with  a  medium-2 
gas  to  produce  x  rays  from  bound-bound  transi¬ 
tions  occurring  during  a  high-temperature  com¬ 
pression  of  an  imploded  pellet.1  ’  The  purpose 
of  these  experiments  is  to  make  full  use  of  the 
x-ray  lines  as  diagnostic  indicators  of  the  tem¬ 
perature  and  density  history  of  the  plasma  con¬ 
tained  by  the  laser- imploded  glass  shell.  Un¬ 
fortunately.  the  densities  being  attained  at  pre¬ 
sent.  typically  10”-10JJ  electrons  per  cubic  cen¬ 
timeter.  are  in  a  somewhat  difficult  regime  for 
their  determination.  Several  techniques  are  avail¬ 
able  to  the  experimentalist,  e.g..  Stark  broa¬ 
dening  of  selected  lines.'1  merging  of  the  series 
limit,7  and  such  selected  line-intensity  ratios  as 
the  intercombination-to- resonance  line  ratio  in 
heliumlike  ions.’  In  the  past,  these  techniques 
have  been  applied  to  laboratory  plasmas  under 
the  assumption  that  these  lines  are  optically  thin. 
However,  the  validity  of  the  density  determina¬ 
tions  strongly  depends  on  whether  these  high- 
density  plasmas  are  indeed  optically  thin  and. 
if  not,  whether  opacity  effects  vitiate  their  use¬ 
fulness  as  density  diagnostics. 

In  this  paper,  a  study  is  made  of  the  radiation 
emitted  by  two  doubly  excited  levels  of  the  helium- 
like  ion  which  show  up  as  satellite  lines  to  the 
hydrogenlike  resonance  line,  and  the  line  ratio 
formed  by  them  as  used  to  determine  the  plasma 
density.  The  sensitivity  range  of  this  diagnostic 
extends  over  several  orders  of  magnitude  in  ion 


density,  but  is  centered  at  approximately  10J1 
ions,  cniJ  for  materials  presently  used  for  micro- 
balloon  seeding,  making  it  a  logical  candidate 
for  application  in  density  determinations  of  com- 
pressed-laser-pellet  experiments.  A  theoretical 
calculation  of  this  line  emission  as  it  is  genera¬ 
ted  in  optically  thin  neon,  aluminum,  and  argon 
plasmas  is  performed.  In  addition,  this  line  emis 
sion  is  also  calculated  in  plasmas  where  photo 
excitation  is  at  least  equally  important  us  col- 
lisional  excitation  in  order  to  assess  the  effect 
of  opacity  on  the  line- intensity  ratio  at  high  den¬ 
sity. 

I  THKORLTH  \L  MODtL 

It  has  been  shown  from  recent  work > 10  that,  in 
the  optically  thin  approximation,  the  line  ratio 
under  consideration  is  independent  of  the  hydro- 
geniike  ground-state  population  and  depends  only 
on  the  electron  density,  electron  temperature, 
and  atomic  parameters  associated  with  the  levels 
being  studied.  We  define  an  abbreviated  level 
structure  as  shown  in  Fig.  1.  and  analytically 
calculate  the  line  ratio  in  terms  of  the  popula¬ 
tions  of  the  two  doubly  excited  states  of  the  he¬ 
liumlike  ion.  the  2v2b  ‘P  and  2b"  ‘P  levels,  and 
their  radiative  transitions,  lsls  'S  —  Z s2p  'P  and 
Is  2p  ’P  —  2pi  'P. 

First,  the  atomic  rate  equations  for  each  of  the 
doubly  excited  states  are  written 

U) 

,J7f  -  cI.v,.v|w.v.v(.vl-lr/ii1.r.v().v1. 

Here.  .V  denotes  the  state  density.  C  is  the  elec¬ 
tron-capture  rate  coefficient.  T  is  the  autoioni- 
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model.  The  electron- capture  rate  coefficient  is 
obtained  by  detailed  balance  of  the  autoionization 
rate,  giving,  for  example,  for  state  2. 

C*  =  r,  C?2/2  gt)  (2  **  */*»« ,  *r.) 1  ■ 2  *  ’*  • ' r* ,  (4 1 

where  ?  is  the  statistical  weight  and  £,  is  the 
energy  of  the  doubly  excited  state  above  that  of 
the  hydrogenlike  ground  state.  The  electron  col- 
lisionai  excitation  rate  coefficient  coupling  the 
two  levels  has  been  calculated  by  the  method  of 
distorted  waves,12  and  the  deexcitation  rate  coef¬ 
ficient  also  obtained  from  detailed  balance  of  the 
excitation  rate  coefficient 


FtG  1.  Simplified  level  structure  for  analytic  treat¬ 
ment  of  the  satellite  line  emission  from  heliumlike  ions. 

zation  rate,  .4  is  the  spontaneous  decay  rate,  X 
is  the  electron-excitation  rate  coefficient,  Y  is 
the  deexcitation  rate  coefficient,  X,  is  electron 
density,  and  the  subscripts  g,  1,  and  2  indicate 
the  Is,  2s2p,  and  2 p1  states,  respectively.  Note 
that  at  the  densities  under  consideration  here, 
the  electron  and  ion  collisions!  rates  coupling 
the  individual  j  components  of  each  of  the  two 
doubly  excited  levels  are  sufficiently  large  that  it 
is  assumed  that  they  will  each  be  populated  ac¬ 
cording  to  their  statistical  weights,  allowing  for 
their  representation  as  the  two  n  levels  in  Fig. 

1  (the  j  degeneracy  has  been  removed).  Under 
the  assumption  that  the  atomic  rates  are  faster 
than  the  typical  hydrodynamic  time  scale  of  the 
plasma  [collisional- radiative  equilibrium  (CRE) 
model],  the  time  derivatives  in  Eq.  (1)  are  set 
to  zero  and  the  solutions  for  .Vt  and  X2  are  given 
by 


.V,  = 


_  -Vi-V v,ct 


XY 


v  .  •Vfl.c,r+-VW3 

9,9,-X^XY 


(25 


where  0,  =  A,  -  T,  *  XX,  and  92=A2+r2*  YX,.  The 
line  intensity  is  now  given  by  .Y,4A£,  where  A£  is 
the  transition  energy.  Forming  the  line  ratio  and 
solving  from  Eq.  (2),  one  can  readily  show  that 


R  = 


X,A,  A£„  ,4,A£„ 


XiAlXEl  AlXEl  F.V 


\l~  YX,(Cl/C2)  +  ‘i,  I 


(3) 


Y  =X  (gi/gj)ei£'  T’  ■  (5) 

where  A£  is  the  energy  separation  of  the  two 
doubly  excited  states. 

As  a  theoretical  check  on  the  calculation,  ap¬ 
proximations  can  be  made  in  Eq.  (3)  to  obtain  the 
values  for  the  ratio  in  the  local  thermodynamic 
equilibrium  (LTE)12  and  corona  limits.14  In  the 
LTE  (high-density)  limit,  it  is  expected  that  col- 
lisiona!  rates  will  dominate  over  competing  ra¬ 
diative  rates.  Setting  the  .4’s  and  T’s  to  zero 
in  the  large  parentheses  in  Eq.  (3),  we  obtain  the 
LTE  limit  (noting  that 
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3  I  T, 


(6) 


In  the  corona  limit,  radiative  processes  are 
dominant  over  competing  collisional  rates.  If  we 
form  the  ratio  X3/Xl  from  Eq.  (2)  and  let  go  to 
Oj  Ti  and  go  to  a,  =  .4,+  T2.  we  get 


-Vj  X,Y  Y  .\-]CiY*X,Cl  a2  ' 


(V) 


Ignoring  the  .V2  XY  term  in  the  numerator,  fac¬ 
toring  out  an  al/XeY  term  from  the  bracket,  com¬ 
bining,  and  expanding  out,  assuming  that  C,a2 
»  X'C2  Y,  we  obtain  the  following  expression  in 
the  corona  limit: 


2  r2(.4t>r,)  ,-i£  r, 

a,±e1  r^.Aj.r.) 


(8) 


where  we  have  eliminated  C,  using  Eq.  (4).  It 
is  expected  that  Eq.  t6)  and  t.8)  will  define  the 
high-  and  low-density  limits  of  the  optically  thin 
line  ratios  studied  here:  comparisons  will  be 
made  in  Sec.  HI. 


II  IONIZATION- RADIATION  MODEL 


and  hence  the  ratio  is  independent  of  Xf. 

For  this  study,  we  have  chosen  to  use  the  ra¬ 
diative  decay  and  autoionization  rates  of  Vainsh¬ 
tein  and  Safronova11  and  average  over  the  indi¬ 
vidual  j  components  so  as  to  be  applicable  to  our 


In  Sec.  I.  it  was  determined  that  the  1  <2t>  'P 
-2/>2JP  to  l.s2s  2S  —  2s2/>  'P  line- intensity  ratio 
was  independent  of  the  populations  of  surrounding 
ground  or  singly  excited  levels.  If  opacity  ef¬ 
fects  are  included  in  the  model,  this  is  not  the 
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case.  When  the  optical  depths  of  the  two  satel¬ 
lite  lines  exceed  1,  photoexcitation  is  expected  to 
significantly  affect  the  populations  of  the  doubly 
excited  states,  coupling  them  to  the  lower  states 
of  the  satellite-line  transitions,  the  ls2sJS  and 
Is 2 f>3P  levels.  Although  previous  authors910 
have  discounted  radiative  absorption  processes 
in  forming  this  ratio  with  the  argument  that  the 
lower  states  of  the  lines  are  excited  states  of 
the  heliumlike  ion  and  are  expected  to  be  well 
down  in  population  from  the  ground  state,  and 
hence  generate  small  optical  depths  for  these 
lines,  we  have  found  this  not  to  be  the  case  at 
densities  and  temperatures  where  the  line  ratio 
is  expected  to  be  a  valuable  diagnostic.  Hence 
some  means  must  be  included  in  the  calculation 
to  take  account  of  radiative  transfer  and  photo¬ 
excitation  when  the  satellite  lines  become  optical¬ 
ly  thick. 

As  a  basis  for  these  calculations,  we  use  an 
atomic  model  based  on  rate  equations  describing 
the  processes  populating  ground  states  and  selec¬ 
ted  excited  states  for  each  element  studied.  The 
processes  in  this  model  include  collisiona!  ioni¬ 
zation;  coilisional.  radiative,  and  dielect ronic 
recombination;  coilisional  excitation;  coilisional 
and  spontaneous  radiative  deexcitation;  and  stim¬ 
ulated  absorption  and  emission.  The  equations 
are  solved  time  independently  for  a  homogeneous 
plasma  of  constant  electron  temperature  and  ion 
density  to  yield  fractional  populations  of  the  sta¬ 
tes.  electron  density,  and  lir.e-emission  inten¬ 
sities  self-consistentiy  with  the  radiation-trans¬ 
port  calculations. 

The  photon-absorption  processes  are  taken  into 
account  by  a  phenomenological  frequency-dif¬ 
fusion  model  which  allows  an  emitted  photon  to 
escape  an  absorbing  plasma  by  successive  col¬ 
lisions  with  ions,  eventually  scattering  into  the 
optically  thin  wings  of  a  broadened  line  profile 
and  escaping  the  plasma.  A  photon-trapping  fac¬ 
tor  can  be  calculated  as  a  function  of  the  optical 
depth  through  the  plasma  at  a  frequency  corres¬ 
ponding  to  the  line-center  frequency  of  the  photon 
transition.  This  factor  can  thus  reduce  the  line 
source  function  depending  on  the  opacity  of  the 
plasma.  It  is  also  included  in  the  rate  equations 
to  allow  for  the  optical  pumping  of  the  upper  le¬ 
vels  of  thick  lines.  The  optical  depth  of  the  line 
is  defined  as 


effective  photon  path  length  through  the  plasma. 

.V  is  the  lower-state  population.  <p  is  the  line 
profile  function,  T,  is  ion  temperature,  and  M 
is  the  ion  mass.  In  order  to  accurately  determine 
the  .\y s  and  hence,  the  optical  depths,  a  number 
of  excited  states  were  included  in  the  atomic  mo¬ 
del  in  addition  to  ail  the  ground  states:  the  n 
=  2, 3 ,4, 5  levels  in  the  hydrogenlike  ion,  the 
2s  lS,  2s  3S,  2p  lP,  2p  3P,  3-singlet,  4-singlet,  and 
5-singlet  levels  in  the  heliumiike  ion,  and  the 
2p,3s,3p,3d,  and  4rf  levels  in  the  lithiumlike  ion. 
A  complete  discussion  of  the  method  of  calcula¬ 
ting  rate  coefficients,  atomic  structure,  and  the 
opacity  model  has  been  presented  in  previous 
studies.15'16 

111.  RESULTS 

A.  Optically  thin  approximation 

In  order  to  obtain  a  clear  understanding  of  the 
role  played  by  collisions  in  determining  this  in¬ 
tensity  ratio,  and  also  to  compare  our  predictions 
with  previous  work,  an  initial  study  of  the  forma¬ 
tion  of  the  level  populations  in  the  optically  thin 
approximation  was  undertaken.  Hence  no  re¬ 
absorption  of  photons  and  subsequent  optical  pum¬ 
ping  was  allowed  in  obtaining  these  lirst  calcu¬ 
lations.  Since  it  has  been  shown  in  Sec.  II  that 
the  line  ratio  is  independent  of  all  but  the  2 s 2/> 
and  2 p2  populations  in  an  optically  thin  plasma, 
the  results  are  merely  the  solution  to  £q.  '31. 

The  ratio  is  depicted  graphically  in  Figs.  2-4  for 
neon,  aluminum,  and  argon  plasmas,  shown  at 
various  temperatures  over  a  wide  range  of  ion 
densities.  As  can  be  seen  from  the  figures,  each 
plasma  produces  the  same  basic  behavior  in  the 
line  ratio:  a  relatively  constant  low-  and  high- 
density  limit,  connected  by  a  rapidly  ascending 
transition  region  over  which  the  ratio  provides 
its  diagnostically  valuable  density  dependence. 


where  the  Doppler  width  is  given  by 

±i>0=(2Ti /Mr2)1  .  (10) 

Here,  f  is  the  dipole  oscillator  strength,  L  is  the 


FIG.  2.  •!  —  JsJ  p  JP>  intensity 

ratio  for  a  neon  plasma  in  the  optically  thin  approxima¬ 
tion. 
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FIG.  3.  (ls2pJP  — 2p“p>/(ls2s3S  — 2s2p3P>  inten¬ 
sity  ratio  for  an  aluminum  plasma  in  the  optically  thin 
approximation. 

The  2s2p  level  is  preferentially  populated  over 
the  2 p1  level  at  low  densities  since  the  ratio  of 
the  radiative  decay  rate  to  the  autoionization  rate 
for  the  2s2p  state  is  smaller  than  that  of  the  2p2 
state.  However,  the  electron  collisions  coupling 
these  two  states  become  increasingly  important 
in  the  redistribution  of  the  populations  as  the  plas¬ 
ma  electron  density  increases,  resulting  in  an 
increase  in  the  2 p2  level  relative  to  the  2s2p  level 
and  a  corresponding  increase  in  the  line-inten¬ 
sity  ratio.  Finally,  in  the  high-density  limit,  the 
two  levels  reach  a  state  of  statistical  equilibrium 
relative  to  each  other  with  the  population  ratio 
given  by  e'48 '  T«  and  the  line  ratio  given  by  Eq. 

(6). 

Note  that  for  purposes  of  density  diagnostics  in 
laser- imploded  pellets  doped  with  a  small  amount 
of  neon,  the  most  effective  region  for  density  pre¬ 
dictions  is  from  3  x  IO19  to  3  x  io31  ions /cm3, 
while  for  argon-seeded  pellets,  the  region  is  from 
3  x  IQ20  io  5  x  io»  ions/cm3.  Thus,  with  the  ap¬ 
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FIG.  4.  Ils2p1p  —  2/>nP>/!ls2s!S—  2s2p3PI  intensity 
ratio  for  an  argon  plasma  in  the  optically  thin  approxi¬ 
mation. 

propriate  gas  doping,  it  is  conceivable  chat  this 
technique  would  be  an  effective  diagnostic  over 
three  orders  of  magnitude  in  density,  provided 
the  pellet  fuel  can  be  heated  sufficiently  to  attain 
temperatures  where  the  satellite-line  emission 
is  intense  enough  to  be  measured. 

Also  shown  in  Fig.  3  are  the  results  of  an  earlier 
work  by  Vinogradov  et  al .9  The  model  used  in 
Ref.  9  leading  to  their  predicted  line  ratio  was 
applied  to  an  aluminum  plasma  and  the  ratio  cal¬ 
culated  for  an  electron  temperature  of  500  eV.  In 
fact,  Ref.  9  includes  only  the  emission  from  four 
transitions  in  the  2p2-ls2p  manifold  (3P2-3P2, 

for  experimental  reasons, 
while  our  calculation  includes  all  six  of  the  com¬ 
ponents  (the  3Pl-3P2  and  3P0-3P1  in  addition  to  the 
four  above)  from  the  2 p*  state.  The  curve  reflects 
the  calculation  from  Ref.  9,  however,  with  all 
six  lines  included,  and  is  seen  to  be  in  excellent 
agreement  with  our  CRE  model  predictions. 

In  order  to  verify  that  our  calculations  approach 
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TABLE  l.  Atomic  parameters  used  in  the  collisional  radiative  model  to  analyze  the  heli- 
umlike  satellite-line  structure 


Parameter  Neon  Aluminum  Argon 


A,  (A) 

12.307 

7.252 

3.763 

A,  (A) 

12.323 

7.26 

3.766 

A£,  (eV) 

1007,4 

1709.6 

3294.56 

A£,  (eV) 

1006.1 

1707. S 

3292.46 

-A£(eV) 

8.64 

11.7 

18.74 

A 

0.3935 

0.3974 

0.4036 

J  2 

0.2390 

0.2657 

0.2639 

A,  (sec-1) 

5.776X  1013 

1.68X  1013 

6.347X  io13 

A.  (sec-1) 

1.05  x  1013 

3.362  x  1013 

1.241  x  io14 

rj  (sec-1) 

1.353  x  1013 

1.357x  1013 

1.387X  to13 

T,  (sec-1) 

4.237x  LO11 

1.987X  10" 

1.352X  io13 

RLT£  (500  eV) 
R:or  (500  eV) 


1.73 

o.lOO 


1.953 

0.2+4 


1.S85 

1.031 
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the  correct  limits,  the  values  for  the  satellite¬ 
line  ratio  in  the  LTE  and  corona  approximations 
were  determined  in  accordance  with  the  para¬ 
meters  given  in  Table  1  and  Eqs.  (6)  and  (8).  Com¬ 
parison  between  the  values  obtained  in  Table  I  at 
electron  temperatures  of  500  eV  and  Figs.  2-4 
generated  by  our  model  also  indicates  good  agree¬ 
ment. 

B.  Optically  thick  plasma 

At  low  plasma  densities,  optical  depths  of  emis¬ 
sion  lines  are  small,  as  defined  by  Eq.  (9),  and, 
consequently,  the  photon-absorption  effects  are 
small  relative  to  competing  collisional  and  spon¬ 
taneous  radiative  rates  and  can  be  neglected.  As 
the  plasma  density  increases,  the  photoexcitation 
rates  become  comparable  or  even  exceed  com¬ 
peting  processes,  resulting  in  significant  changes 
in  the  population  densities  of  the  excited  states 
and  the  subsequent  bound-pound  radiation  emission 
from  these  levels.  In  previous  studies,  the  effects 
of  these  processes  on  the  determination  of  the 
hydrogeniike  satellite- line  ratio  have  been  ignored 
with  the  justification  that  the  lower  state  of  each 
of  the  relevant  transitions  is  itself  an  excited 
level,  and  hence  the  populations  of  these  levels 
should  not  be  large  enough  in  the  density  range 
of  interest  to  contribute  to  a  large  optical  depth 
in  plasmas  of  the  physical  dimensions  commonly 
encountered  in  the  laboratory. 

In  order  to  test  the  validity  of  this  assumption, 
we  performed  calculations  of  the  line  ratio  using 
our  radiation-ionization  dynamics  model  with 
frequency-diffusion  radiation  transport  in  a  way 
similar  to  those  for  the  optically  thin  approxi¬ 
mation  presented  above.  The  implementation  of 
an  expanded  level  structure  in  the  model  is  neces¬ 
sary  if  the  population  densities  of  the  lower  states 
of  each  transition  are  to  be  determined  accura¬ 
tely.  self-consistently  with  other  states  of  the 
neiiumlike  ion.  as  well  as  with  the  excited  and 
ground  states  of  neighboring  ions. 

As  was  mentioned,  the  line  source  function  in  an 
opaque  plasma  is  attenuated  by  a  factor  dependent 
on  the  optical  depth  and  on  the  mechanism  by 
which  the  emission  line  is  broadened.  In  this 
study,  we  have  chosen  to  describe  the  optically 
thick  line  shape  by  a  Voigt  line  profile,  taking 
into  account  the  collisional  and  radiative  broa¬ 
dening  of  both  the  upper  and  lower  levels  com¬ 
prising  the  transition  of  interest.  At  very  high 
densities.  Stark  profiles  are  necessary  to  de¬ 
scribe  the  line  shape  accurately.  This  problem 
is  currently  being  addressed  by  Jacobs  and  Davis, 
but  no  profiles  have  been  calculated  as  of  this 
writing. 
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FIG.  5.  'lsJp'P  —  2p:*Pi'ils2s3S—  ls2p  :P>  inten¬ 
sity  ratio  for  a  neon  plasma  of  radius  1U0  -  ir. — with 
opacity  effects. 

The  results  of  the  calculations  with  opacity  ef¬ 
fects  included  are  presented  in  Figs.  5-7  for  a 
neon,  aluminum,  and  argon  plasma,  respectively 
The  studies  were  performed  for  plasmas  with  a 
radius  of  100  ^m  over  a  range  of  ion  densities 
and  electron  temperatures  identical  to  those  used 
for  studying  this  line  ratio  in  the  optically  thin 
approximation  iFigs.  2-4). 

Upon  comparison  of  the  three  figures,  certain 
aspecs  of  the  calculation  for  each  different  plas¬ 
ma  are  seen  to  be  similar,  and  can  be  discussed 
in  terms  of  a  general  plasma  radiating  in  this 
temperature  and  density  regime: 

(i)  The  line  ratio  exhibits,  initially,  a  decrease 
relative  to  the  opticailv  thin  result,  which  occurs 
at  an  ion  density  somewhat  higher  than  that  at 
which  departure  from  the  coronal  result  takes 
place.  Evidently,  the  \<2p-2p2  line  intensity 
must  be  affected  by  opacity  at  a  lower  density 
than  the  ls2.< -2s2/>  intensity,  despite  the  tact 
that  the  former  has  an  oscillator  strength  rough¬ 
ly  two-thirds  of  the  latter  transition.  This  effect 


FIG.  'I.  'IsC/i  ‘r  —  -p:iP<  u.v.'v  >s-  2.-.4/1  !P>  mien- 
si  tv  ratio  tor  in  Aluminum  plasma  of  radius  -m  — 
with  opacity  effects. 
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can  be  understood  in  terms  of  the  lower- state 
population  densities  as  determined  by  the  col- 
lisional  radiative  model.  At  the  densities  at  which 
the  opacity  effects  first  manifest  themselves,  the 
lower  state  of  the  former  transition,  the  ls2o  3P 
level,  is  more  populated  then  the  lower  state  of 
the  other  transition,  the  ls2s3S  level,  owing  to 
the  strong  collisional  mixing  of  these  levels  at 
these  large  electron  densities.  The  ratio  of  the 
statistical  weights  is  3:1,  but  the  ratio  of  these 
populations  ranges,  typically,  from  2:1  to  3:1 
above  1019  ions/cm3,  as  is  determined  by  the 
actual  collisional  and  radiative  rates  employed 
in  our  model.  Referring  to  Eq.  (9),  therefore, 
the  optical  depth  at  line  center  for  the  ls2p  -2p* 
transition  will  be  approximately  1.3- 1.9  times 
greater  than  that  of  the  ls2s  ~2s2p  transition, 
and,  consequently,  will  experience  photon  re¬ 
absorption  at  a  lower  density,  resulting  in  the 
initial  decrease  in  the  line  ratio.  Hence,  ac¬ 
curate  collisional  and  radiative  rates  are  neces¬ 
sary  for  a  model  to  accurately  predict  the  ion 
density  at  which  deviations  from  the  optically 
thin  approximation  will  first  appear. 

(ii)  The  line  ratio  will  exhibit  this  initial  de¬ 
crease  owing  to  opacity,  in  general,  at  a  higher 
density  for  a  higher-temperature  plasma.  For 
example,  in  Fig.  6,  the  aluminum  plasma  will 
yield  a  line  ratio  which  begins  to  flatten  out  at 
1021  ions  cm2  for  the  300-eV  case,  but  this  doesn’t 
occur  until  3  *  1021  ions/cm3  for  the  800-eV  plas¬ 
ma.  This  effect  is  simply  a  manifestation  of  the 
ion  temperature  dependence  of  the  optical  depth, 
which  can  be  seen,  in  Eq.  (9),  to  fall  off  as  2”l/3. 
Deviations  from  this  effect  can  occur  at  low  tem¬ 
perature,  however,  as  is  seen  in  Fig.  5  for  a 
neon  plasma  at  100  eV  and  in  Fig.  7  for  an  argon 
plasma  at  500  eV.  The  departure  from  the  opti¬ 
cally  thin  result  occurs  at  higher  densities  than 
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FIG.  ?.  •\s2piP  —  Zp13P)/<ls2s‘S  —  2s2p>P)  inten¬ 
sity  ratio  for  in  lrgon  plasma  of  radius  100  um — with 
opacity  effects. 


it  does  for  the  higher-temperature  plasmas.  What 
is  occurring  at  lower  temperatures  is  a  major 
change  in  the  abundance  of  the  heliumiike-ion 
populations.  Over  a  wide  range  of  temperature 
(at  constant  ion  density)  the  heliumlike  ion  is  the 
dominant  ion  population  in  the  plasma16;  hence 
the  Is 2 s 3S  and  ls2p3P  densities  will  be  large. 

At  low  temperatures,  however,  these  densities 
fall  rapidly,  yielding  to  lower  stages  of  ionization; 
consequently,  the  optical  depths  of  the  two  trans¬ 
itions  also  decrease  rapidly,  delaying  the  onset 
of  the  opacity  effects  until  higher  densities,  in 
spite  of  the  lower-temperature  plasma  effect  de¬ 
scribed  above. 

(iii)  The  line  ratio  will  converge,  at  high  den¬ 
sities.  to  an  LTE  limit  significantly  different 
from  that  predicted  by  the  calculation  when  done 
in  the  optically  thin  approximation.  Although  the 
upper-state  populations  of  the  two  transitions  will 
eventually  relax  to  the  same  LTE  ratio  (1:1)  in 
the  high-density  limit  regardless  of  the  opacity 
effects  (collision-dominated  plasma  regime),  the 
radiation  emission  will  be  quite  different  from 
the  optically  thin  prediction.  This  can  be  demon¬ 
strated  by  using  an  approximate  analytic  model,13 
in  which  the  radiation  is  attenuated  by  a  factor 
y,  where1718 

y  3=(a/T0)l/I ,  (ii) 

where  a  is  the  Voigt  damping  parameter,  pro¬ 
portional  to  the  sum  of  all  collisional  and  radiative 
rates  depopulating  the  upper  and  lower  levels  of 
the  transitions.  From  Eq.  (6)  the  LTE  line  ratio 
in  the  optically  thick  case  now  becomes 

RLT.E(thick)=RLTE(thin)(yJ/yl)  .  (12) 

Now,  from  Eq.  (9),  t0  is  proportional  to  the  os¬ 
cillator  strength  and  the  lower- state  densities, 
giving 

yj/hi  1  *.  (13) 

where 

U4) 

(u  and  l  refer  to  upper  and  lower  states)  and,  in 
the  LTE  limit. 

*^1535  /‘^1s3P  ~  t Sls2p  •  (15) 

Noting  that  the  wavelengths  are  approximately 
equal  and  inserting  the  appropriate  statistical 
weights,  one  obtains 

RLTE(thick)  =f?LTE(thin)(  a1A1  /dj.Aj)1  2.  (16) 

At  high  densities,  the  parameter  a  is  dominated 
by  collisional  rates;  we  have  calculated  the  re¬ 
quired  electron-collisional  excitation  rates  ide- 
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TABLE  II.  Electron  collisional  excitation  and  deexcitation  rate  coefficients  which  dominate 
the  collision  broadening  of  the  satellite  lines. 


T  ransition 

Neon  ill)"9) 

Aluminum  (10 -9) 

Argon  (10-9) 

300  eV  X 

5.7 

2.7 

1.43 

Y 

5.87 

2.81 

1.52 

X 

7.7 

3.5 

0.463 

y 

2.6 

1.24 

0.173 

500  eV  X 

4.53 

2.13 

1.14 

Y 

4.61 

2.18 

1.18 

X 

6.1 

2.82 

0.47 

y 

2.1 

9.65 

0.168 

300  eV  X 

3.63 

1.71 

0.914 

y 

3.67 

1.73 

0.936 

X 

4.67 

2.27 

0.460 

y 

1.64 

7.70 

0.16 

excitation  is  obtained  by  detailed  balance)  using 
our  distorted-wave  approximation,12  and  the  dom¬ 
inant  rates  are  those  depicted  in  Fig.  1  as  X,  Y , 
x ,  and  y  and  are  listed  in  Table  H.  Over  the  range 
of  temperatures  studied  here,  the  ratio  a2/ a t  is 
determined  to  be  nearly  constant  at  high  densities 
(see  Table  II)  with  values  of  0.63  for  neon,  0.64 
for  aluminum,  and  0.8  for  argon.  Using  the  values 
given  in  Table  I  and  Eq.  (16)  we  obtain,  for  the 
LTE  line  ratio  in  the  optically  thick  case,  1.05 
for  neon,  1.105  for  aluminum,  and  1.2  for  argon, 
which  are  seen  to  agree  quite  favorably  with  the 
model  results  in  Fig.  5-7.  Small  discrepancies 
occur  in  the  case  of  argon  since  (i)  the  ratio  a2  /at 
is  not  strictly  constant  with  temperature  owing 
to  the  dominance  of  the  doubly  excited-state  mix¬ 
ing  rates  (see  Table  n),  and  (ii)  the  argon  plasma 
has  not  yet  reached  a  state  of  complete  LTE  at 
1024  ions/cm3  (see  Fig.  7). 

Although  several  features  of  the  line  ratio  in  the 
various  plasmas  are  similar,  there  are  still  some 
contrasting  features  apparent  in  the  optically 
thick  results.  The  most  obvious  physical  effect, 
although  only  a  slight  one,  is  the  increase  in  the 
density  at  which  opacity  effects  occur  with  atomic 
number  Z ,  as  seen  in  Fig.  5-7.  This  finds  ex¬ 
planation  in  the  actual  scaling  of  electron  col- 
lisionai  excitation  rate  coefficients,  which  de¬ 
crease  in  magnitude  at  threshold  with  Z ,  A  de¬ 
crease  in  the  populating  rates  at  a  given  density 
corresponds  to  a  decrease  in  excited-level  densi¬ 
ties  and,  consequently,  the  optical  depths  of  the 
transitions  of  interest  here.  Thus  at  1021  ions/cm3 
and  a  nominal  electron  temperature  of  300  eV, 
the  l«2p  ‘P  state  of  aluminum  will  be  more  popu¬ 
lated  than  the  same  state  of  a  similar  argon  plas¬ 
ma.  Hence,  the  optical  depth  of  the  ls2p3P-2pa 
transition  will  be  “thicker”  in  the  aluminum  plas¬ 


ma  than  in  the  argon  plasma,  and  opacity  effects 
will  occur  at  a  lower  density. 

In  Fig.  5,  the  optically  thick  line  ratio  for  the 
neon  plasma  at  100  eV  displays  rather  anomalous 
behavior  at  high  densities.  What  is  occurring  here 
is  a  high-density  effect  decreasing  the  abundance 
of  the  heliumlike  neon  ion  until  the  satellite  lines 
are  no  longer  optically  thick  and  the  ratio  rises 
to  its  optically  thin  LTE  limit,  an  effect  similar 
to  the  one  described  earlier  due  to  low  plasma 
temperature.  Simply  stated,  an  increase  in  ion 
density  at  constant  temperature  usually  drives  a 
plasma  toward  a  state  of  lower  ionization.15  Once 
again,  the  heliumlike  ion,  owing  to  its  ''closed- 
electron-shell”  configuration,  maintains  a  rela¬ 
tively  large  population  density  over  a  wide  range 
of  densities.  However,  at  high  enough  densities, 
particularly  at  low  temperatures,  the  ion  abun¬ 
dance  shifts  toward  lower  charge  states  rapidly, 
generating  a  corresponding  drop  in  excited-state 
populations  and  resulting  in  a  line-ratio  behavior 
as  seen  in  Fig.  5.  The  same  behavior  would  be 
seen  for  aluminum  and  argon  if  curves  for  lower 
temperatures  had  been  of  interest  in  this  work. 
However,  the  onset  of  this  behavior  in  a  labor¬ 
atory  plasma  would  usually  be  accompanied  by  a 
significant  decrease  in  the  line  intensities,  making 
this  diagnostic  inaccurate  at  these  densities  due 
to  its  weak  appearance  in  the  spectra. 

The  most  striking  difference  between  the  opti¬ 
cally  thick  results  for  the  three  plasmas,  how¬ 
ever,  is  seen  as  a  single  shoulder  or  “hump” 
in  neon,  the  absence  of  a  shoulder  in  aluminum, 
and  a  double  shoulder  in  argon.  Although  many 
factors  affect  the  line- ratio  behavior  as  have  been 
discussed  in  detail  above,  the  Z  scaling  of  the 
value  of  the  LTE  optically  thick  limit  and  the  den¬ 
sity  at  which  LTE  is  reached  are  most  important 
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in  determining  this  shape.  For  the  argon  plasma, 
in  Fig.  7,  the  “double  hump”  occurs  as,  first, 
the  ls2/>  -  2/>I  intensity  is  attenuated  by  opacity, 
decreasing  the  ratio;  second,  the  ls2s  —  2  s  2p 
transition  becomes  thick,  attenuating  that  line 
intensity  and  thus  increasing  the  ratio;  third,  the 
plasma  becomes  collision  dominated,  the  line 
ratio  tends  toward  its  LTE  limit,  and  this  limit 
is  lower  than  the  value  the  ratio  has  attained  be¬ 
fore  entering  the  LTE  regime,  causing  the  ratio 
to  decrease  and  creating  the  second  hump.  In  the 
aluminum  plasma  in  Fig.  6,  LTE  is  reached  be¬ 
fore  the  second  shoulder  can  manifest  itself,  and 
the  LTE  value  for  the  ratio  is  approximately 
equal  to  the  value  the  ratio  has  attained  before 
LTE  is  reached,  allowing  for  only  a  very  mild 
shoulder  structure,  most  apparent  at  1200  eV. 

For  the  neon  plasma,  in  Fig.  5,  the  LTE  state 
is  also  reached  at  relatively  lower  densities,  but 
the  LTE  value  of  the  line  ratio  is  much  lower 
than  the  value  the  ratio  has  attained  before  en¬ 
tering  the  LTE  state,  hence  a  single  pronounced 
“hump”  is  evident. 

Since  the  occurence  of  these  shoulders  and  the 
value  of  the  LTE  limit  of  the  line  ratio  seem  to 
be  strongly  dependent  on  the  actual  atomic  para¬ 
meters  of  the  particular  plasma  of  interest,  pre¬ 
diction  of  a  general  shape  would  seem  difficult 
in  light  of  the  many  ways  the  parameters  in  Table 
I  scale  with  Z .  Unfortunately,  this  seems  to  in¬ 
dicate  that  each  element  requires,  to  a  certain 
extent,  individual  treatment  if  this  line  ratio  is 
to  be  applied  successfully  to  diagnose  laboratory 
plasmas  in  which  photoexcitation  is  an  important 
process, 

IV  DISCUSSION 

We  have  studied  the  heliumlike-satellite  line- 
intensity  ratio  with  our  sophisticated  collisional- 
radiative  ionization-dynamics  model  in  order  to 
assess,  within  the  context  of  an  optically  thick 
plasma,  the  value  of  this  ratio  as  a  plasma  den¬ 
sity  diagnostic  in  the  high-densitv.  high-tempera¬ 
ture  plasma  regime  encountered  in  recent  laser- 
pellet  interaction  experiments.  The  behavior  of 
the  Line  ratio  was  first  examined  under  the  as¬ 
sumption  of  an  optically  thin  plasma.  It  was  found 
that  the  density  profile  of  the  line  ratio  was  in¬ 
dependent  of  ion-state  populations  in  this  approxi¬ 
mation.  In  addition,  it  was  found  that  the  den¬ 
sity-sensitive  range  of  the  line  ratio  occurred 
at  higher  densities  with  increasing  atomic  num¬ 
ber,  extending  the  range  of  applicability  of  this 
diagnostic,  depending  on  the  choice  of  the  fill 
gas.  The  results  were  compared  with  other  in¬ 
dependent  work  and  with  limiting  values  deter¬ 


mined  from  simple  analytic  models,  and  in  all 
cases  the  agreement  was  quite  favorable. 

Armed  with  confidence  in  the  model  and  under¬ 
standing  gained  from  the  optically  thin  analysis, 
we  proceeded  to  investigate  the  effects  of  stimu¬ 
lated  absorption  and  emission  on  these  plasmas 
by  “turning  on"  the  opacity  calculation.  The  most 
obvious  result  from  this  study  was  that,  indeed, 
the  satellite- line  ratio  from  plasmas  of  this  type 
was  affected  by  opacity,  despite  the  fact  that  the 
lower  levels  of  the  two  transitions  are  excited 
states  of  the  heliumlike  ion.  Not  only  did  the  de¬ 
viations  from  the  optically  thin  result  occur  at 
densities  where  the  diagnostic  is  most  sensitive 
and  useful,  but  many  of  the  line  ratios  were  mul¬ 
tiple  valued  with  density,  making  accurate  diag¬ 
nostics  difficult.  As  further  analysis  illustrated, 
the  actual  shape  and  overall  behavior  of  the  line 
ratio  versus  density  was  a  sensitive  function  of 
the  atomic  parameters  pertaining  to  the  indivi¬ 
dual  plasma.  This  finding  led  to  two  inescapable 
conclusions:  (a)  accurate  theoretical  modeling 
of  this  diagnostic  in  the  “thick”  case  depends  on  a 
comprehensive  level  structure  within  the  CRE 
model  and  accurately  determined  cross  sections 
characterizing  the  collisional  processes  linking 
the  level  structure,  and  (b)  owing  to  the  complex 
interaction  of  the  various  atomic  parameters  in 
determining  the  level  populations  and  radiation 
field,  and  the  various  Z  dependencies  of  these 
parameters,  simple  Z  scaling  of  the  line- ratio 
behavior  in  a  reabsorbing  plasma  is  not  likely; 
rather,  each  plasma  will  have  to  be  characterized 
individually. 

In  light  of  the  results  obtained  for  the  intensity 
ratio  with  photoexcitation  processes  included  in 
the  calculation,  it  would  appear  that  opacity  ef¬ 
fects  may  cause  some  departure  from  optically 
thin  results  at  densities  presently  encountered 
in  pellet-fusion  experiments.  The  extent  of  this 
departure,  of  course,  will  depend  on  the  initial 
fill  pressure  of  the  admixture.  By  using  very 
low  fill  pressures,  the  density  of  the  additive 
gas  at  peak  compression  will  be  below  the  density 
at  which  opacity  effects  commence  to  affect  the 
ratio;  however,  the  trade-off  is  in  the  reduced 
number  of  emitting  ions,  resulting  in  emission 
lines  with  intensities  ’which  may  be  too  low  to 
measure  accurately  with  existing  spectroscopic 
methods. 

To  determine  the  degree  of  the  self-absorption 
effects  in  typical  microballoon  experiments,  a 
parameter  study  has  been  done  for  neon  and  argon 
where  the  source  size  has  been  scaled  with  the 
density  as  .V*1  1 .  Plasma  radii  were  determined 
by  using  the  law  of  mass  conservation  and  a  set 
of  initial  conditions  chosen  to  model  actual  micro- 
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FIG.  3.  {\s2p  3P  —  2p2  3P)/(ls2s  *P)  intensity  ratio 
for  a  neon  plasma  corresponding  to  initial  conditions  of 
a  130-pm-diam  microballoon  with  1  atm  of  neon  and  9 
atm  of  DT. 


FIG.  9.  Hs2p*P  —  9pJ3Pt'as2s JS—  2s2p}P)  inten¬ 
sity  ratio  for  an  argon  plasma  corresponding  to  initial 
conditions  of  a  230-jim-diam  microballoon  with  0.4  atm 
of  argon  and  30  atm  of  DT. 


balloon  experiments.  For  the  neon  study,  a  bal¬ 
loon  diameter  of  180  |im  and  fill  pressures  of 
1  atm  of  neon  and  9  atm  of  deuterium-tritium 
(DT)  (10%  neon  gas,  by  volume)  were  chosen  as 
typical  initial  conditions.5  The  argon  study  was 
modeled  after  recent  experiments19  which  em¬ 
ployed  230-p.m-diam  glass  shells  with  fill  pres¬ 
sures  of  0.4  atm  of  argon  and  30  atm  of  DT.  While 
the  DT  gas  does  not  affect  the  neon  satellite-line 
radiation  directly,  the  increase  in  electron  den¬ 
sity  due  to  the  ionization  of  the  DT  atoms  will 
affect  the  collisional  dynamics  which  determine 
the  level  populations  of  the  impurity  ions,  and 
must  be  taken  into  account.  This  results  in  an 
increase  in  the  electron  density  by  factors  of 
approximately  2  and  7  for  neon  and  argon,  re¬ 
spectively,  compared  to  that  of  the  seed  gas  alone 
at  the  temperatures  for  which  the  study  was  done. 

In  Fig.  8,  the  results  of  the  neon  study  are  shown 
for  electron  temperatures  of  300,  500,  and  700 
eV,  along  with  the  optically  thin  result  at  300  eV 
plotted  versus  neon-ion  density.  Similarly,  the 
argon  results  at  600.  800.  and  1200  eV  with  an 
optically  thin  curve  for  600  eV  are  shown  in  Fig. 

9  versus  argon-ion  density.  The  effect  of  photo¬ 
excitation  is  seen  to  cause  departure  from  the 


thin  result  at  implosion  diameters  of  between  20 
and  30  pm  in  the  case  of  neon  and  at  a  diameter 
of  about  45  pm  for  argon,  representative  in  both 
cases  of  compression  volumes  attainable  by  ex¬ 
isting  systems.  Clearly,  the  optically  thin  re¬ 
sults  are  not  accurate  at  higher  compressions. 

The  final  point  of  this  discussion  which  we  would 
like  to  emphasize  is  the  necessity  for  compre¬ 
hensive.  detailed  calculations  of  spectral  signa¬ 
tures  when  attempting  to  determine  parameters 
from  these  laboratory  plasmas.  Opacity  effects 
on  the  radiation  emission  need  not  render  the 
spectrum  useless  for  diagnostic  purposes;  multi¬ 
valuedness,  as  seen  in  Fig.  8  and  9,  need  not 
mitigate  line- intensity  ratios  for  applicability  in 
deciphering  the  plasma- radiation  signature.  No 
single  density  or  temperature  indicator  should 
be  relied  on,  but  rather,  searching  for  consis¬ 
tency  among  several  radiation  diagnostics  in  the 
spectrum  is  required  in  order  to  draw  accurate 
conclusions  regarding  the  general  plasma  state. 
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The  radiation  emission  from  dense,  high-temperature  aluminum  plasmas  has  been  studied  to  appraise  the 
influence  of  photoexcitation  on  several  spectral  diagnostics  currently  used  to  determine  average  plasma 
parameters.  The  plasma  state  is  described  by  a  set  of  collisional-radiative  rate  equations  characterizing  the 
competing  atomic  processes  occurring  in  a  homogeneous  plasma  volume  of  constant  density  and 
temperature.  Radiation  is  transported  through  the  plasma  by  the  method  of  frequency  diffusion,  a  technique 
which  employs  collision  times  of  the  photons  as  they  diffuse  through  the  wings  of  a  broadened  line  profile. 
Line  spectra  and  line-intensity  ratios  for  several  pairs  of  emission  lines  of  the  K  shell  are  calculated  in 
optically  dense  plasmas  and  compared  with  those  obtained  from  the  model  within  the  optically  thin 
approximation.  The  radiation-trapping  effects  are  shown  significantly  to  modify  the  behavior  of  the  optically 
thm  constant-density  or  constant-temperature  line  ratio  curves,  which  become  double  valued  in  some  cases. 
Comparisons  of  state  populations  obtained  from  this  model  and  the  local-thermodynamic-equilibnum  model 
are  also  shown. 


I.  INTRODUCTION 

It  has  been  suggested  from  recent  experiments1'2 
in  which  moderate-/  plasmas  have  been  heated 
to  high  temperatures  at  high  density,  that  the  re¬ 
sulting  radiation  emission  may  play  in  important 
role  in  the  dynamic  evolution  of  the  plasma  tem¬ 
perature  and  density  profile.  The  extent  to  which 
the  radiation  may  represent  a  significant  energy- 
loss  mechanism  is  a  moot  question,3  but  clearly, 
resonance  reabsorption  of  emitted  photons  in  an 
optically  thick  medium  will  both  redistribute  the 
plasma  energy,  spatially  and  spectrally ,5  and  can 
alter  significantly  the  radiation  signature  at  the 
diagnostic  detectors. 

In  an  effort  to  provide  theoretical  insight  into 
some  of  the  effects  of  radiation  and  ionization  on 
the  plasma  development,  several  atomic  models 
have  recently  been  developed5"5  which  take  into 
account  the  various  processes  occurring  in  high- 
temperature  plasmas  and  which  ultimately  affect 
the  radiation  emission.  While  a  complete  theor¬ 
etical  analysis  of  the  radiation  dynamics  and  its 
interplay  with  other  basic  plasma  processes  re¬ 
quires  a  full-scale,  detailed,  time -dependent 
model  coupling  ionization  dynamics  with  magneto¬ 
hydrodynamics,  we  feel  that  the  atomic  processes 
leading  to  the  radiation  emission  are  sufficiently 
complex  that  a  better  understanding  of  the  radia¬ 
tion  physics  would  be  more  realizable  if  the  model 
were  limited,  initially,  to  the  basic  ionization  and 
radiation  processes,  in  the  absence  of  convection. 

In  this  work,  we  present  results  from  our  ra¬ 
diation-ionization  model  as  applied  to  the  problem 
of  diagnostics  in  a  dense  high-temperature  alu¬ 
minum  plasma,  with  the  notion  that  systematic 
trends  in  the  emission  features  calculated  from 


this  model  will  lead  to  conclusions  about  the  na¬ 
ture  of  the  radiating  plasma.  These  results  will 
be  compared  with  those  obtained  from  different 
physical  models  as  well  as  with  other  investiga¬ 
tions.5"11 

II.  MODEL  PREDICTIONS 

The  ionization-dynamic  model  employed  in  this 
study  is  similar  to  our  earlier  model. 5,12,13  A 
general  set  of  time -dependent  atomic  rate  equa¬ 
tions  are  solved  for  a  stationary  volume  element 
of  plasma,  homogeneous  in  temperature  and  total 
ion  density.  Each  equation  describes  the  popula¬ 
tion  and  depopulation  of  a  particular  ion  (atom) 
species  as  determined  by  rate  coefficients  char¬ 
acterizing  the  various  collisional  events  the  par¬ 
ticles  experience.  Such  an  equation,  in  the  ab¬ 
sence  of  convection,  would  be 

(i) 

Ul  J  , 

where  Wif  is  the  rate  describing  a  transition  from 
state  i  to  j  (in  sec"1).  Each  ground  state  and  ex¬ 
cited  level  in  the  model  is  represented  bv  an  equa¬ 
tion  of  the  form  given  by  (1).  The  collisional  and 
radiative  processes  are  described  by  the  rate 
Wir  The  model  currently  includes  the  following 
atomic  processes:  Collisional  ionization,  col¬ 
lisional  (3-body),  radiative  (2-bodv),  and  dielec- 
tronic  recombination,  collisional  excitation,  col¬ 
lisional  and  spontaneous  radiative  de-excitation. 
The  collisional  ionization  rates  were  calculated 
using  a  prescription  suggested  by  Seaton.15  and 
the  radiative  recombination  rate  coefficients  were 
calculated  by  the  detailed  balance  of  photoioniza- 
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tion  rates15  employing  effective  Gaunt  factors.16 
The  dielectronic  recombination  rate  coefficients 
are  those  calculated  by  Jacobs  et  al.,1’’  while  the 
spontaneous  decay  rates  were  taken  from  the  NBS 
tables16  or  from  oscillator  strength  extrapola¬ 
tions.13  The  electron  collisional-excitation  rate 
coefficients  were  calculated  in  the  Born  approxi¬ 
mation19  for  the  hydrogenlike  transitions  and  by 
the  distorted -wave  method30  for  the  heliumlike 
and  lithiumlike  transitions.  The  rates  for  the 
collisional  processes  from  higher  to  lower  energy 
states  were  calculated  by  the  detailed  balance  of 
the  corresponding  upward  collisional  rate  coef¬ 
ficients. 

The  processes  of  photoexcitation  and  stimulated 
emission  are  taken  into  account  by  a  phenomeno¬ 
logical  transport  model  which  allows  photons  to 
escape  from  an  optically  thick  plasma  by  scat¬ 
tering  into  the  optically  thin  wings  of  the  absorp¬ 
tion  line  profile.  The  photon-loss  rate  depends 
upon  the  optical  depth  of  the  plasma, 


T*  = 


ire2  . 
mac  r'! 


Is, 

V  rrAu0 


(2) 


where  f  is  the  line  oscillator  strength,  .V.  is  the 
lower -state  density,  and  <p  is  the  line -profile 
function;  the  Doppler  width  is  given  by 

=(2Ti/Mc2)l/i  v0,  (3) 


where  f,  is  the  ion  temperature,  M  is  the  ion 
mass,  and  i/„  is  the  line-center  frequency.  The 
effective  photon  path  length  across  the  plasma  is 
given  by  L;  in  this  way  an  average  measure  of 
the  spatial  extent  of  the  optically  thick  plasma 
is  taken  into  account.  We  can  write  an  equation 
for  the  collisional  processes  undergone  by  pho¬ 
tons  in  going  from  state  j  to  state  i,  viz., 


d.W*  1  g.X .  /i  i  \ 

-(,r, 


(4) 


where  A  is  the  Einstein  coefficient,  ,V»  is  the  pho¬ 
ton  density,  and  f,  and  tt  are  scattering  and  es¬ 
cape  times  for  the  photon,  respectively.  The 
terms  on  the  right-hand  side  of  (4)  describe  spon¬ 
taneous  decay,  stimulated  emission,  stimulated 
absorption  and  escape.  We  can  rewrite  Eq.  (1) 
to  include  the  photon  transport  terms, 


dS 

It 


1  =  Z  -  E  l*vv, 


Assuming  that  the  photon  processes  occur  on  a 


time  scale  faster  than  the  other  collisional  pro¬ 
cesses,  the  time  derivative  is  set  to  zero  in  (4), 

•V*  is  solved  for  and  substituted  into  (5),  removing 
the  A lt  term  from  Wjt.  This  yields  the  probability 
that  a  photon  born  at  plasma  center  will  escape 
to  the  outside  ( XA  times  probability);  an  inte¬ 
gration  over  r  is  then  performed  to  average  over 
all  photons  born  throughout  the  plasma  volume 
(with  L  replaced  by  r  in  the  optical  depth)  to  ob¬ 
tain 


(6) 


We  incorporate  this  factor  in  the  rate  equations 
and  in  the  radiation  calculation, 


P „  =hvNHA  jty .  (7) 

Hence,  both  the  attenuation  of  the  radiation  by 
opacity  and  the  optical  pumping  effects  on  the  level 
populations  are  included  self-consistently.  To 
evaluate  y,  expressions  for  la  and  f,  have  been 
derived  as  a  function  of  the  broadening  mechanism 
influencing  the  line  profile.  For  a  naturally 
broadened  line,  the  plasma  is  assumed  to  be  opti¬ 
cally  thin,  r<l,  and  t,  =0.  For  Doppler  and  Voigt 
line  profiles,  respectively  the  expressions  can  be 
used  for  ft /fJ(13, 21,22 


lnT0)1/2(l  *  l/'2lnr0) 


(8) 


where  r0  is  the  optical  depth  at  line  center  (v  =  v0) 
and  a  is  a  damping  parameter  given  by 


an  =  ,  (9) 

4irAv0 

with  C  representing  collisional  rates.  Hence,  the 
Voigt  profile  is  broadened  by  both  collisional  and 
radiative  processes  from  emitting  and  absorbing 
states. 

Preliminary  studies  of  the  radiation  emission 
using  these  two  broadening  mechanisms  indicate 
that  under  most  plasma  conditions,  a  Doppler  line 
profile  will  result  in  more  attenuation  of  line  ra¬ 
diation  than  will  a  Voigt  profile,  primarily  because 
photons  can  escape  in  the  extended  wing  structure 
of  the  Voigt  profile. 

In  order  to  test  its  validity,  this  method  of  ra¬ 
diation  transport  has  been  applied  to  a  two-level 
atom  for  an  infinite  slab  plasma  medium23  and 
the  fluxes  and  radiation  field  departure  coeffi¬ 
cients  compare  quite  favorably  with  Hearn's  nu¬ 
merical  solution  of  the  exact  radiation  transport 
equation  in  this  geometry.24 

Aluminum  was  selected  for  study  because  of  our 
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FIG.  1.  Aluminum  level  structure  used  in  the  CRE 
model. 

previous  experience  and  understanding  of  aluminum 
models  and  the  availability  of  experimental  re¬ 
sults.25* 27  Basic  to  any  theoretical  calculation  of 
this  kind  is  the  atomic  level  structure;  for  aluminum 
we  have  included  the  14  ground  states,  4  excited  levels 
in  Alxill  («  =  2,j3, 4,  5),  7  levels  in  A1  XII  (2s  3S,  2  p 
3P,  2s  lS,  2 plP,  3-singlet,  4-singlet,  and  5- 
singlet),  and  5  levels  in  Alxi  (2 p,  3s,  3 p,  2d, 

4 d).  In  addition,  we  include  3  double-excited 
levels  (2s2p’iP,  2 p22P,  ls2p12D)  in  order  to  study 
some  prominent  satellite  lines  in  the  spectra  of 
high-temperature  aluminum  plasmas;  the  level 
structure  is  depicted  in  Fig.  1.  Note  that  we  have 
represented  individual  nl  components  of  the  hydro¬ 
genlike  ion  and  the  higher  Rydberg  levels  of  the 
heliumlike  ion  as  single  degenerate  n  states.  The 
collisional  rates  by  both  ions  and  electrons  which 
couple  the  =0  transitions  for  these  levels  have 
been  assessed  m  a  previous  work15  and  have  been 
found  to  be  sufficiently  large  over  a  wide  range  of 
densities  and  temperatures  that  the  various  l 
states  are  populated  according  to  their  relative 
statistical  weights.  Hence,  for  all  Al.\i,'T  states 
and  for  n*  3  in  A1XII ,  we  have  previously  avera¬ 
ged  the  rates  which  describe  processes  initiating 
or  terminating  on  an  individual  nl  level  over  all 
possible  values  of  l  for  a  given  n,  according  to 
the  statistical  weights.  Thus,  the  collisional -ex¬ 
citation  rate  from  the  Is2  state  to  the  n  =  3  singlet 
level,  for  example,  is  obtained  from  the  indivi¬ 
dual  rates  coupling  the  Is2  state  to  the  3s,  3 p,  and 
2d  singlet  levels. 
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HI.  RESULTS 

V  Line  spectrum  from  aluminum  plasmas 

One  reason  for  performing  calculations  of  this 
type  is  to  predict  level  populations  and  radiation 
emission  from  a  plasma  characterized  by  elec- 


FIG.  2.  Aluminum  abundance  curves  vs  electron  tem¬ 
perature  for  an  ion  density  of  ll)25  cm*3  (optically  thini. 

tr on  temperature,  ion  density,  and  size  As  an 
example  of  the  information  obtained  from  our  col- 
lisional-radiative  equilibrium  (CRE)  model  about 
state-density  populations,  an  abundance  curve 
showing  the  distribution  of  ionization-state  popu¬ 
lations  as  a  function  of  temperature  at  a  constant 
ion  density  of  102°  cm*3  in  the  optically  thm  ap¬ 
proximation  is  shown  in  Fig.  2.  The  higher  peaks 
corresponding  to  A1  rv  and  A1 XII  are  due  to  the 
large  ionization  energy  of  a  closed -shell  con¬ 
figuration.  These  peaks  are  usually  broader  in 
temperature  indicating  the  propensity  of  the  atom 
to  remain  in  such  a  configuration.  In  Fig.  3,  the 
“thin”  abundances  of  A1XI-A1XIY  are  shown  as  a 
function  of  ion  density  for  an  electron  temperature 
of  300  eV.  Note  that  over  a  wide  range  of  densi¬ 
ties,  the  closed-shell  A1XII  ion  is  again  the  do¬ 
minant  species. 

A  typical  x-ray  spectrum  as  calculated  by  our 


FIG.  .!.  Aluminum  abundance  curves  vs  ion  density  for 
an  electron  temperature  of  :100  eV:  CRE-thin  i—  t ,  CRE- 
thick  (-  .-i,  LTE  ( - 1. 
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FIG.  4.  Aluminum  line  emission  spectrum,  3—  i  A, 
for  an  ion  density  of  102°  cm'1,  electron  temperature 
of  500  eV  and  radius  of  100  |in, 

model  is  shown  in  Fig.  4  and  the  corresponding 
vacuum  ultraviolet  (vuv)  spectrum  in  Fig.  5,  for  a 
plasma  at  300  eV,  a  density  of  102'3  ions/cm3  and 
a  radius  of  100  jam  for  a  Voigt  profile.  All  lines 
in  Fig.  4  decay  to  the  Is  or  Is2  ground-state  levels, 
whereas,  in  Fig.  5,  all  transitions  terminate  on 
an  n  =  2  level.  The  AlXII  Balmer  series  lines, 
e.g.,  He  (2p‘P-3‘),  are  actually  a  manifold  of  the 
p-s  and  p-d  transitions,  but  since  levels  with 
principle  quantum  number  greater  than  2  are  con¬ 
sidered  as  degenerate  in  our  model,  only  their 
combined  intensities  are  shown  in  Fig.  5.  Also 
notice  the  satellite  structure  to  the  long -wave¬ 
length  side  of  the  hydrogenlike  and  heliumlike  re¬ 
sonance  lines;  these  satellites  are  due  to  the  tran¬ 
sitions  ls2s  3S-2s2p3P,  ls2p  ,P-2pi  3P,  and  ls22p 
2P-ls2p12D,  and  are  important  both  as  relatively 
strong  emitters  and  as  valuable  plasma  diagnos¬ 
tics. 

In  order  to  illustrate  the  significance  of  photo¬ 
excitation  on  the  radiation  emission  of  dense  plas¬ 
mas,  we  have  indicated  the  line  intensities  for  a 
plasma  with  identical  equilibrium  conditions,  but 


FIG.  5.  Aluminum  line  emission  spectrum  24—  oil  A 
for  an  ion  density  of  1020  cm'  5,  electron  temperature  of 
500  eV  and  radius  of  100  urn. 


calculated  in  the  optically  thin  approximation  by 
horizontal  lines  through  the  intensity  spikes  in 
Figs.  4  and  5.  The  differences  are  seen  to  range 
from  a  few  percent,  as  in  the  case  of  the  inter - 
combination  (IC)  line,  to  more  than  an  order  of 
magnitude  as  with  the  Hs  line.  In  general,  one 
would  expect  lines  with  large  optical  depths  to 
experience  a  decrease  in  intensity  due  to  the  large 
photoexcitation  cross  section.  This  is  ,  in  fact, 
the  case  for  the  heliumlike  Lyman  lines  and  li¬ 
thiumlike  lines  shown  in  the  spectra,  including 
the  H e(La)  satellite  line.  The  intercombination 
transition  is  unaffected  because,  in  spite  of  the 
fact  that  its  lower  level  is  a  highly  populated 
ground  state,  it  has  a  rather  small  spontaneous 
decay  coefficient  (2.5  x  1010  sec'1),  and  hence,  a 
relatively  small  optical  depth  (  r,  =0.23)  compared 
to  that  of  the  heliumlike  resonance  line  (  t,  =83). 

The  hydrogenlike  and  heliumlike  Balmer -series 
lines  show  a  universal  increase  in  intensity  when 
radiation  trapping  is  included  in  the  calculation. 
Since  their  lower  states  are  excited  states  of  the 
ions  (n=  2  levels),  their  optical  depths  are  small, 
and  little  photon  reabsorption  takes  place.  In  ad¬ 
dition,  the  Lyman-series  lines  are  so  optically 
thick  that  the  power  levels  are  “optically  pumped” 
by  the  stimulated  absorption  process  and  the  net 
effect  of  the  increase  in  upper -state  populations  is 
a  corresponding  increase  in  intensities  of  the  op¬ 
tically  thin  Balmer -series  lines  when  opacity  is 
included.  Note,  however,  the  Hell,)  line;  al¬ 
though  it  is  a  Lyman-series  transition,  and  has 
an  optical  depth  of  about  2.4,  the  increase  in  the 
5 p  population  due  to  pumping  more  than  compen¬ 
sates  for  the  attenuation  of  the  line  intensity  by 
photoabsorption. 

An  apparent  contradiction  to  the  above  argument 
concerning  lines  with  large  optical  depth  is  seen 
to  occur  in  Fig.  4,  where  it  is  shown  that  the  hy¬ 
drogenlike  Lyman-series  lines  are  universally 
stronger  in  the  thick  calculation  than  in  the  thin, 
even  though  the  optical  depths  of  the  a,  3,  and  y 
lines  are  28.7,  4.6,  and  1.6,  respectively.  In 
fact,  these  lines  are  somewhat  attenuated  by  photo- 
absorption,  but  an  actual  increase  in  the  thick 
intensities  is  caused  by  a  significant  enhancement 
of  the  .41  xm  ground  state  by  collisional  ionization 
from  the  AlXII  excited  states,  whose  populations 
are  now  optically  pumped.  In  order  to  verify  this 
effect,  a  calculation  was  performed  in  which  the 
radiation  transport  was  "turned  off"  in  the  hydro¬ 
genlike  ion.  The  intensity  of  the  H(  La)  line  was 
then  calculated  to  be  over  2  x  1013  W  cm3,  thus 
indicating  that  the  states  of  Alxin  receive  an  in¬ 
crease  in  population  from  .41  XII  states  in  the  op¬ 
tically  thick  case  which  exceeds  the  decrease  in 
line  intensity  from  the  hydrogenlike  Lyman  lines 
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due  to  photon  absorption.  Hence,  we  have  a  case 
where  it  is  possible  that  an  optically  thick  tran¬ 
sition  may  radiate  at  a  higher  intensity  than  would 
be  calculated  in  an  optically  thin  approximation, 
due  to  the  strong  collisional  coupling  between  ions 
and  significant  optical  pumping. 

B.  Diagnostic  line  ratios  in  (he  optically  thin  approximation 

Spectral  line  intensity  ratios  are  a  valuable  diag¬ 
nostic  tool  employed  by  spectroscopists  to  de¬ 
termine  the  average  temperature  and  density  of 
the  emitting  region  of  a  plasma.  Although  abso¬ 
lute  line  intensities  are  directly  linked  to  the  ra¬ 
diation  source  function,  they  are  difficult  mea¬ 
surements  to  make,  requiring  extremely  accu¬ 
rate  calibration  of  Che  detector.  Taking  the  ratio 
of  two  emission  lines  eliminates  some  of  the  more 
difficult  aspects  of  this  calibration;  in  addition, 
if  the  two  lines  lie  in  the  same  region  of  the  spec¬ 
trum,  technical  problems  due  to  nonlinearity  of 
film  or  crystal  response  are  lessened.  We  have 
chosen  to  study  the  behavior  of  five  line  ratios 
commonly  used  in  the  plasma  community  to  diag¬ 
nose  laboratory  and  astrophysical  plasma  sources. 
The  ratio  of  the  hydrogenlike  resonance  line 
(ls-2p)  to  the  heliumlike  resonance  line  (ls2-ls2p 
lP)  and  the  lithiumlike  satellite  line  (ls22p-ls2p2 
2D)  to  the  heliumlike  resonance  line  are  both  used 
to  determine  electron  temperature.  Similarly,  the 
ratio  of  the  hydrogenlike  resonance  line  to  the 
heliumlike  I3  line  (ls-’-ls3/5  lP)  and  the  hydrogen¬ 
like  L3  (Is -3/51  to  the  heliumlike  Z.,  (ls2-ls5p‘P) 
are  temperature  diagnostic  ratios.  Finally,  the 
ratio  of  the  heliumlike  intercombination  (IC)  line 
(ls2-ls2p  3P)  to  the  heliumlike  La  is  used  to  de¬ 
termine  density. 

In  general,  a  good  temperature  diagnostic  ratio 
is  relatively  insensitive  to  changes  in  density  while 
the  converse  is  true  for  a  good  density  diagnostic 
ratio.  We  have  plotted  the  density  dependence  of 
several  of  the  ratios  against  ion  density  from  cal- 


FIG.  6.  ls22p(2P)-Ls2p2(2D)/ls:-Ls2p{‘P)  line  ratio 
vs  ton  density;  optically  thin  t—  > ,  optically  thick  with 
radius  of  100  u  m  ( - 1. 
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FIG.  T.  ls-2p/ls'-ls2*>(‘Pl  line  ratio  vs  Ion  density; 
optically  thin  (->,  optically  thick  with  radius  of  100  nrr. 

( - !. 

culations  done  with  the  CRE  model  las  shown  in 
Figs.  6-10);  these  calculations  have  been  done 
assuming  an  optically  thin  plasma,  temporarily 
neglecting  reabsorption,  in  order  to  see  the  ef¬ 
fects  of  particle  collisions  on  the  ratios  as  the 
plasma  goes  from  a  coronal  to  a  collisional  re¬ 
gime.  The  AlXI  satellite  to  AlXII  L  a  ratio  in 
Fig.  6  is  relatively  flat  up  to  a  density  of  1021 
for  the  temperatures  displayed.  Above  this  den¬ 
sity,  the  ratio  increases  due  to  a  combination  of 
two  effects:  The  increasing  ratio  of  AlXII  ground  - 
state  population  to  A1  XIII  ground-state  population 
and  the  strong  collisional  coupling  between  the 
A1XIU  ground  state  and  the  ls2p  :P  level.  At  low 
densities,  both  the  Is 2/5  AlXII  state  and  the 
ls2p22D  AlXI  state  are  strongly  coupled  to  the 
AlXII  ground  state  by  electron  collisions:  con¬ 
sequently  the  ratio  is  nearly  constant  with  density. 
As  the  density  increases,  the  collisional  recom¬ 
bination  from  the  A1  XIII  ground  state  becomes  the 
dominant  mechanism  populating  the  Is 2 p  level, 
uncoupling  this  level  to  some  degree  from  the  A1 
XII  ground  state.  At  the  same  time,  the  fractional 
population  of  the  A1XIII  ground  state  decreases 
more  rapidly  with  density  than  the  AiXH  ground 
state  isee  Fig.  3).  This  leads  to  a  relative  in- 
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FIG.  S.  lino  ratio  vs  ion  densitv; 
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FIG.  9.  ls-2/>/ls:-ls5/>('Pi  line  ratio  vs  ion  density; 
optically  thin  i-i,  optically  thick  -nth  radius  of  100  um 

crease  in  the  ratio  of  the  l  >  2/> 3  to  ls2 fi  popula¬ 
tions  and  a  corresponding  increase  in  the  line  in¬ 
tensity  ratio.  Notice  that  the  increase  occurs  at  a 
slightly  lower  density  for  the  higher  temperature 
curves.  This  is  due  to  the  greater  fractional  popu¬ 
lation  of  the  hydrogenlike  ground  state  at  higher 
temperatures;  this  enhances  the  coupling  of  the 
U2plP  state  to  this  ground  state  making  its  popu¬ 
lation  less  strongly  coupled  to  the  heliumlike 
ground  state,  yielding  the  density-sensitive  re¬ 
sult  the  ratio  displays  at  higher  ion  densities. 

The  hydrogenlike  to  heliumlike  Lt  ratio  in 
Fig.  7,  the  hydrogenlike  La  to  heliumlike  L}  ratio 
in  Fig.  8,  and  the  hydrogenlike  Li  to  heliumlike 
ratio  in  Fig.  9  all  display  a  definite  density 
effect  at  the  temperatures  studied.  As  seen  from 
Fig.  3,  the  distribution  of  the  heliumlike  ground 
state  is  flat  at  low  densities  and  drops  off  at  high¬ 
er  densities;  the  hydrogenlike  ground  state  is 
more  peaked  with  a  maximum  near  l&-  ions.  cm3. 
The  line  ratios  in  Figs.  7,  3,  and  9  reflect  these 
distributions  rather  plainly,  with  definite  maxima 
occurring  in  the  1021-1023  ions  cm3  range.  It  is 
interesting  to  note  that  while  these  ratios  possess 


FIG.  10.  1s--1ju'/m3P)/15:-1s2/>(‘F1  line  ratio  vs  ion 
density;  optically  thin  (-),  optically  thick  with  radius 
of  100  pm 
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this  undesirable  density  dependence,  a  knowledge 
of  the  density  to  withm  an  order  of  magnitude  al¬ 
lows  one  to  predict  an  electron  temperature  from 
the  curves  with  an  uncertainty  of  usually  less  than 
100  eV.  Since  most  experimental  plasmas  are 
able  to  be  diagnosed  easily  to  densities  within  this 
range,  these  line  ratios  are  still  valuable  indica¬ 
tions  of  the  plasma  temperature. 

The  intercombination-to-resonance  line  ratio 
in  Fig.  10  is  a  well-known  density  diagnostic33 
used  with  reasonable  success3 '•3')  for  laboratory 
plasma.  The  intercombination  line  is  usually  a 
very  weak  line  compared  to  the  heliumlike  re¬ 
sonance  line  due  to  a  relatively  small  radiative 
decay  rate.  As  Z  increases,  however,  this  coef¬ 
ficient  increases  to  within  less  than  a  few  orders 
of  magnitude  of  the  decay  rate  of  the  resonance 
line,  and  a  reasonable  lme  ratio  measurement  can 
usually  be  made  for  Z  s  11.  At  low  densities,  both 
the  ls2p  lP  and  Is 2/5  'P  levels  are  populated  by- 
electron  collisions  from  the  ground  state  and  de¬ 
populated  by  spontaneous  decay,  yielding  a  rather 
constant  ratio  at  corona  densities.  As  the  density- 
increases,  however,  collisional  effects  coupling 
the  four  n  =2  levels  become  increasingly  impor¬ 
tant.  The  effect  is  twofold;  The  density  of  the  ra¬ 
pidly  decaying  singlet  P  level  is  increased  from 
contributions  from  the  other  states,  while  the 
highly  populated  slowly  decaying,  triplet  P  level 
suffers  depopulation  through  contributions  to  other 
»i=2  states.  The  ratio  thus  drops  rapidly  until  a 
density  is  achieved  at  which  the  local-thermody- 
namic -equilibrium  (LTE)  state  distributions  are 
reached  and  the  ratio  becomes  constant  again. 

This  value  is  given  simply  by 

R  =  igAHtriplet)  (io) 

(g.A )  (singlet) 

Using  statistical  weights  and  Einstein  coefficients 
of  9  and  2.54  x  1010  sec"1  for  the  'P  level,  and  3 
and  2.33  x  1013  sec"1  for  the  P  level  gives  us  a 
value  for  R  of  2.69  x  10"3,  in  excellent  agreement 
with  Fig.  10  at  ICP3  ions  cm3.  Although  the  curves 
display  some  temperature  dependence,  an  un¬ 
certainty  of  300  eV  in  electron  temperature  will 
result  in,  at  most,  an  error  in  the  density  de¬ 
termination  of  less  than  a  factor  of  2,  reasonable 
good  accuracy  for  this  plasma  parameter.  Also 
plotted  in  Fig.  10  is  the  intercombination -to-re- 
sonance  line  ratio  as  calculated  bv  Vinogradov 
rt  air3  for  an  aluminum  plasma  at  an  electron 
temperature  of  360  eV.  Although  the  ratio  is  only 
given  up  to  an  ion  density  of  103"  cm"3  in  that  studv, 
the  agreement  between  the  two  calculations  is  very- 
good.  The  small  differences  between  the  values 
obtained  from  that  work  and  ours  are  attributed  to 
the  fact  that  the  collisional-excitation  and  col- 
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FIG.  11.  ls:2p(2P)-ls2p:(2i»/ls:-ls2p(lP>  line 
ratio  v 5  electron  temperature;  optically  thin  t—  i . 

optically  thick  for  30-pm  ( - )  and  500-<im  ( — i 

plasma. 


lisional-ionization-rate  coefficients  used  in  the 
two  studies  have  been  calculated  by  different  me¬ 
thods.13  In  addition,  the  radiative  decay  rate  used 
in  that  work  for  the  Is  2p(  3P)-ls3  transition  differs 
slightly  from  the  one  used  in  our  model.31  How¬ 
ever,  the  similarities  between  the  predictions  of 
the  two  models  is  encouraging. 

The  temperature  dependence  of  the  four  temper¬ 
ature-sensitive  ratios  have  been  plotted  in  Figs. 
11-14  for  an  optically  thin  plasma  (solid  lines)  at  a 
density  of  1020  ions/cm3.  The  satellite-to-re- 
sonance  line  ratio  in  Fig.  11  displays  the  usual 
T^-ig-E'T,  behavior  resulting  from  an  analysis  of 
the  Is2,  ls2p‘P,  ls2p22D  3-level  system.  Ba¬ 
lancing  the  population  of  the  doubly  excited  lith¬ 
iumlike  state  by  electron  capture  with  the  depopu¬ 
lation  of  this  state  by  autoionization  and  spontane¬ 
ous  radiative  decay  yields  the  dielectronic  re¬ 
combination  (DR)  rate32 


FIG.  12.  ls-2‘>/!s2-ls2pilP>  line  ratio  vs  electron 
temperature;  optically  thin  i-i,  optically  thick  for  30- 
um  i - 1  and  300- pm  ( — i  plasma. 


FIG.  13.  ls-2/>/ls:-ls3/>(‘P)  line  ratio  vs  electron 
temperature;  optically  thin  c—  i ,  optically  thick  for  30-pm 
( - 1  and  500-pm  { — )  plasma. 


where  g$  and  are  statistical  weights  of  the  dou¬ 
bly  excited  and  heliumlike  ground  states,  .4r  and 
A a  are  the  radiative  decay  and  autoionization  rates 
for  the  doubly  excited  state,  and  Es  is  the  energy 
separation  between  the  Is2  and  Is 2 p:'D  states.  A 
value  for  the  line  ratio  can  be  obtained  by  di¬ 
viding  aDR  by  the  electron  collisional -excitation 
rate  from  the  Is2  to  the  ls2 plP  level.  An  ap¬ 
proximate  expression  for  this  rate  due  to  Van 
Regemorter33  is 

.Y  =  1.7  x  10-3/(J;A£-'  T,‘  1  '2e  %  .  (12) 

where  f  is  the  oscillator  strength.  is  a  ther¬ 
mally  averaged  Gaunt  factor  and  AE  is  energy 
separation  of  the  transition.  The  division  of  the 
two  rates  yields  the  well-known  Bhalla-Gabriel 


FIG.  if.  l<-!/>/l.<;-ls-3.a<'p>  line  ratio  vs  electron 
temperature;  optically  thin  i-',  optically  th  ck  for 
30-um  ( - t  and  300-pm  ( — t  plasma. 
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formula34  for  this  line  ratio, 


R  =1.212  x  10- 13 


g,/ig.  Te 


f  13) 


This  result  was  first  derived  for  use  in  astro- 
physical  phenomena  and  has  been  a  valuable  diag¬ 
nostic  when  the  two  excited  states  are  populated 
predominantly  from  the  heliumlike  ground  state; 
that  is,  in  low-density  regimes.  As  was  stated 
earlier,  at  higher  densities,  the  ls2p‘P  level  be¬ 
gins  to  be  affected  by  contributions  from  the  Is 2s 
!S  and  Is  states,  and  the  population  density  di¬ 
verges  slightly  from  Its  coronal  value.  For  com¬ 
parison,  we  plotted  the  analytical  result  for  this 
ratio  alongside  the  CRE  result  in  Fig.  11.  One 
curve  employs  a  0.2  value14  for  the  Gaunt  factor 
and  formula  (12)  for  the  excitation  rate,  while  the 
other  uses  the  distorted -wave  result  for  the  ex¬ 
citation  rate  instead  of  formula  (12)  in  obtaining 
the  ratio.  At  an  electron  temperature  of  500  eV, 
our  model  predicts  the  2s  -2p  rate  to  be  50%  of 
the  dominant  ls-2p  rate,  even  at  the  nominal  den¬ 
sity  of  1020  ions/cm3  for  which  Fig.  11  was  ob¬ 
tained.  The  neglect  of  the  complete  collisional 
coupling  at  these  plasma  parameters  would  then 
result  in  a  predicted  radiation  intensity  from  the 
A1XII  La  line  which  is  30%  lower  than  the  CRE 
model,  a  significant  error  in  the  emission  of  the 
most  prominent  line  in  the  x-ray  region  oi  the 
aluminum  spectra.  Hence,  when  employing  this 
ratio  for  accurate  temperature  determination  of 
laboratory  plasmas  at  higher  ( >  1 01'’ '  densities, 
the  analysis  requires  the  complete  n  =2  level 
system. 

The  apparent  agreement  between  the  curve  using 
the  0.2  Gaunt  factor  with  Eq.  (13)  and  the  CRE 
calculations  at  lower  temperatures  is  the  result 
of  the  missing  collisional  couplings  and  the  in¬ 
accurate  rate  coefficient  for  the  ls-2p  transition 
canceling  each  other  in  the  0.2  Gaunt  factor  curve. 

The  satellite-to-resonance  line  ratio  does  pre¬ 
sent  the  plasma  spectroscopist  with  some  in¬ 
herent  difficulties,  however.  The  sensitivity  of 
the  ratio  is  diminished  somewhat  at  higher  tem¬ 
peratures,  the  resonance  line  is  quick  to  be  af¬ 
fected  by  opacity,  and  some  experimental  diffi¬ 
culties  occur  when  trying  to  resolve  the  satellite 
line  completely.  On  the  other  hand,  the  major 
contribution  to  both  upper  levels  responsible  for 
this  ratio  comes  from  the  same  level,  the  helium  - 
like  ground  state.  When  the  two  upper  levels  of 
the  ratio  originate  from  separate  ion  ground 
states,  the  danger  exists  that  these  ions  will  ap¬ 
pear  as  strong  emitters  in  different  spatial  re¬ 
gions  of  the  plasma,  and  the  ratio  will  yield  a 


meaningless  value  for  the  temperature  if  a  tem¬ 
perature  gradient  exists  in  the  plasma.  The  va¬ 
lidity  of  these  model  predictions  would  then  be 
limited,  applying  only  to  experimental  spectra 
which  is  spatially  resolved. 

The  line  ratios  in  Figs.  12-14  are  similar  in 
that  they  compare  a  hydrogenlike  line  intensity 
with  a  heliumlike  line  intensity.  The  basic  fea¬ 
ture  of  all  the  curves  is  the  monotonic  increase 
in  the  ratio  as  the  A1XIII  ground  state  becomes 
more  populated,  relative  to  the  A1XII  ground  state. 
Of  the  three  line  ratios,  the  hydrogenlike  L}  to  the 
heliumlike  L,  in  Fig.  14  is  probably  the  best  one 
in  terms  of  accuracv.  The  two  transitions  which 
make  up  the  ratio  lie  at  very  similar  energies 
t6.055  and  6.13  A),  their  spontaneous  decay  rates 
are  of  the  same  order  (4.77  x  10'-’  and  1.53  x  10u 
sec"1),  opacity  is  not  expected  to  affect  either 
line  in  laboratory -size  plasmas  until  higher  den¬ 
sities  are  achieved,  the  ratio  is  of  order  1  over  a 
range  of  electron  temperatures  i400-800  eV)  com¬ 
monly  encountered  in  laboratory  aluminum  plas¬ 
mas,  and  the  sensitivity  of  the  ratio  in  that  range 
is  very  good. 


C.  Line  ralios  from  optically  thick  plasmas 

In  order  to  account  for  reabsorption  of  emitted 
line  radiation  for  optically  thick  plasmas,  we  at¬ 
tenuate  the  escaping  photon  flux  by  a  phenomeno¬ 
logical  reduction  in  the  radiative  source  function. 
This  reduction  factor,  >,  discussed  earlier,  is 
calculated  to  account  for  the  number  of  ion  scat¬ 
terings  the  photon  undergoes  before  it  diffuses  in 
frequency  to  the  optically  thin  wings  of  the 
broadened  line  profile,  where  it  can  escape.  While 
this  model  can  accommodate  the  radiation  trans¬ 
fer  for  a  natural,  Doppler,  or  Voigt  line  profile, 
the  calculations  presented  here  were  performed 
using  Voigt  line  profiles. 

We  have  plotted  the  line  ratio  behavior  versus 
density  for  a  plasma  with  a  radius  of  100  ,-m  in 
Figs.  6-10  (broken  curves).  Comparing  these 
with  the  optically  thin  curves,  it  is  found  that  the 
thin  and  thick  line  ratios  are  identical  at  .V.  =10‘J 
cm"’,  indicating  that  this  plasma  is  probably 
transparent  to  r  rays  at  these  densities.  Al¬ 
though  this  is,  in  fact,  the  case  here,  the  true 
indication  of  whether  photon  absorption  is  occur¬ 
ring  for  a  particular  line  can  be  found  from  the 
optical  depth,  given  by  Eq.  (2). 

An  optical  depth  at  line  center  c  =  i\)  of  less 
than  0.1  usually  indicates  that  self -absorption  is 
negligible.  However,  larger  optical  depths  can 
occur  accompanied  by  significant  photoexcitation 
in  the  plasma,  yet  no  effect  in  the  line  intensity 
be  apparent.  When  the  plasma  is  radiation  dom- 
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inated,  the  rate  of  stimulated  absorption  may  be 
quite  large  due  to  large  r,  but  if  the  downward 
radiative  processes  (spontaneous  and  stimulated! 
dominate  over  other  collisional  processes  which 
destroy  the  upper  state,  absorbed  photons  will  be 
reemitted  again  and  again.  This  process  will  con¬ 
tinue  until  the  photons  redistribute  in  frequency  to 
the  thin  wings  of  the  line  or  in  space  to  the  edge 
of  the  plasma,  where  they  escape.  Hence,  an  in¬ 
dication  of  whether  photon  absorption  will  affect 
the  line  intensity  is  given  by  the  quenching  para¬ 
meter  for  a  line  whose  upper  state  is  i  and  lower 
state  is  k, 


(H) 


where  Ctj  represents  collisional  rates  and  A,,  is 
the  Einstein  coefficient.  When  is  small,  col¬ 
lisional  quenching  of  the  upper  state  is  not  appre¬ 
ciable,  and  the  emitted  line  intensity  will  be  un¬ 
affected  by  significant  reabsorption.  A  quenching 
parameter  of  near  one,  on  the  other  hand,  indi¬ 
cates  that  any  photon  reabsorbed  is  likely  to  be 
collisionally  quenched,  resulting  in  a  large  re¬ 
duction  in  the  line  intensity.  However,  when  r  is 
large,  the  upper-state  population  will  be  altered 
due  to  the  large  photoexcitation  rate,  regardless 
of  whether  p „  is  large  or  small. 

Comparing  the  general  features  of  the  curves, 
all  of  which  include  a  heliumlike  Lyman-series 
line  as  part  of  the  line  ratio,  a  common  trend  is 
apparent.  The  curves  all  display  a  shoulder  struc¬ 
ture  beginning  at  densities  above  1020  ions/ cm3. 

The  onset  of  this  shoulder  signals  the  point  at 
which  the  AIXII  Lyman  lines  begin  to  be  affected 
by  opacity.  The  La  line  is  usually  the  first  of  the 
lines  considered  here  to  be  attenuated  by  reabsorp¬ 
tion  of  photons  since  its  radiative  decay  rate  is  the 
largest  and  its  lower  state  (which  determines  the 
optical  depth)  is  the  AIXII  ground  state,  the  level 
which  is  the  most  dominant  fractional  density  over 
a  wide  range  of  temperatures  and  densities.  For 
example,  the  intercombination  line  has  an  Ein¬ 
stein  coefficient  which  is  3  orders  of  magnitude 
lower  while  the  satellite  line  and  the  A1  XIII  L a 
line  are  due  to  transitions  which  end  on  different, 
less  populated  levels  over  a  significant  tempera¬ 
ture  range.  The  curves  maximize  and  then  turn 
over  as  the  other  line  comprising  the  ratio  begins 
to  suffer  opacity  reduction.  In  Figs.  6  and  10, 
where  both  the  upper  levels  of  each  ratio  are  popu¬ 
lated  predominantly  from  the  Is2  level  and  gen¬ 
eral  state  abundance  variations  can  be  folded  out, 
the  shoulder  appears  at  a  higher  density  for  the 
high -temperature  curves.  This  is  characteristic 
of  the  r*1'2  dependence  of  the  optical  depth  (see 


Eq.  (3)|;  at  high  temperatures,  the  optical  depth 
will.become  large  only  at  higher  densities.  Hence, 
the  onset  of  opacity  effects  will  occur  at  higher 
densities  for  a  ratio  at  1200  eV  than  one  at  500  eV. 
In  Figs.  7-9,  where  a  hydrogenlike  and  heliumlike 
line  are  compared,  the  shoulder  at  500  eV  is  lar¬ 
ger  than  the  one  at  1200  eV,  emphasizing  the  do¬ 
minance  of  the  AIXII  ion  at  the  lower  temperature. 
However,  the  most  striking  feature  of  many  of  the 
curves,  and,  in  terms  of  diagnostics,  a  rather 
unfortunate  one,  is  the  double -valued  nature  of 
the  line  ratios.  This  aspect  of  these  curves  could 
render  them  useless  as  plasma  diagnostics  unless 
either  the  plasma  density  is  relatively  well  known 
or  several  of  the  diagnostics  are  used  simultane¬ 
ously  in  plasma  parameter  determinations. 

As  an  analytic  check  on  the  model  calculations, 
line -intensity  ratios  in  the  LTE  (high  density) 
limit  with  opacity  effects  included  have  been  cal¬ 
culated  for  the  H(  La)  to  He(  La)  ratio  and  the  IC 
to  He(Ia)  ratio.  To  first-order  approximation 
with  a  Voigt  profile  used  in  the  radiation  transport 
model,  the  line  ratio  can  be  expressed  as 


R  = 


(Tj  \‘'2 

.V..A.AE,  \r, ) 


'15) 


The  main  assumption  here  is  that  the  broadening 
is  nearly  the  same  for  the  two  transitions  forming 
the  ratio.  The  optical  depth  is  proportional  to 
the  radiative  decay  rate,  the  g  factors  of  upper 
and  lower  states,  and  the  lower -state  population 
density;  hence  Eq.  (15)  becomes 


„  _  •’v~2‘X£2  / Aigi:  /gz  .v;\  1  2 
.v,a£,  W,*;  gi-% ! 


'16! 


where  the  prime  refers  to  the  lower  state  of  the 
transition.  In  the  case  of  the  ratio  of  the  two  re¬ 
sonance  lines,  the  Saha  equation10  can  be  used  to 
calculate  the  ratio  of  the  n  =2  and  2 p  lP  levels: 

.\(n=2)  _  2g„-.2  /niek  TAS  ■ 

.V(2pi )  Weg:p  1,  2t n-  / 

A  similar  expression  can  be  used  to  get  the  ratio 
of  the  two  lower  states,  the  AIXII  and  A1XII1 
ground  states.  Inserting  the  appropriate  atomic 
parameters,  the  ratio  is  calculated  to  be  about 
0.26  at  500  eV,  which  compares  quite  well  with  the 
model  prediction  of  0.3  in  Fig.  7,  even  though  the 
levels  considered  here  are  not  in  a  state  of  LTE 
with  each  other  at  a  density  of  1023  ions  cm3.  Den¬ 
sity  criteria  for  LTE  will  be  discussed  in  detail 
in  Sec.  IIID. 

The  LTE  limit  for  the  IC  line-to-resonance  line 
ratio  is  simpler  to  calculate  since  (1)  the  lower 
state  of  both  transitions  is  the  Alxil  ground  state 
and,  (2)  the  strong  collisional  coupling  between 
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the  2 p  singlet  and  triplet  levels  will  equilibrate 
their  populations  to  nearly  a  1:3  ratio  at  even 
moderate  [  ~1023  ions/cm3)  densities.  Hence,  Eq. 
(16)  becomes  (note  that  the  energies  are  nearly 
equal), 


Inserting  the  appropriate  parameters  in  (18)  yields 
a  value  of  0.05  for  the  LTE  ratio  of  the  IC  and 
He(£a)  lines  in  the  optically  thick  case,  while  the 
curves  calculated  by  the  model  tend  to  center 
around  0.04,  again,  quite  favorable  agreement. 

It  is  often  believed  that  the  opacity  effects  on  the 
line  emission  intensities  in  laboratory  plasmas 
might  be  small  enough  to  use  diagnostics  developed 
for  optically  thin  plasmas  if  the  size  of  the  emitting 
region  of  the  laboratory  plasma  was  small  or  the 
densities  were  not  extremely  high.  To  illustrate 
the  opacity  effects  on  laboratory  plasmas,  we  have 
plotted  the  temperature  dependence  on  the  line 
ratios  in  the  optically  thick  case  in  Figs.  11-14  for 
a  plasma  of  1020  ions/cm3  and  radii  of  50  and  500 
inn,  typical  of  a  laser -produced  plasma  and  an 
exploding  wire  plasma.  In  all  ratios  but  the  H(  La) 
to  He(  the  transport  of  the  radiation  has  a 
significant  effect  on  the  values  that  would  be  de¬ 
termined  for  the  electron  temperature.  In  fact, 
the  error  incurred  by  using  optically  thin  calcula¬ 
tions  for  these  plasmas  can  be  as  much  as  200-300 
eV  for  the  500  >im  plasma  and  100-200  eV  for  the 
50-jim  plasma.  Curiously  enough,  the  ratios  may 
also  be  useful  in  determining  the  size  of  the  emit¬ 
ting  region. 

Looking  more  closely  at  the  optically  thick  cal¬ 
culations  plotted  in  Figs.  12-14,  the  results  for 
the  larger -radius  plasma  yield  values  for  the  line 
ratios  which  are  larger  at  low  temperatures  and 
smaller  at  high  temperatures  than  those  predicted 
by  the  optically  thin  calculation.  This  is  a  re¬ 
flection  of  the  shift  in  abundances  from  the  A1XI1 
ion  to  the  A1XIII  ion  as  electron  temperature  in¬ 
creases.  At  500  eV,  the  Is2  level  is  more  popula¬ 
ted  than  the  Is  level;  consequently,  the  optical 
depths  of  the  heliumlike  Lyman  lines  are  greater 
than  those  of  the  hydrogenlike  Lyman  lines,  and 
the  greater  attenuation  of  line  radiation  from 
A1X1I  excited  states  is  manifested  in  the  en¬ 
hancement  of  the  line  ratios.  At  1600  eV,  the 
opposite  is  true  of  the  relative  abundances  of  the 
Is2  and  Is  levels;  therefore,  the  line  ratios  will 
be  decreased  somewhat  compared  to  those  calcu¬ 
lated  in  the  optically  thin  approximation. 

D.  Comparison  vuith  LTE  model 

Since  the  ionization -radiation  model  described 


here  is  characterized  both  by  particle  collisions 
and  photon  processes  in  an  equilibrium  mode,  it 
was  termed  a  collisional-radiative  model.  An¬ 
other  model  frequently  used  to  calculate  the  plasma 
atomic  state  is  the  LTE  model.13  The  models 
differ  in  that  the  LTE  assumption  presupposes 
that  a  radiative  process  is  negligible  when  com¬ 
pared  to  a  competing  ccllisional  process.  The 
CRE  model,  however,  takes  both  processes  into 
account,  explicitly.  An  abundance  curve,  calcula¬ 
ted  using  an  LTE  model,  is  shown  in  Fig.  3  (bro¬ 
ken  curves)  for  a  plasma  at  300  eV,  alongside  the 
CRE  results.  As  would  be  expected,  the  higher 
ionization  stages  approach  the  LTE  limiting  values 
at  higher  densities;  lithiumlike  states  attain  their 
LTE  values  at  an  ion  density  of  2  x  1022  cm"3,  while 
this  occurs  for  the  heliumlike,  hydrogenlike,  and 
stripped  ions  at  densities  of  1023  ,  3  x  1023,  and  6 
x  1023  ions/cm3,  respectively. 

Because  of  its  simplicity,  the  LTE  ionization 
model  has  had  widespread  applications  in  the  plas¬ 
ma  community  for  the  purpose  of  diagnostics  and 
effective  charge  determinations  in  laboratory  plas¬ 
mas.  To  test  the  validity  of  the  application  of  the 
LTE  model,  simple  expressions11  have  been  of¬ 
fered  to  determine  whether  plasma  conditions  are 
within  the  constraints  of  the  LTE  regime.  The 
basis  of  these  expressions  is  the  notion  that  a 
collisional  rate  depopulating  a  level  must  be  a 
factor  of  10  greater  than  a  competing  radiative 
rate  from  this  level  for  the  state  population  to  be 
within  10%  of  its  LTE  value.  We  will  attempt  to 
show,  here,  that  this  can  be  too  strict  a  density 
requirement. 

If  we  consider  the  rates  which  couple  the  Alxrv 
ion  to  the  AIXIII  ground  state,  the  competing  pro¬ 
cesses  are  collisional  (3 -body)  recombination  and 
radiative  (2 -body)  recombination.  The  collisional 
recombination  rate  coefficient  is  given  by,’  in 
cmVsec, 

3  =8.05  x  1 0-28(g,/2gJfx-27;‘  (19) 

while  an  approximate  expression  can  be  used  for 
the  radiative  rate  coefficient,35  in  cm3  'sec, 

a  s  5.2  x  10"l,Z\  (0.43  InA  »0.469\" 1  ’3)  ,  (20' 

while  g  is  the  statistical  weight,  £  is  the  number 
of  outer-shell  electrons,  x  is  the  ionization  en¬ 
ergy,  Z  is  ion  charge,  \  =x/Te,  and  u  and  1  refer 
to  upper  and  lower  states.  The  ratio  3.Va  a  be¬ 
comes  equal  to  10  at  an  electron  density  of  about 
1.7  x  1026  cm"3  yielding  an  ion  density  of  about 
6  x  1025  cm'3.  Hence,  these  two  ground  states 
should  attain  the  population  density  ratio  close  to 
that  predicted  by  the  Saha  equation  for  a  plasma 
at  300  eV  above  an  ion  density  of  6  x  1025  ions,  cm1. 
A  similar  calculation  for  the  hvdrogenlike  and 
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heliumlike  ground  states  yields  an  ion  density  of 
about  9  x  10“  cm'5  at  300  eV. 

A  similar  estimate  can  be  applied  to  excited 
states,  where  the  collisional-deexcitation  rate 
coefficient  must  exceed  the  largest  competing  ra¬ 
diative  decay  rate  by  about  a  factor  of  10  for  the 
upper  and  lower  states  to  be  in  LTE  with  each 
other.  If  we  consider  the  ls-2p  transition  in 
AlXHI  and  employ  the  collisional -excitation  rate 
coefficient  from  our  distorted-wave  calculation, 
the  ion  density  required  for  these  levels  to  be  in 
LTE  with  respect  to  each  other  is  around  1  x  10“ 
cm"3.  A  similar  calculation  with  the  lsi-ls2p‘P 
transition  indicates  that  an  ion  density  of  at  least 
1.5  x  1024  cm'5  is  required  for  these  two  levels  to 
equilibrate  to  their  LTE  values. 

The  values  obtained  from  these  simple  expres¬ 
sions  are  noticeable  higher  than  those  predicted 
by  the  ORE  mcdel  as  seen  in  Fig.  3.  It  would  ap¬ 
pear  that  the  use  of  the  above  expressions  neglects 
the  rather  important  effect  of  the  interplay  between 
the  levels  which  lie  in  energy  between  the  ground 
states.  Since  the  density  at  which  states  equili¬ 
brate  to  their  LTE  values  generally  increases  with 
energy  separation,  many  close-lying  levels  will 
reach  their  LTE  population  ratios  at  much  lower 
densities  than  the  cases  considered  above  which 
are  typified  by  large  energy  differences.  In  an 
apparent  "bootstrapping”  process,  the  coupling 
between  these  levels  and  the  ground  states  is 
enough  to  cause  significant  equilibration  of  the 
ground  states  toward  LTE  at  densities  signifi¬ 
cantly  lower  than  those  predicted  by  a  simple  test 
involving  an  isolated  2-level  system.  Hence,  a 
collisional-radiative  treatment  with  some  excited - 
state  structure  included  is  necessary  to  accurately 
ascertain  at  what  density  and  temperature  a  plasma 
can  adequately  be  described  by  an  LTE  model. 

Although  the  restrictions  on  the  density  required 
for  LTE  are  seen  to  be  eased  somewhat  by  the 
results  of  this  study  with  the  CRE  model,  care 
must  be  taken  when  making  theoretical  predictions 
based  on  diagnostics  from  an  LTE  calculation. 

At  densities  where  K- shell  excited  states  have 
equilibrated  to  their  LTE  population  ratios  but 
the  ground  states  have  not,  the  absolute  popula¬ 
tions  of  the  excited  states  (and  consequently  the 
radiation  emission)  will  be  incorrect,  since  they 
are  determined  relative  to  the  still  non-LTE 
ground -state  densities. 

It  is  well  known  that  photoexcitation  in  an  opaque 
plasma  will  force  it  toward  a  state  of  LTE  at  lower 
densities  than  would  be  predicted  in  the  optically 
thin  approximation.  To  illustrate,  opacity  effects 
were  included  in  the  calculation  for  a  100-nm, 

300 -eV  plasma,  and  the  fractional  distributions 
plotted  in  Fig.  3  (dotted -dashed  lines)  also.  Notice 


that  the  optically  thick  CRE  results  are  always 
between  the  optically  thin  CRE  and  the  LTE  pre¬ 
dictions,  and  that  there  is  a  substantial  change 
in  the  populations  in  a  direction  closer  to  those 
for  LTE  over  a  density  range  of  102o-lO“  ions/ 
cm3.  Notice  also,  however,  that  the  plasma  be¬ 
comes  collision  dominated  above  5x  10“  ions/ cm3 
and  the  CRE  calculation  with  opacity  effects  in¬ 
cluded  does  not  differ  significantly  from  the  CRE 
thin  prediction  of  the  density  at  which  the  plasma 
finally  attains  the  LTE  state. 

IV  DISCUSSION 

This  work  has  investigated  some  of  the  funda¬ 
mental  aspects  of  line  emission  from  a  hot  dense 
aluminum  plasma.  In  particular,  we  have  ex¬ 
amined  the  radiation  signature  by  modeling  the 
plasma  ion  dynamics  using  a  collisional-radiative 
mode!  which  includes  opacity  effects.  Line  ratios 
were  calculated  and  then  compared,  in  two  cases, 
with  older  theoretical  predictions,  and  contrasted 
with  the  result  predicted  when  opacity  effects  are 
taken  into  account.  We  have  shown  that  what  may 
appear  to  be  a  well-behaved  and  understood  diag¬ 
nostic  in  an  optically  thin  approximation  may ,  in 
fact,  be  a  rather  complex  double -valued  function 
dependent  on  both  plasma  temperature  and  density 
when  the  optical  depth  of  a  photon  through  the  plas¬ 
ma  becomes  large.  An  important  result  of  this 
work  which  we  would  like  to  stress  is  not  that 
these  line  ratios  may  become  ineffective  diag¬ 
nostic  tools  when  they  are  taken  from  a  spectra 
emitted  from  an  optically  dense  plasma.  Rather, 
it  now  becomes  apparent  that  a  comprehensive 
treatment  of  the  radiation  diagnostics  is  required. 
Close  agreement  between  a  number  of  separate 
temperature  or  density  diagnostics  can  resolve  the 
dilemma  of  double  valuedness.  The  comparison  of 
parameters  determined  from  line  ratios  and  those 
obtained  from  other  methods  such  as  series  mer¬ 
ging,36  Stark  profiles.37  or  slope  of  the  continuum 
can  also  be  an  effective  technique.  However,  the 
entire  spectrum  must  be  analyzed  in  detail  to  pro¬ 
perly  make  use  of  the  plasma  diagnostic  techniques 
available. 

The  results  of  this  model  have  also  been  com¬ 
pared  with  the  LTE  result  and  it  has  been  seen 
that  the  plasma  cannot  be  fully  described  by  an 
LTE  model  until  rather  large  densities  have  been 
obtained.  Although,  in  an  opaque  plasma,  excited 
states  may  approximate  their  LTE  populations  at 
a  significantly  lower  density  than  in  an  optically 
thin  plasma,  due  to  optical  pumping  through  re- 
absorption  of  photons,  the  ground  states  will  still 
maintain  populations  between  the  LTE  and  coronal 
limits  at  this  lower  density,  only  slightly  affected 
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by  the  opacities  considered  here.  Thus,  though  the 
excited -state  populations  may  be  correctly  deter¬ 
mined  relative  to  the  next-higher  ground-state 
density,  line  ratios  calculated  from  lines  which 
originate  in  excited  states  from  different  ioniza¬ 
tion  stages  will  be  in  error  since  the  ground  states 
are  not  accurately  given  by  their  LTE-determined 
values.  It  is  also  worth  noting  that  there  is  a  disa¬ 
greement  between  abundance  curves  quoted  here 
and  those  shown  in  Ref.  9,  which  depicts  highly 
overpopulated  Aim  and  A1XI  ions.  This  result 
appears  unphysical  in  light  of  the  previous  dis¬ 
cussion  of  electron  shell  configurations  and  their 
effect  on  the  ionization  dynamics.  The  lithium¬ 
like  ion  with  an  ionization  potential  of  4-12  eV 
should  ionize  easily  to  the  closed -shell  A1XII  ion 
and  should  never  constitute  a  major  partial  density 
fraction  over  a  wide  range  of  temperatures  at  an 
ion  density  of  only  1022  cm"’  as  is  shown  in  Ref.  9. 
This  is  in  agreement  with  work  done  previously 
by  other  authors38- 39  on  aluminum  as  well  as  other 
materials  of  nearly  the  same  atomic  number.  In 
addition,  the  study  in  Ref.  9  includes  a  graph 
showing  that  their  CRE  model  is  in  agreement 
with  the  LTE  result  for  the  Alxrv  population  at  an 
ion  density  of  about  5  x  1022  cm'3,  a  factor  of  10 
lower  than  in  Fig.  3.  It  is  our  opinion  that  this 
disagreement  is  due  to  their  collisional-ionization 
rate  from  A1XIU .  An  approximate  form  for  the 
rate  is  given  in  Ref.  6,  in  cm3/sec, 

S*£x-*T,lV<'*r«  ,  121) 

where  {  is  the  electron  number  in  the  outermost 
shell.  This  would  predict  a  rate  for  AlXU  which 
is  about  twice  as  large  as  that  for  A1XIII  (ioniza¬ 
tion  potentials  are  within  l(Fo  of  each  other).  This 
is  in  disagreement  with  Fig.  5  of  Ref.  9  which 
shows  the  Alxm  rate  to  be  about  two  orders  of 
magnitude  greater  than  the  A1X1I  rate  at  T€  =300 
eV.  The  AlXU  collisional-ionization  rate  used 
in  this  study  is  in  good  agreement  with  the  one 
quoted  in  Ref.  9,  hence  it  is  the  questionably  large 
AlXIII  ionization  rate  which  is  the  cause  of  their 
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Some  discussion  of  Che  approximations  made  in 
constructing  this  model  is  warranted  in  order  to 
properly  assess  the  significance  of  the  analysis 
presented  here.  The  equations  governing  the  phv- 
sics  described  above  are  solved  numerically  by  a 
zero-dimension  (single-zone;  computer  model; 
this  allows  only  for  treatment  of  plasmas  which  are 
homogeneous  in  temperature  and  density  at  the 
time  of  the  radiation  pulse.  In  addition,  the  rate 
equations  are  solved  time-mdependently;  hence, 
the  assumption  is  made  that  the  width  of  the  ra¬ 
diation  pulse  is  small  compared  to  the  hydro - 
dynamic  time  scale  of  the  plasma.  The  atomic 
level  structure  used  was  selected  in  order  to  de¬ 
scribe  the  most  important  K- series  emission  lines, 
based  on  analysis  of  experimental  spectra  and  the 
states  contributing  most  significantly  to  its  forma¬ 
tion,  either  quantitatively  or  qualitatively  (within 
the  temperature  range  of  interest  of  this  work). 
Thus,  we  have  retained  some  structure  which  is 
less  important  energetically,  for  diagnostic  pur¬ 
poses.  We  have  not  included  detailed  structure 
in  most  of  the  /.-shell  ions,  leaving  that  for  later 
investigations.  Finally,  several  atomic  processes 
have  been  omitted  such  as  multiphoton  excitation, 
photoionization,  and  inner-shell  collisions.  The 
model  retains  the  flexibility  to  include  these  pro¬ 
cesses  in  future  work,  but  their  contribution  to¬ 
ward  the  determination  of  the  level  populations  is 
considered  of  less  importance  in  plasmas  studied 
in  this  paper. 
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abstract 


Radiation  amission  over  a  vide  range  (5-5000  %.)  of  wavelengths  in  a 
high  density  aluminum  plasma  has  been  studied  theoretically.  The  atomic 
model  includes  explicit  treatment  of  a  large  number  of  states  and 
transitions  characterizing  the  line  emission  from  the  A1  III  -  A1  XIII 
ions.  The  state  populations  are  calculated  using  a  set  of  atomic  rate 
equations  describing  a  plasma  in  collisional-radiative  equilibrium.  The 
ionization  state  of  the  plasma  is  calculated  self-consistently  with 
line  and  continuum  radiation  transport  using  a  1-dimensional  multi-cell 
structure  for  the  plasma.  A  new  radiation  transport  scheme  has  been 
implemented  in  this  model  based  on  probability  of  escape  for  both  line 
and  continuum  photons.  The  total  L  and  K-shell  radiation  emission  is 
described  as  a  function  of  plasma  temperature ,  density,  and  source  size, 
and  is  compared  to  predictions  made  from  corona  and  local  thermodynamic 
equilibrium  models.  It  is  shown  that  over  a  wide  range  of  plasma  para¬ 
meters,  neither  of  these  two  limiting  models  accurately  describes  the 
plasma  emission  characteristics.  Several  line  intensity  ratios  from 
L-shell  lines  are  studied  and  shown  to  be  useful  diagnostic  indicators 
of  the  plasma  temperature  and  density.  Finally,  theoretical  spectra  are 
calculated  and  discussed  in  the  context  of  their  comparability  with 
experimental  spectra  in  order  to  make  contact  with  real  plasma  radiation 
signatures  and  exploit  the  predictive  capabilities  of  the  model. 


■ -«**-** 


I.  Introduction 

The  x-ray  amission  from  high  density  aluminum  plasmas  has  received  much 
attention  in  recent  years  due  to  interest  in  results  from  exploding  wire 
experiments  ,  laser  plasmas ,  and  particle  beam— foil  interactions^^. 
Theoretical  analysis  of  the  radiation  emission  from  these  dense  plasmas  is 
useful  for  understanding  both  the  plasma  energetics  and  determination  of  plasma 
parameters  from  spectral  diagnostics.  Several  models  have  been  developed  in  an 
effort  to  understand  various  aspects  of  aluminum  plasma  emissions,  concen¬ 
trating  on  K-shell  ions ^  some  L -shell  and  M-shell^"^  ions. 

In  this  study  we  include  excited  states  from  all  ions  of  the  K  and  L  shells  of 
aluminum  as  well  as  A1  III  in  the  M  shell.  In  addition,  the  ionization  dynamics 
and  radiation  transport  calculations  are  solved  self-consistently  for  all  in¬ 
cluded  level  populations  producing  line  and  continuum  photons  in  a  one-dimensional 
plasma  model.  The  multicell  aspect  of  this  model  also  allows  us  to  treat 
stationary  plasmas  in  which  there  exists  temperature  and  density  gradients 
characteristic  of  imploding  wire  array  plasmas. 

This  work  will  be  presented  in  four  sections.  First,  a  detailed  description 
of  the  level  structure  and  rate  coefficients  employed  in  this  model  is  presented 
along  with  a  brief  description  of  the  ionization-radiation  calculation  as  well 
as  the  radiation  transport.  Second,  the  gross  characteristics  of  both  the  line 
and  continuum  radiation  are  studied  as  they  vary  with  the  plasma  electron  temper¬ 
ature,  ion  density  and  plasma  radius.  The  results  are  then  compared  with  those 
obtained  from  local  thermodynamic  equilibrium  (LTE)  and  corona  models  and  define 
the  plasma  conditions  necessary  for  the  valid  application  of  these  models.  Third, 
we  explore  the  details  of  the  bound-bound  radiation  spectrum  to  describe  various 
line  intensity  ratios  from  the  L  shell  ions  which  can  effectively  be  applied  to 
diagnose  plasma  parameters  in  temperature  regimes  lower  than  that  in  which  K-shell 
ions  are  abundant.  Finally,  in  the  fourth  section  we  present  theoretically 
generated  emission  spectrum  and  discuss,  in  detail,  their  implications  with 
respect  to  plasma  diagnostics.  Although  these  calculations  have  been  done  for 
an  aluminum  plasma,  the  interplay  between  particle  collisions  and  photon  transport 
in  determining  level  populations  and  radiation  emission  are  expected  to  be  typical 
of  many  dense  plasmas.  Therefore,  we  have  attempted  to  be  as  general  as  possible 


0  3 


in  our  analysis  of  the  results,  realizing  that  many  of  the  trends  occurring 
in  this  study  will  be  typical  of  L-shell  emission  from  other  materials.  Some 
features  of  the  data  will  nevertheless  be  strongly  deoendent  on  the  detailed 
rates  and  structure  representative  of  aluminum  in  particular,  and  ve  have 
made  an  effort  to  point  out  these  occurrences  when  not  obvious. 

II.  ATOMIC  MODEL 

A.  Level  Structure 

Integral  to  any  theoretical  prediction  of  the  radiation  emission 
from  a  plasma  which  spans  a  significant  temperature  range  in  its  evolution  is 
the  detailed  system  of  states  included  in  the  atomic  model.  In  a  previous 
work  on  K— shell  emission*'^,  we  chose  to  include  several  higher  lying  states 
whose  emission  was  diagnostically  useful  but  were  not  energetically  signifi¬ 
cant.  In  this  study,  it  was  necessary  to  include  many  states  for  energetic 
purposes  alone.  3ecause  of  the  smaller  energy  differences  between  ground  and 
excited  levels  in  L  shell  ions  as  compared  to  those  of  the  K-shell,  excited 
level  populations  constitute  a  non-negligible  fraction  of  the  total  population 
at  ion  densities  typical  of  exploding  wire  plasmas  and  one  would  like  to  include 
many  levels  in  each  ion  to  insure  completeness  of  the  radiative  description.  On 
the  other  hand,  several  criteria  restrict  the  number  of  levels  one  can  practi¬ 
cally  include  In  a  model:  (i)  the  availability  of  data  for  higher  Rydberg 
members  of  L-shell  aluminum  ions,  (ii)  the  unwieldiness  of  solving,  self-con- 
sistently  with  radiation  transport,  a  set  of  rate  equations  for  several  hundred 
levels,  (iii)  tine  restrictions  on  computer  runs  to  solve  this  large  set  of 
equations  for  a  meaningfully  extensive  set  of  parameters.  With  these  restrictions 

in  mind,  the  number  of  levels  to  include  per  ion  was  decided  upon  after  scanning 

(13) 

several  sources  of  XUV  spectral  data  and  identifying  the  strongest  lines. 

These  spectra  showed  that  in  =  0  transitions  of  the  form  2s^2pn  -  2s2prr!"^  and 
n+2 

2s2?'  -  2?  were  very  evident  in  many  of  the  shoes,  despite  the  high  density 

of  the  plasmas.  In  addition,  very  few  lines  were  identified  above  princiole 
quantum  number  3  for  all  but  the  alkali  ions,  A1  XI  and  A1  III.  We  chose, 
therefore,  to  include,  in  addition  to  the  ground  states,  the  n  =  2  and  n  **  3 
levels  in  all  ions  in  as  complete  a  manner  as  possible,  as  well  as  a  few  higher 
n  states  in  the  two  alkali  ions.  The  complete  level  structure  employed  in  this 
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study  is  shown  Table  I  along  with  the  state 

The  energies  o:  the  states  were  taken  from 

(14) . . 

sccnar  anc  cna  caoj.es  oi  Marten  ana  _a^ 


bound  transitions  were  taken  from  reference 

(16) 

Palumbo 


energy  measured  from  ground) . 

the  level  diagrams  of  Bashkin  and 
0-5  i 

uoas  ;  wave j.en gens  or  oounc- 
14  and  the  tables  of  Kelly  and 


In  a  previous  study  on  the  K-shell  emission it  was  found  that  at 
high  densities,  many  of  the  individual  l  states  of  a  given  n  level  may  be 
in  Boltzmann  equilibrium  with  each  other  due  to  heavy  particle  coll:'  sions . 
Specifically,  if  the  energy  separation  between  level  nl  and  n V  is  sma^l  com¬ 
pared  to  the  plasma  thermal  energy,  icn-ion  collisions  will  cause  the  levels 
to  equilibrate  at  densities  lower  than  electron-ion  collisons  would.  Since  the 
primary  purpose  of  this  investigation  is  to  study  high  density  plasma  radiation 
(typically,  ion  densities  greater  than  10  cm  ),  the  assumption  has  been 

made  that  many  of  the  close-lying  states  have  equilibrated  and  the  ratio  of 
their  populations  is  given  by  the  ratio  of  their  statistical  weights.  The 
individual  n l  levels  are  then  grouped  as  a  single  "degenerate"  n  level  with 
the  appropriate  statistical  weight.  In  addition,  any  rata  populating  this  n 
level  is  an  average  (over  initial  states)  of  each  rate  populating  the  indi¬ 
vidual  ni  states.  This  method  of  treating  levels  of  similar  electron  configu¬ 
ration  and  separated  by  only  a  few  electron  volts  in  energy  allows  the  model 
to  carry  a  large  number  of  excited  levels  while  reducing  computation  time  that 
would  have  been  used  calculating  a  known  relationship  between  many  of  the 
levels.  The  brackets  in  Table  I  indicate  the  levels  which  were  appropriately 
combined  in  our  model. 


B.  Rate  Coefficients 

With  the  exception  of  electron  collisional  excitation  rate  coeffi¬ 
cients  and  spontaneous  decay  rates,  the  methods  of  calculating  the  various 
rate  coefficients  used  in  this  study  have  been  described  in  detail  in  previous 
papers^’ ^ \and  only  brief  reference  will  be  made  to  these.  Every  state  is 
coupled  to  the  next  (energetically)  highest  ground  state  by  collisional  ioni¬ 
zation,  photo-ionization,  collisional  recombination  and  radiative  recombination. 

( 1  Q) 

The  collisional  ionization  races  are  calculated  by  Seaton's  prescription  "  ,  and 
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Che  photoionization  rates  are  calculated  in  the  hydrogenic  approximation 

(2i  ) 

using  affective  free-bound  Gaunt  factors  .  Collisional  and  radiative  raccn- 


Lmaci.cn  rates  are  then  calculated  as 


addition,  adjacent  ground  scaces  are  coupled  by  dielectronic  recombination  ~ 

found  by  summing  capture  rates  over  a  manifold  of  states  above  the  ionization 

limit  which  then  decay  via  cascade  to  the  ground  state. 

Excited  levels  of  a  given  ion  are  coupled  to  other  selected  excited 

levels  and  the  ground  state  by  electron collisional  excitation  and  de-excitation, 

stimulated  amission  and  absorption,  and  spontaneous  radiative  decay.  In  the 

X-sheil  ions,  A1  XXI  and  Al  XXII,  and  in  A1  XI,  the  collisional  couplings  axe 

complece,  including  forbidden  and  spin-flip  transitions  as  well  as  those  which 

(22) 

are  dipole  allowed.  For  Al  XIII,  the  Coulomb -Bom  approximation  was  used. 


while  a  distorted  wave  calculation  with  exchange 


was  used  to  calculate  the 


coefficients  for  Al  XI  and  Al  XII.  Comparisons  of  these  two  methods  have  been 
made  for  several  transitions  in  Al  XIII^^\  and  excellent  agreement  was  obtained. 

The  prohibitively  large  number  of  levels  to  be  collisionally  coupled 
in  the  remaining  ions  precluded  the  indiscriminate  use  of  the  distorted  wave 
calculation  for  these  transitions.  Instead,  a  calculation  employing  the  semi- 


classical  impact  parameter  method  with  thermally  averaged  Gaunt  factors ' 
was  used  to  couple  levels  betwean  which  electric  dipole  transitions  are  allowed, 
while  forbidden  transitions  were  neglected.  When  dipole  oscillator  strengths, 
required  by  the  impact  parameter  code,  were  not  available,  a  3ates-Damgaar d 
calculation v  was  used  to  generate  them. 

Many  of  the  ions  contain  excited  levels  from  several  spin  systems  and 
it  was  necessary  to  collisionally  couple  these  systems  (usually  between  the 
lowest  lying  excited  state  in  the  system  and  the  ground  state);  the  distorter 
wave  calculation  was  used  for  these  transitions  to  take  account  of  the  spin 
If  these  couplings  are  not  included,  the  state  populations  in  the  spin  sv-- 
not  coupled  to  the  lower  ground  state  are  found  to  float  somewhat  a:  1 
trcn  densities,  being  coupled  only  to  the  ground  state  above 
ionization  and  radiative  recombination.  While  collisional  ior.—  . 
the  dominant  de-populatir.g  mechanism  for  these  states  at  :cr~n-_ 
collisional  excitation  from  the  ground  state  dominates  ?  ve r  - 
bination  as  the  main  populating  process;  hence,  if  the-’ 
rates  are  not  included,  incorrect  populations  tor  tr¬ 
oy  the  model . 


averaged  Gaunt  factors 
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Dipole  oscillator  strengths  and  spontaneous  radiative' decay  rates  for 

the  L-shell  ions  were  extracted  from  several  sources.  The  very  complete  tabu— 

O  ,  _  frA  \ 

in cion  by  lindgard  and  Nielsen  was  the  source  of  these  races  for  1 ichiun— 

like  Al  XI  and  sodlualike  Al  III.  Many  of  the  oscillator  strengths  describing 

the  transitions  2s22pn  -  2s2pn+1  and  2s2pQrfl  -  2pQ+‘2  for  beryllium-like  Al  X 

through  oxygenlike  Al  VI  were  found  in  the  calculations  by  Fawcett  ^  ^ .  The 

following  references  were  also  consulted  as  sources  for  oscillator  strengths: 

Martin  and  Simons®0^  for  Al  X,  Shamey®1  ^  and  Weiss®2  *  for  Al  IX,  McIntyre 

et  al.®3^for  Al  VIII,  Buchet  and  Buchet-Poulizac®4  ^  for  Al  V  -  VII, 

Kern  ah  an  et  al.  ^  ^  for  Al  III  -  VII,  Mewe®^^for  Al  V,  Shorer^  ^  and 

G8  ) 

Loginov  and  Gruzdev'’  for  Al  IV.  Finally,  the  tables  of  transition  probabilities 
(39  ) 

by  Wiese  et  al.  were  used  for  many  of  the  L-shell  f-values  not  located  in  the 
above  sources.  As  stated  earlier,  all  oscillator  strengths  not  contained  in 
the  references  were  generated  from  a  Bates-Dangraad  calculation. 

A  brief  discussion  at  this  point  about  the  effect  of  neglecting  for¬ 
bidden  transitions  in  the  L-shell  ions  is  warranted.  At  low  density,  this  is  a 
good  approximation  if  the  radiation  emission  is  the  quantity  one  is  most  con¬ 
cerned  with.  At  these  densities ,  the  plasma  is  in  the  coronal  regime  and  excited 
states  are  populated  predominantly  by  electron  collision  from  the  ground  state. 
Consequently,  only  low-lying  n  -  2  levels  will  have  significant  populations  and 
the  bulk  of  the  line  radiation  will  be  due  to  transitions  from  these  states. 

Since  almost  all  of  the  transitions  involving  a  2s-2p  excitation  are  dipole- 
allowed,  the  n  »  2  populations  and  hence  the  radiation  will  be  calculated 
accurately.  Since  many  of  the  n  _>  3  levels  are  forbidden  in  their  coupling  to 
Che  ground  state,  their  populations  may  not  be  as  accurate,  but  the  radiation 
from  these  states  is  negligible  at  low  densities.  At  high  densities,  collisional 
ionization  and  recombination  between  excited  states  and  the  next  highest  ground 
state  dominate  the  redistribution  of  the  populations  among  the  excited  states; 
since  these  rates  are  included  in  the  model,  all  population  densities  are 
accurately  determined  in  the  LIE  regime.  The  collisional -radiative  regime  between 

corona  and  LIE  is  where  collisional  excitation  processes  are  most  important  in 
determining  the  populations  of  states  with  n  _>  3.  Although  forbidden  tran¬ 
sitions  will  definitely  affect  the  populations  of  and  radiation  from  these 
levels,  we  place  the  maximum  error  in  populations  due  to  the  neglect  of  these 


rates  at  a  maximum  of  a  factor  of  2  and  an  uncertainty  of  at  most  SOX 
in.  the  total  radiation  emission  (based  on  our  previous  work  vita  coliisional. 
couplings  in  the  aluminum. Krshell) .  The  inclusion  of  forbidden  transitions 
and  their  effect  en  the  L-shell  line  emission  will  be  the  subject  of  a  future 
study. 


C.  Ionization  Dynamics  and  Radiation  Transport 

The  scheme  for  solving  the  ionization  dynamics  of  an  equilibrium 
plasma  has  been  described  in  detail  in  earlier  work^  ^ .  Briefly,  an  atomic 
rate  equation  of  the  form 


-  E 

j 


wij  »i 


CD 


can  be  written  for  each  ground  and  excited  state  included  in  the  model,  where 
is  a  rate  for  populating  state  i  and  W  is  a  rate  for  depopulating  state 
i.  In  this  study,  the  approximation  is  made  that  the  plasma  radiation  pulse 
is  short  compared  to  characteristic  hydrodynamic  time  scales, and  hence,  state 
populations  are  constant  during  the  duration  of  the  pulse.  This  allows  for 
a  steady-state  solution  and  the  time  derivative  can  be  set  to  aero  in 
equation  (1).  In  addition,  no  hydrodynamic  equations  are  as  yet  solved  together 
with  the  rate  equations;  thus,  the  solution  to  the  set  of  equations  (1)  re¬ 
presents  a  plasma  frozen  in  space  and  time  at,  for  example,  the  point  of  peak 
plasma  emission. 

The  radiation  transport  employed  here  represents  an  improvement  on  our 
phenomenological  frequency  diffusion  technique  reported  on  earlier  C6»l?)^ 
based  on  a  calculation  of  the  probability  of  escape  of  photons  born  in  the 
plasma.  The  probability  that  photons  born  in  a  given  zone  of  the  plasma  and 
re-absorbed  in  a  different  zone  (or  allowed  to  escape  to  the  detector)  is 
calculated,  taking  into  account  both  stimulated  emission  and  stimulated  absorption 
processes.  These  probabilities  are  calculated  assuming  that  the  line  emission 
profile  is  a  Voigt  and,  therefore,  both  Doppler  and  coliisional  broadening 
contribute  to  the  line  width. 
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The  probability  of  escape  for  a  Voigt  line  profile  is  defined  by 
40  ,  . 

aummer  and  is  given  oy. 


_ 

Pe 


ft! 


where  tq  is  the  line-center  optical  depth  and  a^  is  the  Voigt  damping 


parameter  defined  by. 


£  +  l  X..  +  £  A.,  +  £  X., 

k  i*  k  ik  k  ik  k  ik 

4  ir  i  v. 


Here,  A  and  X  represent  all  radiative  and  eollisional  rates  depopulating  the 
upper  level,  i,  and  the  lower  level,  j,  of  the  transition,  while  Av^  is  the 
Doppler  width  of  the  line. 


r-i? 

IvJ  v< 


where  is  ion  temperature,  M  is  ion  mass,  and  ^  is  the  photon  frequency 
at  line  center.  Using  the  expression  in  equation  (2) ,  the  probability  that 
a  photon  emitted  in  cell  m  will  traverse,  unabsorbed,  to  a  cell  boundary  can 
be  calculated.  Thus,  using  the  difference  of  such  probabilities,  one  can 
determine  the  probability  that  a  photon  emitted  in  zone  m  will  be  absorbed 
in  "one  m* .  By  calculating  zone-to-zone  coupling  coefficients,  C(m,mr), 
in  this  way  for  each  emission  line,  one  computes  the  probability  of  escape 
from  zone  m  by, 


P  (m) 
z 


1  -  C(m,m) 


where  C(m,a)  is  the  coupling  coefficient  relating  the  photon  absorption 
probability  before  reaching  the  first  zone  boundary.  The  photon  pumping  rate 
in  zone  m  due  to  photons  emitted  in  ocher  zones  and  absorbed  in  zone  m  is 
then  given  by, 


(6) 


•k 

where  5  is  che  density  of  the  upper  state  of  i.uc  transition,  A  is  the  radia-_ 
eive  decay  rate,  A V  is  the  zonal  volume,  and  the  coupling  coefficients, 
C(m,m’),  contain  terms  accounting  for  both  stimulated  absorption  (photo¬ 
excitation)  and  stimulated  emission.  The  quantities  calculated  in  eqns.  (5) 
and  (6)  are  employed  in  the  atomic  rate  equations  since  the  reabsorption  of 
photons  can  significantly  alter  the  ionization  state  and  excited  state  distri¬ 
butions  in  the  plasma^.  Finally,  we  define  the  probability  of  escape  of  a 
line  photon  from  zone  m  to  the  spectrometer. 


and  the  line  intensity  that  exists  the  plasma. 


i  .  2 


(m)  A  P  (m)  hv  AV(m)  . 
esc  o 


In  a  time-dependent  calculation,  1^  would  be  integrated  in  time  to  describe 

the  temporal  evolution  of  the  spectral  line  as  the  state  densities  and  opacity 

of  the  plasma  change.  A  more  detailed  discussion  of  the  method  of  calculating 

the  coupling  coefficients  and  the  probability  of  escape  transport  technique 

41 

rzm,  be  found  in  the  work  of  Apruzese  et  al  .  In  addition  to  transporting  line 
photons,  we  have  developed  a  method  for  transporting  bound-free  continuum 
photons  as  well.  Using  the  photoionization  cross  sections  calculated  by  a 
method  described  earlier,  a  probability  of  escape  for  photons  originating  from 
radiative  recombinations  is  calculated  based  on  the  recombination  edge  profile 
and  implemented  in  the  rate  equations  in  a  way  similar  to  that  for  the  line 
photons.  A  detailed  description  of  this  method  will  be  described  in  a 
following  work. 
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III.  RADIATION  POWER  IN  PLASMAS 


The  main  emphasis  of  this  study  is  to  quantify  both  the  gross  and  detailed 
features  of  the  radiation  emission  from  the  L-shell  iocs  of  dense  a  t  ’m-inim 
plasmas.  Previously,  it  was  necessary  to  characterize  the  gross  nature  of  the 
emission  by  either  a  corona  model  in  the  low  density  limit,  or  an  LTE  model 
in  the  high  density  limit,  in  this  section,  we  present  collislonal-radiatlve 
equilibrium  (CRE)  results  from  a  study  of  the  line  and  continuum  radiation  as 
they  vary  with  plasma  temperature,  density  and  size.  In  addition,  comparisons 
are  made  between  predictions  of  this  CRE  calculation  and  the  other  two  limiting 
models  to  establish  regions  of  validity  for  their  application. 

In  order  to  assess  the  collisional  properties  of  the  plasma  and  allow 
for  a  more  equitable  comparison  between  models,  a  study  of  the  total  line 
emission  as  a  function  of  temperature  and  ion  density  was  conducted  with  the 
initial  assumption  that  the  plasma  was  optically  thin.  The  results  are  shown 
in  Fig.  1,  where  the  total  bound-bound  emission  power  coefficient  (Power  den¬ 
sity  divided  by  electron  density  divided  by  ion  density)  is  plotted  as  a 

11  23  —3 

function  of  electron  temperature  for  various  ion  densities  from  10  -  10  cm  . 

Plotting  the  power  coefficient  allows  one  to  easily  discern  the  departure  from 
corona  equilibrium,  since  that  model  predicts  a  coefficient  which  is  independent 
of  density.  The  radiation  from  L-shell  ions  is  seen  in  Fig.  1  to  peak  about 
60-120  eV  while  the  K-shell  emission  maximizes  at  1  keV.  The  curves  for  10 ^ 
and  10^  cm  ^  are  Identical,  indicating  that  the  plasma  can  truly  be  described 
by  a  corona  model  at  those  densities.  Deviations  from  coronal  equilibrium 
become  noticeable  above  10^  cm  ^  for  L-shell  radiation  and  above  10^  cm  ^ 
for  K-shell  radiation.  At  low  densities,  the  maximum  radiative  emission  at  the 
L-shell  exceeds  that  of  the  K-shell  by  an  order  of  magnitude,  while  at  high 
densities,  the  K-shell  radiation  is  much  stronger.  This  is  due  to  the  increase 
in  electron  collisional  quenching  of  L-shell  photons  as  the  density  is  in¬ 
creased  and  the  L-shell  ion  states  tend  toward  LTE.  In  fact,  at  very  high 
density,  the  L-shell  peak  disappears  completely,  leaving  a  functional  dependence 
suggestive  of  the  power  radiated  by  a  black-body,  viz 

P  -  oT4  (9) 

where  c  is  the  Stefan-Boltzmann  constant  and  T  is  plasma  temperature.  The 
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deviation  from  corona  equilibrium  for  the  K-shell  first  appears  at  about  10 
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ions/ cm  ;  this  occurs  at  a  much  higher  density  for  the  K-shell  because  typical 
energy  separations  between  ground  and  excited  states  in  the  K-shell  are  much 
larger, and  her.ce,  collisional  de-excitatian  rates  are  much  smaller  than  chose 
for  L-sheli  ions.  (The  corona  approximation  is  no  longer  valid  when 
collisonal  de-excitation  becomes  comparable  to  radiative  decay  in  depopulating 
excited  states) . 

The  curves  shown  in  Fig.  1  are,  of  course,  due  to  the  manifold  of  L  and 
K-shell  ions  emitting  line  radiation  which  is  peaked  at  temperatures  corre¬ 
sponding  to  roughly  the  peak  in  their  ground  state  abundances.  The  Individual 
ionic  contributions  to  the  total  line  emission  are  shown  in  Fig.  2,  where 

radiated  power  density  vs.  T  (electron  temperature)  for  an  optically  thin 
19  3  e 

plasma  at  10  ions/cm  is  plotted.  Notice  that  the  boronlike  A1  IX  peak 

falls  somewhat  below  the  A1  VIII  peak  disrupting  what  was  a  monotonic  Increase 

in  peak  intensities  with  increasing  ionic  charge.  This  could  be  due  to:  some 

inherent  atomic  structure  differences,  such  as  a  reduction  in  intensity  as  a 

closed  subshell  is  attained;  a  possible  error  in  an  oscillator  strength;  or 

some  level  structure  being  left  out  of  the  model.  Since  the  peak  would  need 

an  increase  of  about  30Z  to  bring  it  in  line  with  the  ionic  trend,  it  is  possible 

that  the  radiation  from  some  n  -  3  states  not  included  in  the  model  could  account 

for  the  difference.  According  to  model  predictions,  the  strongest  radiating 

lines  at  70  eV  (A1  IX  peak  temperature)  in  A1  IX  are  2s22p(2P)  -  2s23s,  -  3p, 

-3d  (doublet),  2s2p2(4?)  -  2s2p3d(4P,  4D) ,  and  2s22p(.2P)  -  2s2p2(2S,  2P,  *n. 

in  that  order.  The  2s2p3s,  -  3p,  (quartet)  states  were  not  included,  but  it's 

possible  that  they  could  have  accounted  for  a  significant  portion  of  the  A1  IX 

2  4 

emission  due  to  their  transitions  to  the  2s2p  (  ?)  level. 

It  was  stated  earlier  that  in  the  coronal  regime  the  bulk  of  the  plasma 
line  radiation  emanating  from  the  L-shell  ions  is  due  to  An  -  0  transitions 
within  the  n  »  2  states.  To  illustrate  how  that  changes  as  density  Increases, 

the  individual  line  power  densities  from  A1  VIII  (carbonlike)  are  shown  vs.  T 

19  -3  * 

for  an  optically  thin  plasma  at  10  cm  in  Fig.  3.  The  transitions  in  A1  VIII 

shown  in  this  figure  are  listed  in  Table  II  along  with  their  approximate  wave¬ 
lengths.  Note  that  while  n  ■  2  transitions  still  dominate,  a  large  portion  of 
the  radiation  now  comes  from  decays  from  n  ■  3  levels.  The  peak  intensities 
for  these  transitions  occur  at  slightly  higher  temperatures  than  the  n  *  2 
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transitions  due  to  thslr  larger  excitation  energies.  The  extjent  to  efcich 
the  higher  Rydberg  states  contribute  to  the  total  ionic  line  emission  at 
high  plasma  dan3icias  ran  seriously  affect  the  validity  of  the  atomic  so dal.  . 

If  not  enough  level  structure  is  included  in  the  model*  the  theoretical  calcu¬ 
lation  may  seriously  underestimate  the  gross  radiation,  features.  As  an  example, 

one  would,  probably  assume  that  n  ■  4  and  n  ■  5  levels  would  be  important  at 

20  -3 

densities  above  ^  10  cm  based  on  Fig.  3;  hence,  the  two  curves  in  Fig.  1 

21  23  -3 

for  ion  densities  of  10  and  10  cm  could  be  in  error  due  to  the  omission 

of  these  end  even  higher  levels.  This  would,  in  fact,  be  the  case  if  it  were 

not  for  lowering  of  the  ionization  potential  in  dense  plasmas.  Using  the  method 

(42) 

prescribed  by  Grlem  ,  the  ionization  lowering  was  determined  for  carbonlike 
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Al  VI 11  in  a  50  eV  plasma.  At  an  ion  density  of  10  cm  ,  the  lowering  Is 

21  -3 

only  about  5.0  eV;  however,  at  10  cm  ,  it  is  nearly  45.0  eV  and  the  reduced 

ionization  limit  falls  approximately  10  eV  above  several  of  the  n  -  4  excited 

states  in  this  ion,  pushing  higher  energy  states  into  the  continuum.  Shile 

these  estimates  are  only  approximate,  it  is  doubtful  that  gross  errors  are 

present  in  the  radiation  calculations.  Moreover,  it  Is  found  that  lithiumlike 

19  -3 

Al  XI  Is  the  strongest  L-shell  radiator  above  10  cm  ,  peaking  at  .about 

100  eV,  and  excited  levels  up  to  n  *  6  have  been  included  for  this  ion  to 

assure  approximately  correct  calculation  of  the  peak  intensity  at  high  den- 
21  3 

sltles  (at  10  ions/ cm  and  100  eV,  the  calculated  ionization  potential  of 
Al  XI  lies  Just  above  the  n  ■  5  level). 

As  plasma  density  increases,  the  level  populations  and,  thus,  the  radi¬ 
ation  emission  approach  values  that  would  be  predicted  by  an  LIE  model.  This 
trend  is  illustrated  in  Fig.  4,  where  we  have  used  an  LTE  calculation  with  the 
identical  level  structure  as  In  our  CRE  model  to  predict  the  optically  thin 

line  emission;  a  comparison  is  made  with  the  CRE  results  for  three  ion  densities. 
21  3 

At  10  ions/cm  Che  L-shell  line  emission  is  very  nearly  in  LTE  while  the 

X-shell  is  not.  At  10^  ions/ cm ^ ,  the  L-shell  is  in  LTE  and  the  K-sbell  shows 

(42) 

only  small  differences  from  LTE.  An  approximate  hydrogenic  formula  can 
be  used  to  determine  the  density  at  which  excited  states  go  into  LTE; 

H.  L9xl017  (|S)3  (^)1’  aT3  U0> 
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where  AE  is  Che  transition  energy  and  E^  is  the  ionization  potential  of 

hydrogen.  Using  the  n  •  3  state  of  A1  XI  at  120  eV  yields  an  electron  density 

2?  -3  .  21  -3 

ox  1.3  x  10  ~  cm  or  an  ion  density  or  -  2  x  IQ  an  which  gives  very 

reasonable  agreement  with  the  result  in  Fig.  4,  at  120  eV. 

In  the  interest  of  comparison  with  corona  and  LIE  models,  CRE  results 
have  been  presented  for  optically  thin  plasmas  only.  Photon  reabsorption, 
however,  will  play  a  more  important  role  for  L-shell  line  emission  than  for  the 
K-shell  in  dense  laboratory  plasmas.  In  order  to  study  the  effect  of  photo¬ 
excitation  and  stimulated  emission  on  the  radiation,  the  plasma  size  was 
varied  for  a  number  of  calculations,  since  the  photoexcitation  rates  are 
dependent  on  the  optical  depth  which  is,  in  turn,  dependent  on  the  plasma 

physical  extent.  Shown  in  Fig.  5  is  the  bound-bound  power  density  vs.  T 
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for  a  plasma  at  10  ions /cm  and  a  radius  which  varies  from  SO  pm  to  5000um 

(the  thin  curve  is  also  plotted).  Note  that  at  these  parameters,  the  K-shell 

is  only  slightly  affected  by  opacity  while  the  L-shell  emission  is  reduced  at 

peak  (120  eV)  by  factors  of  1/3,  1/6,  and  1/9  of  the  thin  result  as  the 

radius  increases  by  an  order  of  magnitude.  Since  the  optical  depth  is  also 

proportional  to  the  population  density  of  the  lower  state  of  the  transition, 

it  is  expected  that  these  opacity  effects  will  be  magnified  at  higher  densities, 

while  at  lower  densities,  the  emission  will  tend  toward  the  optically  thin 

values. 

In  addition  to  line  photons,  the  radiation  from  free-bound  transitions 
(due  to  radiative  recombination)  is  also  transported.  The  effects  of  opacity 
on  the  continuum  are  illustrated  in  Fig.  6  where  the  continuum  power  density 
is  shown  vs.  for  the  same  plasma  parameters  as  in  Fig.  3.  Note  that  in  the 
region  of  L-shell  emission  the  effect  is  reversed!  Since  reabsorption  and 
transport  of  line  photons  results  in  strong  optical  pumping  of  the  excited 
state  populations,  the  increase  In  ground  state  populations  via  collisional 
ionization  results  in  a  subsequent  increase  in  radiative  recombination  photons 
which  more  than  compensates  for  the  attenuation  of  continuum  photons  by  photo- 
ionization  processes  (photoionization  is  almost  negligible  at  these  plasma 
parameters).  Hence,  the  continuum  radiation  power  density,  being  strongly 
tied  to  the  excited  state  densities,  Increases  as  the  radius  of  the  plasma 
Increases.  This  anomalous  behavior  is,  however,  strongly  density  dependent. 
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and  photon  attenuation  by  photoionization  processes  will  overtake  the  pumping 
affects  described  above  as  the  plasma  density  rises.  Also  notice,  upon  com¬ 
parison,  that  while  K-sheli  line  emission  dominates  over  continuum  emission  -- 
by  about  a  factor  of  5,  the  continuum  radiation  from  L-shell  ions  is  comparable 
to  or  exceeds  the  line  photons  for  laboratory-size  plasmas  at  this  density. 

Having  examined  both  bound  and  continuum  contributions  to  the  emission 
individually,  we  are  now  in  a  position  to  understand  the  total  radiation 
emission  characteristics,  shown  in  Fig.  7,  for  identical  plasma  parameters 
as  in  Figs.  5  and  6.  The  most  striking  feature  is  the  relatively  small 

difference  between  total  radiated  power  density  for  the  50,  500,  and  5000  urn 
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plasmas  over  all  ranges  of  temperature  at  10  ioos/cm  .  Only  the  optically 
thin  approximation  differs  significantly,  and  only  by  a  factor  of  6  at  the 
L-shell  peak  at  130  eV.  A  consideration  of  the  behavior  of  line  and  continuum 
radiation  from  the  previous  two  figures  shows  quite  clearly  the  compensating 
effect  of  their  contributions. 

An  approximation  that  is  often  used  to  characterize  the  emission  from 

dense  plasmas  is  the  black-body  emission  rate,  given  by  eqn.  (9).  Converting 

the  flux  to  a  power  density,  the  black-body  (BB)  radiation  emission  is  also 

plotted  in  Fig.  7  for  the  three  relevant  plasma  radii.  While  the  slope  of  the 

BB  curves  is  similar  to  the  CHE  curves,  large  quantitative  errors  result  from 
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using  this  approximation  to  describe  a  plasma  at  10  ions/cm  except  at  tempera¬ 
tures  lower  than  15  eV  for  the  5000  pm  plasma.  It  is  conceivable,  however,  that 
one  could  use  the  BB  approximation  to  describe  the  plasma  emission  by  employing 
a  BB  temperature  lower  than  the  actual  plasma  electron  temperature.  For  example, 

a  black-body  radiation  function  would  approximately  describe  the  emission  of  a 
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5000  pm  plasma  at  10  cm  if  the  BB  temperature  used  was  about  half  of  the 

actual  temperature;  factors  of  about  1/3  and  1/6  could  be  also  used  for  the 

500  and  50  p a  plasmas,  respectively.  The  approximation,  of  course,  is  only 

19 

good  up  to  the  temperature  where  the  L-shell  radiation  peaks,  since,  at  10 

_3 

cm  ,  the  K-shell  is  substantially  removed  from  the  LTE  state  which  eqn.  (9) 
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characterizes.  The  approximation  is  still  quite  crude  at  10  ions/cm  because 
the  slope  of  the  BB  is  steeper  than  the  CHE  curve.  To  illustrate  the  trend 
toward  black-body  emission  as  density  increases,  a  curve  has  been  drawn  in 


Fig.  4  (dotted  line)  which  represents  the  T4  slope  of  the  3B  curve  adjacent  to 
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the  optically  thin  CSE  calculation  at  10  cm  .  Considering  that  the  compari¬ 
son  is  between  a  detailed  model  of  the  microprocesses  in  a  plasma  and  a  simple 
formula  describing  an  ideal  radiator,  the  agreement  is  significant. 

IV.  XUV  DIAGNOSTICS 

The  motivation  for  modeling  detailed  atomic  physics  and  radiation  pro¬ 
cesses  in  high  density  plasmas  is  twofold:  (1)  In  medium-to-high  Z  plasmas, 
not  only  is  radiation  a  significant  mechanism  for  the  loss  and  re— distribution 
of  plasma  energy,  but  a  large  proportion  of  this  energy  can  be  contained  in  ion 
internal  energy.  Consequently,  the  atomic  model  can  seriously  affect  the 
partitioning  of  energy  in  hydrodynamic  calculations  .  (2)  In  addition  to 

accurate  determination  of  this  partitioning,  a  detailed  atomic  model  has  the 
added  attractiveness  of  being  capable  of  accurate  prediction  of  the  details 
of  emission  which  can  then  be  used  to  diagnose  conditions  in  laboratory  plasmas. 

Merging  of  lines  and  continuum  ,  as  well  as  line  broadening and 
(45-47) 

shift v  ,  can  be  helpful  in  determining  plasma  densities.  Absolute  line 

intensities  can  be  used  to  predict  the  approximate  size  of  the  emitting  region 

(42) 

of  the  plasma.  Line  intensity  ratios  are  often  accurate  indicators  of  both 
temperature  and  density.  Line  ratios  have  been  used  extensively  in  diagnosing  the 
temperature  and  density  of  plasmas  which  emit  strongly  in  the  K-shell,  and  many 
studies^’  50)  ^ave  been  dong  to  quantify  soft  x-ray  line  ratios  for  this  use. 

Few  investigations  have  dealt,  however,  with  the  spectral  diagnostics  of  dense, 
medium  and  high  Z  plasmas  which  lie  in  the  XUV  region^^  .  In  this  section, 
we  investigate  several  possible  line  intensity  ratios  from  L-shell  aluminum  ions 
that  show  promise  as  density  and  temperature  indicators  for  plasmas  typical  of 
exploding  wire  and  laser-created  plasma  experiments. 

A.  Temperature  Determinations 

Appropriate  line  intensity  ratios  for  temperature  diagnosis  in  hot 
plasmas  are  relatively  abundant  in  most  XUV  spectra.  Since  the  ionic  popu¬ 
lations  displays  the  pronounced  "picket  fence"  effect  with  temperature  (as  does 
the  ionic  radiation,  as  seen  in  Fig.  2),  intensities  of  lines  from  adjacent 
ions  form  ratios  which  are  very  sensitive  to  changes  in  plasma  temperature. 
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This  vould  seem  to  suggest  that  any  number  of  pairs  of  lines  from  adjacent 
ions  are  effective  diagnostics;  in  fact,  there  are  further  restrictive  con¬ 
siderations  .  In  choosing  what  we  consider  to  be  good  diagnostics,  three 
additional  criteria  need  to  be  considered:  (1)  A  good  temperature  diagnostic¬ 
line  ratio  should  be  comprised  of  two  lines  which  lie  close  together  in  the 
spectrum.  This  eliminates  many  of  the  experimental  uncertainties  associated 
with  non-linear  response  of  crystal  or  film  due  to  differences  in  wavelength. 

(2)  The  two  lines  should  also  be  of  sufficient  intensity  to  be  easily  measured 
and  their  peak  intensities  should  be  comparable,  hence,  yielding  a  ratio  near 

1  in  the  temperature  range  of  the  plasma  where  the  lines  are  strong  emitters. 

(3)  Good  temperature  diagnostic  ratios  should  be  somewhat  independent  of  density 
variation  (over  a  limited  range  of  interest)  and  should  also  exhibit  little 
change  due  to  opacity  effects.  Of  the  three  conditions,  the  last  is  the  most 
difficult  to  satisfy. 

In  Table  III,  selected  transition  pairs  comprising  temperature  sensitive 
line  intensity  ratios  for  every  adjacent  pair  of  L-shell  ions  are  listed  which 
satisfy  the  first  two  criteria.  Since  individual  j-components  of  the  transitions 
are  not  treated  in  this  model,  the  contributions  over  a  number  of  wavelengths 
need  to  be  summed  in  several  cases  in  order  to  experimentally  obtain  the  total 
intensity  of  the  transition  as  calculated  by  the  model;  these  wavelengths  are 
also  given.  The  line  ratios  for  these  transitions  vs.  electron  temperature  for 
an  optically  thin  plasma  at  10  ions /cm  have  been  calculated  and  are  shown 
in  Fig.  8.  Note  that  these  ratios  provide  useful  diagnosis  over  a  wide  range 
of  temperatures  corresponding  to  the  range  over  which  L-shell  radiation  is 
prominent.  In  an  earlier  work^  on  aluminum  K-shell  diagnostics,  it  was  shown 
that  line  intensity  ratios  formed  from  lines  from  heliumlike  and  hydrogenlike 
ions  provide  insight  into  plasmas  temperatures  over  an  effective  range  beginning 
at  about  300  eV.  In  light  of  this  result  and  from  examination  of  Fig.  8,  there 
exists  a  large  temperature  range  from  about  150  eV  to  about  300  eV  in  which  no 
temperature  diagnostic  is  available.  This  is,  of  course,  the  range  which  lies 
between  the  lithiumlike  and  heliumlike  systems  -  the  barrier  between  the  two 
atomic  shells.  A  line  ratio  was  calculated  for  this  pair  of  ions,  but  A1  XII 
lines  don't  effectively  turn  on  until  temperatures  of  about  250-300  eV  are 
reached,  resulting  in  this  difficult  region  for  diagnosis.  However  this,  in 
itself,  is  informative  in  interpreting  spectra.  Namely,  any  time-resolved 


spectrum  in  which  both  lithiumlike  and  heliumlike  lines  occur  together 
(specifically  25-50  X,  in  aluminum)  «nd  show  comparable  intensities  is  an 
indication  of  a  rather  substantial  temperature  gradient  in  the  plasma  source. 

This  is  aoc  necessarily  true  for  lines  from  adjacent  ions  vi thin  the  L-shell. __ 
The  third  criteria  referred  to  earlier  involves  the  variation  of  the 
temperature-sensitive  ratios  with  density  and  opacity.  If  the  line  ratio  varies 
significantly  with  either  of  these  plasma  parameters,  the  problem  of  interpreting 
the  experimental  spectrum  becomes  a  very  complicated  pu2zle  in  which  the  assess¬ 
ment  of  size,  temperature,  and  density  are  all  strongly  coupled.  In  order  to 
illustrate  these  effects  on  the  line  ratios  listed  in  Table  III,  the 

2p-3s/2s^-2s3p(^P)  ratio  (#7)  was  selected  as  an  example.  This  ratio  was  calcu- 
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lated  from  50  to  200  eV  for  plasmas  of  10  ,10  ,  and  10  ions/cm  in  the 

optically  thin  approximation  and  at  10  ions/ cm  with  opacity  effects  for 

plasma  of  radii  50,  500,  and  5000ym;  the  results  are  displayed  in  Fig.  9. 

The  results  show  primarily,  that  the  line  intensity  ratio  depends  only  slightly 

on  the  plasma  density  between  10^  and  10^  ions/cm^.  The  ratio  tends  to 

diverge  more  as  the  temperature  increases,  so  that  at  170  eV,  near  the  upper 

temperature  limit  for  applicability  of  this  diagnostic  (the  A1  X  line  is  too 

faint  above  this  temperature),  the  ratio  at  10  cm  is  about  1.6  times  the 
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ratio  at  10  cm  .  If  one  were  to  use  the  low  density  curve  as  a  universal 

one,  and  attempt  to  diagnose  a  high  density  plasma,  errors  in  the  predicted 

temperature  of  25-30  eV  could  be  expected.  In  reality,  the  density  is  usually 

known  to  within  an  order  of  magnitude,  and,  using  the  results  of  this  model, 

one  could  anticipate  only  small  discrepancies  in  a  temperature  determination 

using  this  ratio.  Nevertheless,  the  density  effects  on  the  ratio  are  so  small 

within  the  range  of  Fig.  9,  that  practically  speaking,  the  ratio  is  density 

insensitive.  The  reason  for  this  insensitivity  is  not  immediately  apparent. 

Close  inspection  of  the  dominant  depopulating  rates  of  the  upper  levels  of  each 

transition  shows  that  over  the  4  orders  of  magnitude  change  in  density  in  Fig.  9, 

the  levels  traverse  from  the  coronal  region  to  the  collision-dominated 

2 

plasma  regime.  The  radiative  decay  coefficients  for  the  2p-3s  and  2s  -2s3p 

transitions  are  1.5  x  10^  sec  ^  and  4.6  x  10^  sec  \  respectively.  At  10 
5 

ions/ cm  and  100  eV,  the  respective  electron  collisonal  de-excitation  rates 
8  8  3 

are  5.3  x  10  and  2.8  x  10  cm  /sec,  well  below  the  radiative  decay  rates, 
implying  a  coronal  behavior  for  the  upper  level  populations.  Thus,  the  line 
intensities  are  given  approximately  by 
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where  u  and.  1  denote  upper  and  lower  stares  of  the  transition.  Bence,  at 

constant  temperature,  each  line  intensity  scales  as  the  total  ion  density 

21 

squared, and  the  line  ratio  remains  constant  at  coronal  densities.  At  10 

3  12 

ions/ca  ,  however,  the  collisional  de-excitation  rates  are  5.3  x  10  and 

12  3 

2.8  x  10  cm  /sec,  exceeding  the  radiative  decay  rates  by  factors  of  35  and 

6,  respectively.  At  this  density,  the  L  shell  ion  populations  are  approaching 

their  LTE  distribution.  In  LIE,  we  can  relate  the  2s3p  "h?  state  population 
2 

(N  )  to  the  2s  S  ground  state  population  (N  )  by  the  Saha  equation. 
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We  can  also  express  the  relationship  between  the  2s  and  3d  (N^)  levels  of 
A1  XI  via  Boltzmann's  Law  , 


N2  S2 

r  •  r  <-“02™ 

o  o 


Eliminating  N  and  g  ,  we  find, 
o  o 
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Thus,  in  strict  LTE,  the  ratio  of  the  populations  and,  hence,  the  line  inten¬ 
sities  (in  the  optically  thin  approximation)  scale  by  a  factor  of  ^ .  To 
illustrate,  we  have  plotted,  as  isolated  points,  the  line  ratio  at  densities 
of  10  and  10  ions/cm  and  a  temperature  of  100  aV  in  Fig.  9.  It  becomes 

apparent  that  this  line  ratio  loses  its  "density  insensitive"  quality  somewhere 
21  22  3 

between  10  and  10  ions/cm  as  the  levels  make  the  transition  from  the  collisional 
radiative  regime  to  that  of  LTE.  In  fact,  this  is  the  criteria  which  determines 
the  density  sensitivity  of  the  other  line  intensity  ratios  in  Table  III;  as 
two  adjacent  ions  approach  a  state  of  LTE,  the  line  ratio  formed  by  emission 
lines  from  each  of  them  becomes  very  dependent  on  the  plasma  density  and  unsuit¬ 
able  as  a  temperature  diagnostic  unless  the  density  is  known  to  good  accuracy. 
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The  effect  of  photoexcitation  on  the  line  ratio  is  also  shown  in  Fig.  9. 

As  the  temperature  increases,  it  is  generally  the  transition  from  the  lower-charge 
ion  which  tirst  becomes  attenuated  by  the  plasma  opacity,  causing  an  initial  in¬ 
crease  in  the  value  of  the  line  ratio  compared  to  the  optically  thin  value.  As 
the  temperature  continues  to  increase,  the  ionic  populations  change  and  the 
transition  from  the  higher-charge  ion  becomes  optically  thick,  reducing  the 
ratio  somewhat  as  compared  to  the  thin  value;  both  of  these  effects  are  seen  in 
Fig.  9.  Notice  that  the  larger  5000  ym  plasma  exhibits  a  ratio  closer  to  the 
thin  values  than  the  500  ym  plasma  from  about  50  to  150  eV.  This  is  a  function 
of  the  relative  values  of  oscillator  strengths  and  collisional  quenching  rates 
of  the  respective  transitions.  This  behavior,  therefore,  is  particular  to  these 
transitions  at  this  density,  and  is  not  characteristic  of  any  general  trend  in 
radiation  transport  effects.  The  absence  of  radiation  transport  can  lead  to  errors 
in  the  calculation,  as  seen  in  Fig.  9,  but  these  errors would  lead  to  temperatures 
which  are  typically  10  to  20  eV  higher  than  the  correct  calculation  yields.  At 
lower  densities,  the  tendency  toward  an  optically  thin  plasma  reduces  these  effects 
somewhat.  The  way  to  reduce  opacity  effects  in  line  intensity  ratios  is  Cl)  to 
choose  transitions  which  have  small  oscillator  strengths,  or  C2)  which  end  on 
levels  which  are  not  as  highly  populated  as  the  ground  state.  Unfortunately, 
this  results,  in  the  first  case,  in  choosing  transitions  which  may  be  too  weak 
to  detect.  The  second  case  is  not  useful  since  the  L-shell  excited  state  popu¬ 
lations  are  not  significantly  smaller  than  the  ground  states  in  the  density  regime 
of  interest  for  this  investigation.  Hence,  if  a  truly  accurate  determination  is 
required,  a  detailed  calculation  with  radiation  transport  is  necessary. 

3.  Density  Determinations 

The  criteria  discussed  in  the  previous  section  concerning  good  tempera¬ 
ture  diagnostic  line  ratios  apply  to  good  density  diagnostic  line  ratios  by 
simply  interchanging  the  words  "temperature"  and  "density".  Density  diagnostics, 
however,  were  much  more  difficult  to  locate  in  the  spectrum  since  the  variation 
of  line  intensities  with  density  is  not  as  straightforward  as  it  is  with  temper¬ 
ature.  In  addition,  several  of  the  diagnostic  ratios  first  selected,  although 
density-sensitive,  displayed  only  small  changes  in  the  value  of  the  line  ratio 
with  density  variations;  thus,  experimental  error  would  introduce  large 


j 
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uncertainties  in  density  determinations.  We  found  we  were  most  successful  by 

choosing  lines  which  originated  from  different  spin  systems  within  the  same 

ion.  Because  of  the  forbidden  nature  of  the  transitions  of  different  spin-system 

levels  to  the  ground  state,  density  effects  were  often  profound  in  determining 

the  population  of  the  excited  states.  Lines  from  the  same  ion  were  chosen  to 

ensure  that  the  radiation  from  both  lines  was  emitted  from  the  same  location  in 

the  plasma.  A  problem  which  occurs  in  these  diagnostics  for  L-shell  ions  is  the 

narrow  range  of  temperatures  over  which  the  ions  exisc  (see  Fig.  2).  As  the 

plasma  temperature  deviates  slightly  from  the  small  interval  over  which  the 

transitions  radiate  strongly,  the  line  ratios  become  extremely  temperature 

dependent.  Three  pairs  of  lines  were  selected  as  density  diagnostics  over  a 

limited  range  of  temperatures,  however,  and  they  are  listed  in  Table  IV  along 

with  the  relevant  wavelengths.  The  line  intensity  ratios  for  these  transitions 

are  plotted  vs.  ion  density  in  Figs.  10-12  for  an  optically  thin  plasma  for 

several  electron  temperatures  at  which  the  lines  are  strong  radiaters.  In 

Fig.  10,  the  two  transitions  are  from  nitrogenlike  A1  VII;  from  Fig.  2,  the 

line  emission  from  this  ion  peaks  at  about  35  eV  and  the  line  ratio  at  30  and 

18  3 

40  eV  in  Fig.  10  are  reasonably  close  above  3  x  10  ions /cm  .  The  ratio  at 
50  eV  displays  larger  deviations  from  the  other  curves  above  this  density,  since 
the  A1  VII  emission  is  already  15  eV  removed  from  the  peak  emitting  temperature. 
Note  also  the  double-valuedness  of  the  ratio  complicates  spectral  interpretation 
<ven  further.  This  ratio  is  typical  of  the  transition  pairs  we  Investigated 
from  A1  IV  -  A1  VII  ions;  significant  temperature  dependence  over  even  small 
ranges  of  T@.  In  Figs.  11  and  12,  the  line  ratios  show  much  less  temperature 
sensitivity  and  no  double-valuedness  in  the  useful  density  region.  In  Fig.  12, 
the  ratio  changes  very  little  over  a  temperature  range  of  80  eV,  causing ,  for 
example,  an  error  of  at  most  a  factor  of  5  in  the  density  determination  if  the 
experimental  line  ratio  was  2.0.  It  is  assumed,  however,  that  the  temperature 
could  be  determined  to  better  accuracy  than  that  from  diagnostics  presented 
earlier,  and  the  density  could  probably,  be  determined  to  less  than  a  factor  of 
2  error  by  a  comprehensive  analysis  of  the  spectrum.  Agreement  between  many 
aspects  of  the  theoretical  and  experimental  spectrum  is  in  fact  the  only  way  to 
insure  that  the  plasma  has  been  accurately  characterized  by  the  model.  In 
the  following  section,  this  point  is  explored  in  relation  to  the  generation  of 
theoretical  spectra  with  CRE  calculations. 
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CALCULATION  OF  THEORETICAL  SPECTRA 


In  the  previous  sections,  we  presented  the  results  of  L-shell  aluminum 
calculations  which  we  daeoed  useful  in  gaining  understanding  of  cne  radia¬ 
tion  energetics  in  dense  plasmas  over  a  wide  range  of  temperatures,  as  well 
as  providing  diagnostic  information  applicable  to  the  determination  of  plasma 
parameters.  Although  one  can  often  make  a  strong  case  for  an  attempt  to 
characterize  a  plasma  via  an  Isolated  diagnostic,  we  are  of  the  opinion  that 
ionization-radiation  dynamic  calculations  have,  at  present,  reached  a  level 
of  sophistication  that  allows  for  direct  comparisons  of  theoretically  generated 
spectra  and  those  obtained  experimentally.  Moreover,  it  is  our  contention 
that  only  by  this  comprehensive  replication  of  line  and  continuum  radiation 
signatures,  that  anything  definitive  can  be  said  about  the  ability  of  the 
theoretical  calculations  to  characterize  the  experimental  plasma.  In  this 
section,  theoretical  spectra  calculated  by  our  radiation  model  for  a  plasma 
with  a  temperature  gradient  similar  to  that  found  in  Imploding  wire  array 
plasmas  are  presented  and  discussed.  These  spectra  are  not  intended  to  be  an 
attempt  to  model  any  individual  wire  implosion  but  are  meant  to  be  illustra¬ 
tive  of  the  ability  of  the  theoretical  model  to  produce  spectra  suitable  for 
comparison  with  experiment  using  only  plasma  temperature,  density,  and  size 
of  emitting  region  as  input  parameters.  For  examples  of  actual  experimental 

spectra  from  imploding  wire  plasmas  and  other  dense  plasma  sources,  the  reader 

13 

is  referred  to  an  upcoming  paper  by  D.  Nagel. 

Within  the  probability  of  escape  phenomenology  employed  by  this  model,  a 
radiation  spectrum  can  be  calculated  by  summing  the  frequency-by- frequency 
contributions  from  each  plasma  zone  for  line  emission,  recombination  radiation, 
and  bremsstrahlung.  Only  the  bremsstrahlung  is  not  transported  through  the 
plasma  since  the  free-free  opacity  is  small  for  the  density  and  temperature 
regime  of  interest  here.  The  contribution  from  each  zone  is  determined  from 
the  radiative  source  function  and  the  probability  of  escape  to  the  vacuum  region 
beyond  the  outer  boundary  of  the  plasma.  This  zone-by-zone  treatment  produces 
a  spectrum  in  which  lines  are  superpositions  of  a  number  of  different  profiles 
and  the  continuum  is  a  superposition  of  a  number  of  different  slopes,  depending 
on  the  zoually  determined  temperature  and  density  as  well  as  the  transport 
dependent  atomic  state  population  densities.  The  results  are  displayed  as 
spectral  intensity  (in  watts/cm  /angstrom)  radiated  into  4ir  radians  and  plotted 
versus  wavelength. 
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The  temperature  gradient  for  this  example  vaa  chosen  so  as  to  approximately 
describe  an  Imploded  cylindrical  aluminum  plasma  which  as  assembled  on  axis 
and  is  at  the  point  of  maximum  emission  in  its  cine  evolution.  The  plasma  is 
taken  to  be  1  mm  in  diameter  with  a  linear  temperature  gradient  from  80  eV 
in  the  central  region  to  30  eV  near  the  periphery.  A  hot  K-line  emitting  core 
is  usually  present  in  these  plasmas  but  was  neglected  here  since  It  is  our 
purpose  to  illustrate  the  emission  solely  from  the  L-shell  lines  in  this  study; 
however,  the  K-lines  are  easily  included  in  the  calculation  by  altering  the 
temperature  profile  to  reflect  K-shell  emission  temperatures.  A  comparison 
of  theoretical  and  experimental  spectra  from  aluminum  K-shell  ions  is  the  sub¬ 
ject  of  another  study^.  Although  a  density  gradient  most  likely  exists  in 
actual  imploded  plasmas,  none  was  Included  in  this  example  in  order  to  eliminate 
confusion  over  which  effects  are  due  to  a  temperature  gradient  and  which  to  a 
density  profile.  Moreover,  a  temperature  profile  alone  is  sufficient  to  illus¬ 
trate  the  varied  line  emission  from  the  L-shell  Ions  as  would  be  recorded  by 
a  spectrometer  for  an  actual  experiment. 

The  calculated  theoretical  spectra  for  the  above  conditions  are  reproduced 
in  Figs.  13-13.  Approximately  30  mi  of  source  and  instrumental  broadening  has 
been  convoluted  with  the  line  profiles  to  duplicate  similar  conditions  as  would 
be  present  in  experimental  spectra,  since  this  broadening  mechanism  often 
dominates  line  profiles  for  plasmas  of  these  temperatures  and  densities.  In 
Figs.  13  &  14,  the  spectral  intensity  has  been  plotted  on  a  linear  y-axis  as 
would  an  experimental  spectrum;  the  wavelength  range  has  been  separated  into  the 
30-65  5  and  63-105  5  regions  in  order  to  show  better  detail.  The  30  X  cutoff 
just  includes  the  A1  XI  2s-6p  line,  the  shortest  wavelength  line  in  our  L-shell 
model  (excluding  the  satellite  line  from  the  doubly-excited  ls2p2  h)  state). 

If  K-shell  emission  intensity  had  been  comparable  to  the  L-shell  lines,  however, 
several  Balmer  series  lines  from  the  hydrogenlike  and  heliumlike  ions  could  be 
seen  among  the  L-shell  lines  ranging  from  about  25  to  52  X.  Notice  that  we  have 
allowed  several  very  prominent  lines  to  "saturate"  our  theoretical  "spectrometer 
crystal"  in  order  to  display  some  of  the  lower  intensity  lines  of  the  spectra. 

The  jagged  base  line  is  due  to  recombination  edges  forming  the  free-bound 
continuum  (the  bremsstrahlung  is  very  much  lower);  in  actual  experimental  spectra, 
the  noise  level  often  washes  out  much  of  the  finer  detail  of  this  feature,  while 
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the  Absorption  edges  of  filters  frequently  distort  the  true  shape  of  this 
emission  profile.  Identif icacion  of  the  individual  lines  is  accomplished 
by  numbering  in  order  of  increasing  wavelength,  with  Table  V  serving  as 
a  key;  also  given  in  Table  V  is  the  transition  wavelength  and  the  optical 
depth  at  line  center,  calculated  from  the  central  axis  of  the  plasma  to  the 
outer  boundary.  The  wavelengths  quoted  in  this  table  are,  in  some  cases, 
averages  over  a  number  of  transitions  with  similar  initial  and  final  electron 
configurations  but  different  angular  momentum  states;  hence,  they  may  not 
correspond  exactly  to  values  quoted  by  standard  references  for  a  given 
transition. 

From  Fig.  13,  it  is  seen  that  the  strongest  emitters  in  this  plasma  are 
the  lithiumlike  2p-3d  transition  (#10),  the  berylliumlike  232p-2s3d,  -3s 
(singlet)  transitions  (#16),  and  the  boronlike  2p-3d,  -3s  (doublet)  transitions 
(#17).  For  ths  latter  two  lines,  ths  intensity  is  due  almost  entirely  to  the 
transitions  from  tha  3d  levels,  since  both  the  population  of  these  states  and 
radiative  decay  rates  from  them  exceed  those  of  the  corresponding  3s  levels. 

An  interesting  feature  that  became  apparent  when  comparisons  were  made  with 
experimental  spectra  was  the  relative  heights  of  the  A1  XI  2s-5p  (#2)  and  2p-6d 
-6s  (#3)  lines.  In  all  spectra  examined,  the  shorter  wavelength  line  always 
exceeded  the  longer  one  in  intensity.  Tet,  when  optically  thin  calculations 
were  anlayzed,  the  longer  wavelength  line  dominated  at  all  temperatures  and 
over  a  wide  range  of  densities.  It  was  only  when  opacity  effects  were  included 
in  the  analysis  that  the  emission  lines  reversed  their  relative  intensities  and 
agreed  with  experimental  data.  The  phenomenon  is  explained  by  the  larger 
optical  depth  of  the  2p-6d  line;  since  both  the  2s-5p  and  2p-6d  oscillator 
strengths  are  comparable  (about  a  2:3  ratio),  the  difference  in  opacity  lies 
with  the  population  density  of  the  lower  state  of  the  transition.  At  these 
densities,  the  2s  and  2p  levels  are  populated  approximately  according  to  their 
statistical  weights  at  temperatures  where  A1  XI  lines  emit  strongly.  Con¬ 
sequently,  the  2p-6d  optical  depth  will  be  about  2  to  3  times  larger  than  that 
of  the  2s-3p  transition  and  the  former  line  will  be  attenuated  more  as  the 
plasma  becomes  increasingly  optically  dense.  Initially,  it  was  thought  that 
this  feature  could  be  useful  as  an  opacity  diagnostic  to  estimate  the  size  of 


Che  emitting  region;  however,  Che  relative  changes  in  intensities  of  the  two 
Hoes  vi.cn  plasma  sice  was  found  Co  he  too  small  to  derive  sufficiently 
accurate  quantitative  information  to  be  able  to  make  a  definitive  statement  •- 
about  this  plasma  parameter. 

In  Fig.  14,  lines  from  the  lower  ionization  stages  are  prevalent,  with 
A1  VI  -  Al  VIII  dominating  the  emission.  The  lines  from  A1  III  -  A1  V  do 
not  appear  with  sufficient  intensity  to  be  seen  since  these  ions  are  dominant 
at  temperatures  lower  than  the  30  eV  cutoff  in  our  profile  used  in  this 
example.  In  fact,  this  lower  limit  on  the  temperature  gradient  was  chosen 
after  several  examples  of  data  from  aluminum  wire  experiments  were  scanned 
and  few  Intense  contributions  to  the  line  spectrum  from  these  ions  found. 

In  order  to  illustrate  the  gross  characteristics  of  the  radiation  spectrum 
from  this  sample  run.  Fig.  15  was  plotted  with  a  lograithmic  y-axis ,  exhibiting 
low  intensity  lines  as  well  as  the  stronger  ones.  At  the  longer  wavelengths, 
above  about  240  lines  due  to  An«0  transitions  among  the  n«2  levels  of  the 
L-shell  ions  are  prominent  along  with  a  few  weak  lines  from  higher  Rydberg 
state  transitions  in  Al  XI.  The  real  value  of  a  comprehensive  spectrum  as 
shown  in  Fig.  15  is  in  its  ability  to  indicate,  at  a  glance,  which  wavelength 
regions  are  responsible  for  the  bulk  of  the  emission.  In  particular,  the  wave¬ 
length  dependence  of  the  free-bound  continuum  radiation,  which  was  seen  from 
Fig.  6  to  constitute  a  significant  fraction  of  the  total  emission,  is  very 
clearly  illustrated.  This  extensive  picture  of  the  radiation  emission  as  a 
function  of  wavelength  is  extremely  useful  when  knowledge  of  the  emission  from 
particular  "energy  bins"  is  required.  In  particular,  single-temperature  analytic 
models  of  the  plasma  dynamics  often  characterize  the  plasma  radiation  in  a  very 
approximate  way,  such  as  assuming  that  all  the  line  emission  is  radiated  at  a 
single  frequency  and  an  appropriate  optical  depth,  in  attempt  to  mock  up,  in 
a  concise  manner,  the  actual  emission  energetics  while  taking  some  opacity 
effect  into  account.  Based  on  runs  made  with  a  flat  temperature  profile,  we 
find  that  comprehensive  spectral  plots  display  even  more  localized  emission  than 
is  shown  in  Fig.  15  for  the  exemplary  temperature  profile.  Thus,  these  plots 
serve  to  quantify  the  overall  plasma  emission  in  a  way  that  may  easily  be 
implemented  in  a  simple  analytic  calculation. 
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Although  qualitative  comparisons  have  been  made  using  the  spectra 
generated  by  this  model  and  a  few  random  experimental  spectra,  no  attempt 
ha3,  as  ysc,  ba^n  >T_ada  to  quantitatively  compare  the  theory  with  a  specific 
experiment  In  order  to  diagnose  the  plasma  conditions  at  the  time  of 
emission.  Recently,  however,  spectroscopists  at  many  laboratories  have 
realized  the  wealth  of  information  contained  in  the  XUV  region  of  the  emission 
spectrum  from  dense  plasmas,  and  have  begun  to  make  serious  attempts  to  obtain 
good  spectral  data  in  that  wavelength  range  as  well  as  in  the  soft  X-ray 
region  dominated  by  the  K-shell  lines.  It  is  our  hope  that  through  direct 
comparison  of  theoretical  and  experimental  data,  this  model  may  make  adequate 
use  of  that  information  and  shed  light  on  the  continuing  problem  of  the 
characterization  of  dense  plasma  dynamics  and  the  determination  of  basic 
plasma  parameters . 
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Figure  I. 

Figure  2. 

Figure  3. 

Figure  4. 

Figure  5. 
Figure  6. 

Figure  7. 

Figure  8. 
Figure  9. 

Figure  10. 

Figure  11. 

Figure  12. 


T3c.1L  line  amission  power  coefficient  vs.  electron  temperature 
for  densities  from  iO3*-  iQ2-3  ions/ cm3  (optically  chin 
approximation) . 

Ionic  line  emission  power  densities  vs.  electron  temperature  for 
a  plasma  at  10 19  ions/ cm3  (optically  thin  approximation);  total 
L-shell  emission  is  also  shown  (broken  line) . 

Individual  line  emission  power  densities  for  A1  VIII  (caxbonllke) 
vs.  electron  temperature  for  a  plasma  at  10 19  ions/cm3  (optically 
thin  approximation) ■  See  Table  II  for  line  identifications. 

Comparison  of  CRE  (solid)  and  LIE  (broken)  line  emission  power 
coefficients  for  plasma  densities  from  1019  -  10 23  ions  cm3 
(optically  chin  approximation) . 

Total  line  emission  power  densities  (with  opacity  effects)  for 
a  plasma  cylinder  at  10 19  ions/cm3  and  radii  from  50  -  5000  ym. 

Total  free-bound  continuum  emission  power  densities  (with  opacity 
effects)  for  a  plasma  cylinder  at  IQ39  ions/cm3  and  radii  from 
50  -  5000  um. 

Total  line  plus  continuum  emission  power  density  (with  opacity 
effects)  for  a  plasma  cylinder  at  1019  ions/cm3  and  radii  from 
50  -  5000  ym.  Power  density  from  similar  plasmas  in  the  black- 
body  limit  are  also  shown  (straight  curves  in  upper  left). 

Temperature-sensitive  line  intensity  ratios  for  selected  L-shell 
lines  at  10 19  ions/cm3  (optically  thin  approximation). 

Line  intensity  ratio  for  the  2p-3s/2s2-2s3p(1P)  transitions  for 
densities  of  1017,  10 19  and  10 21  lons/cm3  (optically  thin  approxi¬ 
mation)  and  for  a  plasma  cylinder  at  1019  ions/cm3  (with  opacity 
effects)  for  radii  of  50  -  5000  ym. 

Line  intensity  ratio  for  the  2s2p(2P)  -  2s3d(2D)/2s2p(2P)  -  2s3d(2D) 
transitions  vs.  ion  density,  for  several  electron  temperatures 
(optically  thin  approximation) . 

Line  intensity  ratio  for  the  2s22p2(2?)  -  2s^2p3d(2D)/2s2p2(^S)  - 
2s2p23d(5P)  transitions  vs.  ion  density,  for  several  electron 
temperatures  (optically  chin  approximation). 

Line  intensity  ratio  for  the  2s22p2(^S)  -  2s22p23s(^P)/ 

2s22p3(2D)  -  2s  22p  2(  1D)  3s  (  *0)  transitions  vs.  ion  density,  for 
several  electron  temperatures  (optically  thin  approximation) . 
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Figure  Captions  (continued) 


Figure  13.  Aluminum  spectrum  from  30-65  angs crons  for  a  plasma  ac 
10 19  ions/cm3  and  a  linear  temperature  gradient  from 
30  to  80  eV. 


Figure  14. 


Aluminum  spectrum 
10 19  ions /cm3  and 
30  to  80  eV. 


from  65-105  angstroms  for  a  plasma  at 
a  linear  temperature  gradient  from 


Figure  15.  Aluminum  spectrum  from  0-700  angstroms  for  a  plasma  at 
10 19  ions/ cm3  and  a  linear  temperature  gradient  from 
30  to  80  eV. 
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appendix  d 


Ms.  Miriam  Gersten 
Maxwell  Laboratories 
9244  Balboa  Avenue 
San  Diego,  California  92123 

Dear  Miriam: 

I’ve  enclosed  a  brief  calculation  that  compares  the  optimum  radiation 
emission  one  could  expect  from  a  boron  plasma  radiating  in  the  K-shell 
and  an  aluminum  plasma  radiating  in  Che  L-shell.  The  information  sources 
are: 

Boron  -  "Systematic  Trends  in  the  Radiative  Emission  Rates 
for  Low-  and  Medium-2  Elements  in  High-Temperature 
Plasmas",  J.  Davis  and  V.  L.  Jacobs  (to  be  published 
in  JQSRT) . 

Aluminum  -  "Soft  X-ray  and  XDV  Production  from  High  Density, 
High-Temperature  Aluminum  Plasmas",  D.  Duston  and 
J.  Davis  (submitted  for  publication). 

The  boron  K-shell  line  emission  peaks  at  80  eV  and  the  power  coefficient 
there  is  1.3  x  10“28  watt  •  cm3. 


The  aluminum  L-shell  line  emission  peaks  at  120  eV  and  the  power 
coefficient  is  2  x  10“ 27  watt  •  cm3. 

19  3 

For  equal  radiation  output  as  an  aluminum  plasma  at  10  ions/ cm  the 

following  boron  density  is  required: 

•  28 

(1.3  x  10  )  N  N.  (boron)  ■  (2  x  10  )  N  N .  (aluminum) 

ex  ex 

for  boron  N  %  4.5  S.  in  the  K-shell;  for  aluminum  N  ^  10  N  in  the  L-shell 
(lithiunlike)  1 


Hence , 


(boron) 


(2x10  ^)  (10)  (aluminum) 


(1.3x10* 28)  (4.5) 


or. 


N^(boron)  *  5.84  N^( aluminum) 


for  equal  line  emission. 


There  are  3  contingencies,  however: 

1.  Both  calculations  are  done  "optically  thin";  we  expect  the  attenuation 
due  to  opacity  to  affect  both  plasmas  equally  though. 

2.  Continuum  photons  are  not  included;  our  best  guess  indicates  that 
there  will  be  more  A1  L-shell  continuum  than  B  K-shell,  though. 

3.  If  your  A1  plasma  was  cooler  them  120  eV,  the  L-shell  would  not  be 
radiating  optimally,  and  it's  possible  your  boron  plasma  emission 
could  be  comparable. 

The  conclusion  is,  the  lower-Z  stuff  shouldn't  produce  XUV  from  the 
K-shell  as  efficiently  as  medium-Z  does  from  the  L-shell. 

Sincerely, 

A 

/. 

c 

Dwight  Duston 
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DIRECT  SOLUTION  OF  THE  EQUATION  OF  TRANSFER 
USING  FREQUENCY-  AND  ANGLE-AVERAGED  PHOTON 
ESCAPE  PROBABILITIES.  WITH  APPLICATION  TO  A 
MULTISTAGE.  MULTILEVEL  ALUMINUM  PLASMA 
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I Receiled  17  August  1079) 

Abstract— A  formalism  is  developed  which  permits  direct  steady-stale  solution  of  the  transfer  equation 
using  escape  probabilities  averaged  over  angle  and  frequency.  A  matrix  of  probability-based  coupling 
coefficients,  which  are  related  to  the  kernel  function  K,.  is  used  to  obtain  the  source  function  for  a  Doppler 
profile  in  plane  parallel  geometry.  Comparison  is  made  with  exact  solutions,  establishing  the  high  accuracy 
of  the  technique.  The  method  is  extendable  to  different  physical  situations  by  simply  modifying  the 
coupling  coefficients.  As  an  example  of  a  more  realistic  application  of  the  formalism,  we  have  solved  for  the 
ionization-excitation  state  of  a  planar  aluminum  plasma  at  600  eV  in  collisional-radiative  equilibrium.  The 
results  agree  well  with  those  obtained  from  the  conventional  multifrequency-multiangle  formalism.  Ad¬ 
ditionally.  we  have  used  the  technique  to  gauge  the  effects  of  transport  of  lines  connecting  excited  states 
with  each  other. 


I.  INTRODUCTION 

In  the  literature  of  radiative  transfer  theory,  much  attention  has  been  devoted  to  photon 
escape-probability  concepts  as  a  means  of  aiding  in  the  interpretation  of  detailed  solutions  to 
the  transfer  equation  in  astrophysical  and  laboratory  contexts.  In  this  paper,  we  develop  a 
formalism  which  permits  utilization  of  such  concepts  to  obtain  direct  solutions  of  the  steady- 
state  source  function.  When  applied  to  a  plane-parallel.  Doppler-broadened  medium,  the 
technique  eliminates  entirely  the  need  for  either  frequency  or  angle  integration,  while  yielding 
generally  excellent  agreement  with  previous  results  obtained  by  Avrett  and  Hummer.1 

2.  MULTICELL  RADIATIVE  COUPLING  EQUATIONS 
If  a  medium  which  is  finite  in  extent  in  at  least  one  direction  is  divided  into  a  number  .V  of 
smaller  regions,  line  photons  emitted  in  each  region  will  have  a  finite  probability  of  being 
absorbed  in  any  of  the  other  N  -  1  regions,  as  well  as  being  reabsorbed  in  the  local  region  or 
escaping  entirely  from  the  medium.  Let  C„  be  the  probability  that  a  line  photon  emitted  in 
geometrical  region  i  traverses  the  distance  between  regions  i  and  /  and  is  absorbed  in  region  j. 
.V„,  will  refer  to  the  total  upper  level  population  of  region  i.  and  ,-t«,  stands  for  the  spontanous 
transition  probability  (sec1)  for  the  line  in  question.  In  each  region  i.  W„  and  D,  will  stand  for 
the  total  collisional  population  and  depopulation  rates,  respectively,  of  the  upper  level  (.VUI). 
The  above  concepts  and  notation  clearly  imply  the  following  equation  for  the  rate  of  change  of 
the  line  upper  level  population  in  region  i: 

^  =  .V„  W,  *  £  ,VU,.-L,C„  -  Xt.-L,  -  D,).  1 1 ) 

The  first  two  terms  on  the  r.h.s.  account  for  local  collisional  population  of  the  upper  level 
and  radiative  population  of  the  level  by  photons  emitted  from  ali  the  cells  j  into  cell  i. 
respectively.  The  final  term  ts  the  sum  of  spontaneous  and  collisional  depopulation  of  the  upper 
level.  In  this  paper,  we  will  confine  ourselves  to  the  case  where  d.VJdf  =  0.  that  is.  utilization 
of  Eq.  ( I)  to  obtain  the  steady  state  upper  level  population  in  each  region,  which  is  equivalent  to 
finding  the  source  function  when  the  lower  level  population  is  known.  For  a  large  class  of 
problems,  such  as  those  treated  tn  Avrett  and  Hummer,  the  ground  state  (lower  level) 
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population  is  much  larger  than  that  of  the  upper  level  and  the  absorption  coefficient  is  known 
with  high  accuracy.  In  these  cases,  the  steady-state  version  of  Eq.  (1)  requires  for  its  solution 
the  inversion  of  one  J V  x  .V  matrix— to  solve  S  linear  equations  in  S  unknowns — once  the  C„ 
are  calculated  from  the  absorption  coefficients.  When  the  absorption  coefficient  is  not  known  a 
priori,  the  appropriate  version  of  Eq.  ll)  must  be  set  up  for  each  level  considered  and 
populations  for  each  level  must  be  obtained  iteratively  using  the  C„  calculated  from  the 
populations  (hence,  absorption  coefficients)  of  the  previous  iteration. 

Aside  from  the  atomic  constants,  calculation  of  the  C„  from  the  absorption  coefficients  to 
permit  solution  of  Eq.  (1)  forms  the  heart  of  the  mathematical  problem  and  is  now  detailed. 

J  CALCULATION  OF  THE  COUPLING  CONSTANTS 
For  the  case  of  plane-parallel  geometry  and  Doppler  profile  with  a  spatially  constant 
Doppler  width,  it  has  proven  possible  to  develop  very  fast  and  extremely  accurate  algorithms 
for  obtaining  the  C„,  primarily  because  of  the  vast  amount  of  effort  by  other  workers  which  has 
been  devoted  to  probability  concepts  for  this  case. 

Consider  a  line  photon  emitted  at  a  point  in  cell  j,  in  the  direction  of  cell  /,  in  the  medium. 
The  probability  P„  that  the  photon  is  absorbed  in  ceil  i  is  given  by 

P„=P,(T„)-P,(7„-rA7,).  (2) 

where  P,(r)  is  the  angle-averaged  probability  that  a  photon  tranverses  an  optical  depth  r  without 
being  absorbed  or  scattered.  In  Eq.  (2).  r,,  is  the  line-center  optical  depth  from  cell  j  to  the 
boundary  of  cell  i  closest  to  cell  j  and  St,  is  the  optical  depth  of  cell  f  itself.  The  practicality  of 
the  method  here  described  is  obviously  dependent  upon  having  an  efficient  technique  for 
obtaining  the  P/s.  Holstein2  obtained  an  expression  for  the  monodirectional  escape  probability 
P,  valid  for  completely  redistributed  Doppler-profile  line  photons  at  large  optical  depths.  We 
have  numerically  calculated  the  integral 


*«•£'•(?)*•  * 

for  a  large  range  of  optical  depths  and  have  adopted  the  following  algorithm  to  obtain  P, 
quickly  and  accurately. 

(a)  At  70s3.  a  cubic  spline  polvnomial  has  been  fitted  to  the  exact  result  at  20  points 
roughly  equally  spaced  from  =  0  to  7„  =  3. 

(b)  For  r0  >  3,  the  following  analy  tic  expression  is  accurate  to  better  than  for  r,<  5  and 
4%  for  5  5  7tf  <  3  x  10*-. 


P,l7„) 


0.286 

7, A  (In  ( 1.95  7„i) 


At  this  point,  it  should  be  noted  that  in  calculating  the  C„  coupling  cells  of  finite  width.  P,  must 
be  averaged  over  the  cell  originating  the  photons.  This  is  especially  important  when  the 
originating  and  receiving  cells  are  adjacent  since  P,  can  varv  by  an  order  of  magnitude  or  more, 
in  some  cases,  from  the  front  to  the  back  of  the  originating  cell.  This  averaging  process  is  easily 
accomplished  by  analytically  integrating  the  above  expressions  across  the  originating  cell  j  so 
that 


(4) 


In  Eq.  (4),  Qr)y  and  (Ar),  are  the  line-center  perpendicular  optical  depths  of  cells  j  and  i  and  r„ 
is  the  optical  depth  between  ceils  /  and  i.  measured  between  the  two  closest  boundaries.  The 
factor  of  1/2  in  Eq.  (4)  accounts  for  the  assumed  equal  probability  of  photon  emission  in  either 
direction  from  cell  j. 
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4.  RELATION  TO  THE  EXACT  TRANSFER  EQUATION 
Avrett  and  Hummer1  have  written  the  formal  solution  for  the  source  function  Sir)  in  a 
plane-parallel  atmosphere  of  optical  depth  T  as 


Sir) 


=  (l-Po)[T 

Jo 


K\\(t  -  r)|S(f)  dr  t  P,)B. 


(5) 


where  B  is  the  Planck  function  for  the  local  electron  temperature.  PQ  is  the  "quenching 
parameter",  or  probability  per  scattering  that  the  photon  is  lost  from  the  line  and  K \  is  the 
kernel  function.  In  Ref.  1,  Avrett  and  Hummer  developed  a  useful  asymptotic  expansion  for 
the  Doppler  kernel  function  and  showed  that,  for  large  r  (measured  at  line  center), 


K,  (r) 


I 

4r’  Vlrr  In  r)' 


(6) 


Inspection  of  Eq.  (5)  reveals  that  Ki(r)  is  the  analog  of  the  discrete  probability-based  coupling 
coefficients  which  are  used  in  the  present  treatment.  K\(t)  couples  the  regions  of  the  medium 
together  and,  since  the  integral  is  carried  out  over  r.  it  is,  in  analogy  with  Eqs.  (1)  and  (2). 

1  dP 

K.(r0)--±2p  .  ,7) 

-  dr 

The  rate  of  change  of  the  angle-averaged  escape  probability  P,  across  optical  path  r0 
determines  the  efficiency  with  which  photons  are  absorbed  per  unit  optical  depth  after  crossing 
the  path  r0. 

The  analytic  expression 


P<(To)  =  0.286/roV[ln(l.95r,))].  18) 

which  we  have  used  computationally  in  the  calculations  presented  below  should,  when 
differentiated,  approach  Eq.  (7)  closely  in  the  limit  of  large  r. 

Differentiating  Eq.  (8)  yields,  as  r  —  *. 

AlL= _ 4 _  ,9) 

2drr„  6.993r- v'lln  r)’ 

(r  measured  at  line  center) 


which  differs  by  only  1.4%  from  the  exact  expression  (6).  For  most  cases  of  interest,  nearly  all 
of  the  coupling  occurs  within  optical  depths  much  smaller  than  3  x  10J,  where  Eq.  (8)  is  usually 
accurate  to  better  than  4%.  Also,  as  seen  below,  errors  of  a  few  per  cent  in  the  coupling  matrix 
do  not  result  in  large  source  function  errors  when  the  method  is  applied  to  specific  media  of 
finite  optical  depth. 


5  COMPARISON  WITH  THE  TWO-LEV  EL- 4  TOM  MODEL 
To  explore  the  computational  viability  of  the  approach  detailed  above,  we  have  applied  Eq. 
(1)  to  a  range  of  plane-parallel  media  of  varying  optical  depths  and  quenching  parameters.  The 
cases  presented  in  Figs.  1  and  2  reflect  comparisons  of  our  approach  with  exact  solutions 
obtained  by  Avrett  and  Hummer,1  The  physical  interpretation  of  the  solutions  has  been 
thoroughly  discussed  in  Ref.  I  by  Avrett  and  Hummer.  Since  the  media  are  symmetric  about 
the  midplane,  the  calculation  has  been  set  up  by  establishing  25  or  75  cells  in  half  the  medium 
and  coupling  each  cell  "to  itself"  and  to  others  across  the  symmetry  plane,  as  well  as  by  direct 
photon  coupling  within  the  same  half  of  the  medium.  In  each  case,  the  cells  are  spaced 
logarithmically  in  optical  depth  (ItIt  ~  constant)  and  the  optical  depths  of  the  cells  close  to  the 
boundary  are  less  than  unity  to  allow  for  the  anticipated  rapid  change  of  the  source  function 
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log  r  (LINE  CENTER! 

Fig.  I.  The  steady-state  source  function  for  a  Doppler  profile  is  shown  for  an  optical  depth  of  2.82  x  10'  and 
a  Planck  function  8  =  1.  Results  obtained  from  the  present  treatment  are  shown  along  with  the  exact 
solutions  obtained  by  Avrett  and  Hummer  in  Ref.  I  The  optical  depth  is  measured  from  the  center  of  the 
plane-parallel  medium  to  the  edge  and  at  line  center.  Quenching  parameters  of  10  4  (effectively  thick)  and 
10"®  (effectively  thin)  are  assumed  in  the  curses  shown. 


lop  r  l LINE  CENTEBl 

Fig.  2.  As  in  Fig.  I.  except  that  optical  depths  of  28.2  and  2.82  x  |0'  are  shown  and  only  one  quenching 

parameter.  I0'4.  is  considered. 
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line  tranNported  b\  both  direct  and  probabilistic  methods; - .  lines  transported  using  the  probabilistic 

formalism  only. 


near  the  edge  of  the  medium.  We  have  also  computed  results  for  50  cells  per  half-slab,  which 
are  not  plotted  for  reasons  of  clarity. 

To  obtain  solutions  within  a  few  per  cent  of  the  exact  values,  only  25  cells  need  be  used  for 
optical  depths  of  -  10-’.  but  50-75  cells  are  necessary  at  r  -  10*.  These  considerations  apply  for 
both  the  effectively  thick  and  the  effectively  thin  cases.  As  is  seen  in  Fig.  2.  for  r  =  2.8  x  10r.  the 
25-celJ  calculation  yields  source  functions  a  factor  of  2.5  greater  than  the  exact  solution  at  small 
optical  depths,  but  goes  to  the  Planck  function  IS  =  1)  at  the  correct  optical  depth.  The  75-cell 
model  yields  source  functions  about  25%  too  high  in  the  nonthermalized.  effectively  thin  regions  at 
this  very  high  optical  depth. 

6.  APPLICATION  TO  AN  ALUMINUM  PLASMA 

The  method  described  here  appears  to  hold  considerable  promise  for  computationally 
efficient  modeling  of  laboratory  plasmas.  Detailed  multi-angle,  multi-frequency  calculations 
have  shown3  that  the  ionization  dynamics  and  energetics  of  high-temperature  plasmas  are 
significantly  affected  by  the  trapping  of  line  and  continuum  radiation.  Accurate  simulation  of 
such  plasmas  and  computation  of  diagnostically  useful  line  ratios  requires  that  all  optically 
thick  lines  be  transported.  Given  these  considerations,  we  have  applied  the  method  to 
determine  the  -date  of  an  optically  thick  aluminum  plasma  in  coilisional-radiative  equilibrium. 
Specifically,  we  have  chosen  a  planar  plasma  at  a  temperature  of  600 eV.  with  a  half-width  of 
750^  and  a  total  ion  density  of  lx  I01’  ions  per  cm3.  As  a  framework  within  which  the 
technique  is  applied,  we  have  used  the  multistage,  multilevel  aluminum  model,  which  has  been 
extensively  described  elsewhere.4  The  accompanying  energy  level  diagram  (Fig.  3)  indica’es  the 
excited  state  structure  and  transported  lines  of  the  Al  (XII)  and  Al  (XIII)  stages,  which  are 
dominant  at  our  chosen  temperature  and  density. +  Additionally,  the  Al  (XI)  ls;2s-2p.  -3 s,  -3 p. 
-Id.  -4 d,  and  \s2p:  excited  states  are  included,  with  transport  of  the  satellite  line  ls:2p  :P- 
\s2p:  :D.  This  model  has  been  extensively  applied  using  direct  methods  to  solve  the  radiative 
transport  equation  and  the  results  will  be  thoroughly  described  for  a  range  of  plasma  conditions 
in  an  upcoming  paper.  The  salient  features  of  the  results  for  this  particular  plasma  are  now 
described  in  the  separate  sections  below. 

"In  these  calculations,  photoionization  has  been  neglected,  as  hase  radiative  recombination  transitions  except  for  those 
onto  the  bare  proton  to  form  H  n  =  I  and  n  =  '.  onto  Hn  =  1  to  form  He  I  j:i5  and  I  s2p‘P.  and  onto  He  ls;iS  lo  form  Li 
Ij-VS 
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(a)  Comparison  of  Voigt  and  Doppler  profiles 

We  have  performed  a  separate  calculation  for  the  identical  plasma  directly  solving  the 
transport  equation  with  Voigt  profiles  for  each  line  whose  widths  are  determined  by  the  total 
lifetime  of  the  respective  upper  levels.  The  calculation  revealed  very  little  difference  in  state 
populations  from  that  observed  using  a  pure  Doppler  profile.  At  this  density  and  size,  therefore, 
a  Doppler  profile  is  an  adequate  representation  of  the  line  opacity  of  the  plasma.  At  higher  ion 
densities,  the  increased  wings  of  the  Voigt  profile  coupled  with  higher  optical  depths  causes  a 
significant  lowering  of  the  excited  state  densities  compared  to  that  obtained  with  a  Doppler 
profile.' 

(b)  Effect  of  photon  pumping 

The  effects  of  radiative  pumping  of  the  excited  states  on  the  overall  state  of  the  plasma  are 
perceived  most  readily  by  solving  for  the  plasma  state  ignoring  radiative  transfer,  that  is. 
assuming  that  lines  which  are  really  optically  thick  are  instead  thin.  This  is  done  in  our 
calculation  as  a  means  of  providing  an  initial  guess  as  a  basis  for  iteration  to  the  true  optically 
thick  state.  In  Figs.  4  and  5.  we  have  plotted  the  level  densities  as  a  function  of  c  for 
helium-like  (XII)  and  hydrogen-like  (XIII)  aluminum,  respectively.  Both  the  assumed  thin  and 
true  optically  thick  solutions  are  shown.  The  most  obvious  feature  of  Figs.  4  and  5  is  that  the 
excited  state  densities  are  raised  well  above  their  optically  thin  values  by  repeated  line  photon 
scatterings  (re-excitations  of  the  level).  Relative  to  the  ground  state  densities,  the  Al  (XII) 
ts2plP  population  at  plasma  center  is  increased  by  a  factor  of  170  and  that  of  the  Al  (XIII) 
n  =  2  state  is  raised  by  a  factor  of  68.  These  values  are  in  reasonable  agreement  with  the 


Fig.  4.  Excited  and  ground  state  densities  of  Al  (XIII  are  shown  as  functions  of  position  in  the  plasma. 
Results  are  given  for  both  thin  i  ignoring  radiative  transferl  and  thick  calculations.  Also  plotted  are  densities 
obtained  when  transport  of  excited  state  lines  is  neglected,  when  such  densities  differ  appreciably  from 

those  of  the  full  calculation. 
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enhancement  of  excited  state  populations  predicted  by  photon-escape  theory  arguments.”  The 
table  of  line  power  outputs  reveals  that,  at  our  chosen  density  of  1  x  io1*,  the  Al  < XII ) 
ls:-lr2p'P  line,  for  instance,  is  effectively  thin.  That  is.  even  though  the  line  center  optical 
depth  happens  to  be  50  in  this  case,  the  emergent  power  is  attenuated  only  slightly  from  the 
optically  thin  value.  Most  of  the  line  photons,  therefore,  escape  the  plasma  after  many 
scatterings.  For  the  effectively  thin  case.  ,V„  the  number  of  scatterings,  is  equal  to  P/'."  where 
P ,  is  the  mean  photon  escape  probability  for  the  plasma.  To  produce  170  scatterings  for  this 
line  requires  a  P,  of  about  1/170  or.  as  revealed  by  Eq.  (3)  and  the  subsequent  approximate 
expressions,  an  optical  depth  of  25.  This  optical  depth  is  exactly  halfway  between  the  surface 
and  center  of  the  plasma,  confirming  the  usefulness  of  the  mean  escape  probability  concept. 
Similarly,  good  agreement  is  seen  for  the  1—2  line  of  Al  (XIII). 

Figures  4  and  5  also  illustrate  spatial  gradients  in  certain  of  the  excited  states  of  this  plasma, 
which  is  of  course  of  uniform  temperature  and  total  ion  density.  The  gradients  appear  in  the 
excited  states  which  are  controlled  primarily  by  radiative  pumping.  This  pumping  occurs  either 
directly  to  the  excited  state  in  question  or  to  an  adjoining  state  which  is  strongly  collisionully 
coupled  to  the  state.  The  \s2p3P  state  of  Al  (XII)  is  so  strongly  controlled  by  collisions  that  no 
significant  spatial  gradient  appears.  The  gradients  in  the  radiatively-controlled  excited  states  are 
a  result  of  the  lack  of  pumping  radiation  incident  on  the  outer  edge  of  the  plasma.  This  results 
in  a  lowered  excited  state  density  at  the  outer  boundary,  hence  a  lowered  inward  line  photon 
flux  emitted  near  the  outer  boundary.  The  interior  of  the  plasma  (tp  1)  can  thus  sense  the 
boundary  of  the  plasma  and  the  result  is  the  spatial  profiles  of  Figs.  4  and  5.  These  effects  have 
been  discussed  in  detail  elsewhere.'8 

(c)  Effects  of  excited  state-to-excited  state  photo-pumping. 

Also  indicated  in  the  power  output  table  and  plotted  in  Figs.  4  and  5  are  the  effects  of 
photon-pumping  by  lines  connecting  excited  states.  This  is  of  particular  interest  since  the 
transport  of  such  lines  is  usually  neglected.  The  optical  depths  of  such  lines  are  strongly 
dependent  on  photon  pumping  of  their  lower  levels  from  the  ground  state.  Thus,  they  are  not 
easily  estimable  and  may  vary  extremely  sensitively  with  temperature,  density  and  plasma  size. 
We  plan  a  detailed  study  of  these  effects,  which  will  be  subsequently  reported.  We  confine  our 
remarks  here  to  the  specific  plasma  discussed  above.  In  the  particular  plasma  under  con¬ 
sideration.  the  Al  (XII)  !j2p'P-ls3</'D  and  Al  (XIII)  n  =  2  —  n  =3  lines  are  the  only  excited 
state  lines  with  appreciable  optical  depths  (0.61  and  0.38  at  line  center,  respectively).  The 
principal  effect  of  neglecting  transport  of  these  lines  is  to  underestimate,  by  about  a  factor  of  2. 
the  populations  of  the  Al  (XIII)  n  =  3  and  Al  (XII)  ls3d'D  states.  The  n  =  2  level  populations, 
which  depend  primarily  on  pumping  from  the  ground  state,  are  not  significantly  affected.  Nor 
are  the  n  =  4  levels  substantially  increased,  since  lines  arising  from  n  =  2  and  n  =  3  pumping 
n  =4  are  optically  thin  in  both  species.  The  optical  depths  of  the  above  mentioned  2  —  3  lines 
imply  that  about  one  scattering  will  follow  collisionai  excitation  of  the  lines,  resulting  in  the 
doubling  of  the  upper  level  populations  when  the  transport  effect  is  accounted  for.  The  power 
output  of  the  Al  (XIII)  1-3  line  increases  by  43%.  rather  than  doubling,  because  near  the 
surface  of  the  plasma  the  n  =  3  level's  population  only  increases  by  25%.  since  the  pumping 
effects  are  not  as  large  at  the  boundary.  The  level  populations  near  the  plasma  edge  are.  of 
course,  more  important  than  those  in  the  interior  in  determining  power  output.  The  calculated 
outputs  of  the  2-3  lines  in  both  Al  (XII)  and  Al  (XIII)  decrease  when  the  lines  are  transported 
because  of  reabsorption  and  quenching.  The  optically  thick  calculation,  however,  yields  a  vastly 
greater  output  in  these  2-3  lines  than  that  obtained  when  all  lines  (including  resonance  lines)  are 
assumed  thin.  This  is.  of  course,  due  to  resonance  pumping  of  the  n  =  3  level.  It  is  thus  seen 
that  two-level-atom  scattering  theory6  gives  good  results  for  1—2  and  1  —4  lines  in  this  specific 
case  but  that  the  1—3  line  in  Al  (XIII)  is  significantly  altered  in  output  by  radiative  pumping 
from  n  =  2. 

Id)  Advantages  and  limitations  of  the  formalism 

No  significant  differences  were  observed  in  the  results  obtained  using  conventional  multi¬ 
frequency,  multiangle  transport  compared  to  those  yielded  by  the  present  method.  Since  our 
formalism  eliminates  frequency-  and  angle-integration,  it  is  quite  advantageous  numerically.  A 
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further  advantage  of  the  coupling  constant  technique  wav  achieved  specifically  in  its  use  in  a 
multistage,  multilevel  model.  As  the  excited  state  densities  build  up  during  the  "time-stepping” 
iteration  scheme.1'  the  ground  states  change  very  little  over  5-10  steps.  Thus,  the  coupling 
constants,  which  depend  on  the  ground  state  densities  only  for  the  thick  resonance  lines,  need 
be  updated  for  these  tines  (the  vast  majority  of  lines  transported:  onlv  every  5-10  steps 

Two  major  limitations  affect  the  application  of  the  present  scheme.  First,  efficient  numerical 
algorithms  for  the  cell-to-cell  coupling  constants  have  been  obtained  only  for  the  case  of  a 
spatially  constant  line  Doppler  width.  It  is  far  from  obvious  that  similar  algorithms  can  be 
obtained  to  apply  to  the  general  case,  where  gradients  of  varying  magnitudes  and  directions  are 
conceivable.  Second,  the  emitted  radiation  cannot  be  resolved  in  angle.  This  second  limitation 
can  be  mitigated  for  the  steady  state  case  by  using  a  detailed  multi-ray  procedure  to  compute 
the  angular  distribution  of  radiation  once  the  excited  state  densities  are  determined  by  the  more 
efficient  coupling  method.  However,  for  a  time  dependent  calculation  calling  a  rav-tracing 
scheme  every  timestep  to  obtain  angular  resolution  would  clearly  be  as  expensive  as  using  the 
multi-frequency,  multi-angle  technique  to  directly  compute  the  state  of  the  plasma. 

'  SUMMARY  AND  CONCLUSIONS 

An  efficient  and  readily  adaptable  technique  for  radiative  transfer  calculations  has  been 
developed  which  eliminates  the  need  for  angle-and-frequency  integrations.  A  specific  numerical 
algorithm  has  been  developed  for  a  Doppler  profile  in  plane-parallel  geometry  and  applied  in  a 
detailed  analysis  of  multilevel  coupling  in  a  highly  ionized  aluminum  plasma.  Extension  of  the 
method  to  other  geometries  and  line  profiles  te  g.  Voigt  profiles)  can  be  achieved  by  obtaining 
expressions  for  cell-to-cell  photon  escape  and  absorption  probabilities  for  these  cases.  The 
numerical  consequences  of  this  additional  physics  will  appear  solely  in  the  use  of  different 
coupling  coefficients.  Use  of  the  technique  has  made  possible  calculations  which  have  revealed 
in  detail  the  effects  of  excited-state-line  transport  on  line  power  outputs  and  level  populations  in 
an  aluminum  plasma.  In  addition,  the  uses  and  limitations  of  an  even  simpler  one-cell  scattering 
theory  approach  have  been  assessed. 


Sole  added  in  proof  :  The  20-point  cubic  spline  lit  10  P, <  r„i  for  r„  <  3  is  extremely  accurate  for  both  P,  and  ns  gradient 
However,  a  somewhat  less  time-consuming  approximation  to  P, l  r„)  for  r„  s  J  is  given  by 

■a,.,. _ ! _ 

'  I  -  1. 86l607r„ -U.8l7.Wv 

This  expression  differs  from  the  exact  P,  by  a  maximum  of  "r>  for  Os  r„s  J.  The  gradient  of  P,.  however,  is  not  well 
represented  for  r„  <  0.5.  Despite  this  fact,  use  of  the  above  simplified  expression  yields  results  which  are  everywhere  within  a 
few  per  cent  of  the  20-point  spline  results  for  the  Avretl-Hummer  cases  w  hich  are  illustrated  in  Figs  I  and  2.  Both  the  above 
expression  and  the  one  given  for  r  >  3  in  Section  3  are  analytically  integrablefor  use  in  Eq.  i4l.  and  are  equal  at  the  joining  point 
r„  =  3  to  several  significant  figures. 


\cknowiedyemenl— This  work  was  supported  in  part  by  ihe  Defense  Nuclear  Agency. 


REFERENCES 

1.  E.  H  Avrett  and  D.  G.  Hummer.  Vf.VR.A.S  130.  29?  1 19651 

2.  T  Holstein.  Phvr  Ret  72.  1212(1947). 

3.  J  Davis.  K.  G.  Whitney,  and  J.  P.  Apruzese.  JQSRT  20.  353  1 1978). 

4.  J.  Davis,  and  K.  G.  Whitney.  J.  Appl.  Phys.  47.  1426  1 19761. 

5  J  P  Apruzese.  J.  Davis,  and  K.  G  Whitney  .  JQSRT  17.  557  1 1977), 

6.  J  P  Apruzese.  J.  Davis,  and  K.  G.  Whitney.  J.  Ph\s.  B  II.  L643  1 1978). 

*  J.  P  Apruzese.  J.  Davis,  and  K.  G.  Whitney.  J.  Appl.  Ph\s.  47.  4433  ( I9’6i. 

8.  J  P  Apruzese.  J.  Davis,  and  K.  G.  Whitney.  J.  A ppl.  Phys.  48.  66"  1 197”). 

9.  K.  G  Whitney,  J.  Davis,  and  J.  P  Apruzese.  Cooperative  Effects  in  Matter  and  Radiation  i  Edited  by  C.  M.  Bowden.  D 
W.  Howgate.  and  H.  R.  Robl).  p  115.  Plenum  Press.  New  York  i  I97"i 


OSRT  Vol  ;j  Vo  '—i • 


E-9 


